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CHAPTER ONE 



INTRODUCTION 



Lake Huron was the subject of an intense limnological survey during the 
spring, summer, and fall of 1980. This binational study was undertaken by the 
United States and Canada to evaluate Lake Huron water quality throughout the 
ice-free seasons. This study was part of the Great Lakes surveillance program 
which included monitoring the five Great Lakes on a rotational basis. Lake 
Huron was last monitored in 1974-75. The surveillance program had four major 
objectives: 

1. To determine the status of the open and nearshore waters of the Great 
Lakes and to compare the status with the standards, criteria, and 
objectives for the protection of aquatic life in the Great Lakes. 

2. To provide data to characterize the chemical, physical, microbiologi- 
cal, and biological aspects of the environment against which future 
changes may be evaluated. 

3. To compare present data with data collected in the past in order 
to determine if the Great Lakes are changing and how these changes 
may be occurring. 

4. To determine how these changes are related to waste reduction and 
pollution abatement programs. 

Two of these objectives are addressed at length in this report. These two 
objectives are: Comparison of current limnological conditions to historical 
conditions to identify long-term trends (objective 3); and identification of the 
extent of oligotrophic and eutrophicatic waters to locate areas impacted by 
pollutants (objective 4). Objective 4 is particularly important as a means of 

1 



evaluating the improvement mitigating measures have had on Lake Huron water 
quality. This report deals with that aspect of the Lake Huron surveillance 
program concerning nutrient chemistry and physical-chemical variables. Included 
in this material is the distribution of chlorophyll biomass and the relationship 
of chlorophyll distribution to physical-chemical variables. 

The Federal Water Pollution Control Act as amended in 1972 provides 
authority to the United States Environmental Protection Agency (USEPA) to 
conduct surveillance and monitoring of the Great Lakes. In practice. Great 
Lakes surveillance programs are binational efforts between the United States and 
Canada with the Great Lakes National Program Office of USEPA leading the effort 
for the United States. The authority for the Lake Huron study was from the same 
source as for the Lake Michigan intensive survey of 1976-77. This authority is 
detailed in the Lake Michigan intensive survey report (Rockwell £t £l. 1980). 

The 1980 Lake Huron survey was conducted through several United States and 
Canadian research organizations. The results presented below were compiled by 
three organizations. The Great Lakes National Program Office (GLNPO) of USEPA 
and the Canada Centre for Inland Waters (CCIW) of the Canadian Department of the 
Environment conducted the open-lake sample collection and analysis. The Great 
Lakes Research Division (GLRD) of The University of Michigan assisted the USEPA 
in data analysis and report preparation. The results of the US-Canadian 
sampling and data analysis presented below are not the complete results from the 
1980 surveillance program. Other results from the Lake Huron survey include: 
nearshore monitoring conducted by the Michigan Department of Natural Resources 
and the Ontario Ministry of the Environment , phytoplankton species enumeration 
and identification conducted by GLRD, GLNPO, and CCIW, zooplankton species 
enumeration and identification conducted by GLRD, trace metal analysis conducted 



by GLRD, and an intensive survey of Saginaw Bay conducted by USEPA Large Lakes 
Research Station at Grosse lie, Michigan. 

The importance of the Lake Huron 1980 study is enhanced when the inter- 
national aspects of the lake are considered. As an international boundary. 
Lake Huron is jointly managed by the United States and Canadian governments. 
This management is implemented with guidelines set by the International Joint 
Commission (I.J.C.). The I.J.C. requires cooperative studies such as the Lake 
Huron survey to develop realistic recommendations for lake management as well 
as to provide information on the status of the Great Lakes. Thus, the surveil- 
lance program fills a dual role as an evaluation of the status of Lake Huron in 
1980 and as a database for future management strategies. An additional aspect 
of the binational position of Lake Huron is that pollutants entering the system 
are likely to cross an international boundary. Whole lake surveillance provides 
one means of estimating this ndesirable international exchange. 

The design of the 1980 study allowed broad inferences of overall water 
quality in Lake Huron. But, the design has its limitations; and as a conse- 
quence, the inferences from this study are also limited. The most important 
limitation is that the 1980 study is only a one-year surveillance program. 
Observed seasonal trends from 1980 cannot be extrapolated to other years. 
Similarly, long-term trend analysis in water quality is not possible from a one- 
year study. This shortcoming has been somewhat alleviated by comparison of the 
1980 data to a large historical database from previous Lake Huron studies. 
The reader is referred to the report by Rossmann and Treese (1981) for a 
complete bibliography of Lake Huron studies. 



A further limitation of this study is the areal coverage provided by the 
station grid. Lake Buron, Georgian Bay, and the North Channel \vere sampled with 
a fairly coarse station grid. This study did not include Saginaw Bay, which is 
discussed in Smith et al. (1983). But, discharges from Saginaw Bay are a major 
factor affecting the water quality of southern Lake Huron. This report, while 
not considering the water quality in Saginaw Bay itself, does include the impact 
of Saginaw Bay water mixing into open Lake Huron. 



CHAPTER TWO 
CONCLUSIONS 

Water quality in Lake Huron was very good throughout most of the lake 
during 1980. This conclusion was based on an analysis of the spatial and tem- 
poral distribution of nutrients, chlorophyll, and particulate materials. 
Nutrient concentrations in Lake Huron were low compared to the other Great 
Lakes, and the annual ranges of nutrient concentrations were small; these two 
aspects of the limnology of Lake Huron indicated an oligotrophic ecosystem. 

Although the water quality in Lake Huron was characterized as very good, 
local, nearshore nutrient loadings produced occasional regions of poor water 
quality. These regions, which were primarily identified by high nutrient and/or 
chlorophyll concentrations, were both local in areal extent and short lived in 
duration. Regions of poor water quality in Lake Huron were not chronic, but 
rather occasional during periods of high river runoff, high algal productivity, 
etc. Some of these regions of poor water quality were: the Ontario shore of 
southern Lake Huron, the mouth of Saginaw Bay, and Thunder Bay. Water quality 
in the mouth of Saginaw Bay has improved since the last Lake Huron intensive 
survey. Remedial steps taken between 1974 and 1980 significantly reduced 
nutrient loads to lower Saginaw Bay (Dolan et al. 1981). As a result, poor 
water quality in outer Saginaw Bay and adjacent Lake Huron was only occasionally 
observed during 1980. The effects of remedial measures in other locations of 
Lake Huron were not evident in 1980. 

In addition to regions of high nutrient loadings, some regions received 
river runoff which was extremely low in nutrients and dissolved salts. 
In particular, rivers which drained the Canadian Shield into the North Channel 



and northern Georgian Bay had very low chloride and nutrient loads. The effect 
of these rivers was prevalent throughout much of northern Lake Huron during the 
spring thaw. 

For nearshore areas, 1980 Secchi disk depths, total nitrogen, dissolved 
reactive silica, sulfate, and conductivity were greater than those reported 
for 1974. Conversely, the 1980 total phosphorus and chloride concentrations 
were less than those reported for 1974. 

Within the open lake, 1980 total phosphorus and dissolved reactive silica 
were less and dissolved nitrate plus nitrite was greater than those reported 
for 1974. The decrease in total phosphorus followed a series of relatively 
colder winters. There was a 50% decrease in mean dissolved reactive silica 
and dissolved nitrate plus nitrite maximum depletions in 1980 relative to those 
of 1974. Thus Lake Huron's trophic status appears to have improved between 1974 
and 1980. 

A result of the 1974 and 1980 intensive studies of Lake Huron was the 
recognition that distinct zones or segments of differing water quality existed. 
After the 1974 study. Lake Huron was divided into 19 segments of somewhat homo- 
geneous water quality in each segment (IJC 1976). The boundaries between these 
segments did not change among cruises or seasons. An alternative approach to 
segmentation was used with the 1980 data. Boundaries were based on analytical 
mathematical techniques and allowed to change among cruises. The results of the 
1980 analysis showed that 19 distinct segments did not exist in Lake Huron. 
Rather, only six to nine homogeneous water masses were identified in the lake 
during each cruise, implying that fewer stations than those used in 1980 may 
adequately characterize Lake Huron. Furthermore, the distribution of those 
water masses varied considerably among cruises. This variation was linked 



to changes in the physical, including thermal, regime throughout the year. 
The overall inference from this analysis was that the physical environment in 
Lake Huron controlled the location of water masses, while biological events 
controlled the concentration of nutrients within each water mass. Because most 
biological processes were not measured during 1980, the inference that biologi- 
cal events controlled nutrient concentrations was somewhat speculative. 

The inferences from the 1980 study were expanded by the inclusion of his- 
torical data from Lake Huron. The historical data included nutrient and physi- 
cal measurements since 1954. These data were analyzed along with the 1980 
results to determine if long-term trends of water quality exist in Lake Huron. 
The database was limited in the number of variables available for analysis and 
the annual coverage. Nonetheless, seven variables (conductivity, chloride, 
sulfate, soluble silica, nitrate plus nitrite, soluble reactive phosphorus, 
and total phosphorus) were analyzed for long-term trends ranging from 9 to 28 
years. The annual range of historical data limited the analyses to spring and 
summer months only. 

Long-term changes in the chemistry of Lake Huron waters have been identi- 
fied using the results of twelve studies spanning a 26-year period of observa- 
tions. As might be expected of a dynamic system lake, Lake Huron long-term 
trends can be complex. For roughly one-half of the parameters investigated, 
changes over the 26 years were curvilinear or oscillatory. Of these, one-half 
had a linear trend superimposed over the oscillations. For each of the various 
ways we manipulated the data base, sulfate and nitrate were found to increase. 
For the majority of cases, dissolved silica and chloride decreased. Sulfate 
increased at a rate of 0.05 to 0.13 mg/L/yr or 0.3 to 0.8% of the 1980 mean con- 
centration per year. Nitrate increased at a rate of 0.0030 to 0.0089 mg/L/yr 



or 1.1 to 3.1% per year. Dissolved silica decreased at a rate of 0.040 to 
0.049 mg/L/yr or 2.8 to 3.6% of the 1980 mean per year. Chloride decreased at a 
rate of 0.02 to 0.06 mg/L/yr or 0.4 to 1.2% of the 1980 mean per year. All four 
parameters should continue to be monitored. 

For those parameters which showed reversals in trend over the period of 
observation, reversals consistently occurred in the mid-1970s. It is unfortu- 
nate that data were not available for this critical period of time. Since the 
mid-1970s, dissolved reactive silica and soluble reactive phosphorus have 
decreased, and nitrate and sulfate have increased. Changes in long-term trends 
could be or have been attributed to factors such as changes in sampling and 
analytical methods, reduction in anthropogenic inputs from the major point 
sources, and differences in biological utilization of these elements. Despite 
the changes which have occurred in Lake Huron between 1954 and 1980, Lake Huron 
is still an oligotrophic lake. 



CHAPTER THREE 

METHODS 

by David C. Rockwell 

One hundred thirty eight stations were sampled six times from 10 April 
through 8 November 1980 by the Great Lakes National Program Office (GLNPO) 
USEPA, Region V, and the Canada Centre for Inland Waters (CCIW). Two winter 
surveys of 11 stations each were conducted during 20-21 January and 24-25 Febru- 
ary 1981 (stations 1, 3, 5, 7, 9, 10, 13, 15, 16, 19, and 21). The station 
locations are given in Figure 3.1 and latitudes, longitudes in Table 3.1. 
Cruise dates for each major lake segment are given in Table 3.2. Water column 
depths for each station are given in Table 3.3. 

The sampling strategy conformed to the Great Lakes International Surveil- 
lance Plan (GLISP). In homothermous waters four depths were sampled: surface 
(1-m depth), mid-depth, bottom depth less 10 meters, and bottom depth less 

2 meters. In thermally stratified waters the mid-depth sample was replaced by 
three samples: lower epilimnion (2m above the upper knee), mid-thermocline, 
and upper hypolimnion (2 m below the lower knee). Where samples were less than 

3 meters apart, one sample was omitted. Thus a minimum of four and a maximum 
of seven depths were sampled at each station. 

The sampling platforms used were the R/V R. Simons on cruises one, two, 
and four; the CSS Limnos on cruises three, five, and six; and the US Coast Guard 
Cutter Bramble for the winter cruises. 

Figure 3.2 is a flow chart of sample processing on the USEPA R/V R. Simons . 
Samples were collected by means of 8-L Niskin® water sampling bottles. An inte- 
grating sampler made by CCIW was used to sample the upper 20 meters for chloro- 
phyll, particulate organic carbon, and particulate organic nitrogen. 
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Table 3.1. 1980 Lake Huron surveillance stations. 





STATION NO. 


CROSS REFERENCE 


LATITUDE 


LONGITUDE 


1 






43 = 


'05 '24" 


82°23'30" 


2 






43 


11 24 


82 17 54 


3 






43 


15 25 


82 02 18 


4 






43 


19 30 


81 47 18 


5 






43 


32 54 


81 44 42 


6 






43 


28 00 


82 00 00 


90 






43 


24 00 


82 18 00 


7 






43 


20 30 


82 30 24 


8 






43 


34 00 


82 29 06 


9 






43 


38 00 


82 13 00 


91 






43 


42 00 


82 01 00 


10 


Master 


43 


45 12 


81 46 54 


11 






43 


57 24 


81 47 12 


12 






43 


53 24 


82 03 24 


92 






43 


48 30 


82 22 00 


13 






43 


45 12 


82 34 06 


14 






43 


56 30 


82 40 00 


15 






44 


00 00 


82 21 00 


16 


Master 


44 


07 54 


82 45 00 


17 


HU 


0050 


44 


06 00 


82 52 00 


18 


SB 


0051 


44 


07 25 


83 10 15 


19 


HY 


0049 


44 


09 00 


82 58 00 


20 


HY 


0048 


44 


13 00 


83 05 00 


21 


HU 


0047 Master 


44 


16 00 


83 12 00 


22 


SB 


0049 


44 


12 40 


83 22 40 


23 


HY 


0046 


44 


20 00 


83 18 00 


24 


HU 


0024 


44 


16 00 


82 55 00 


25 


HU 


0026 


44 


23 00 


83 16 00 


26 


HU 


0025 Master 


44 


20 00 


83 05 00 


27 






44 


11 54 


82 30 12 


28 






44 


12 18 


81 40 36 


93 






44 


06 00 


82 07 00 


29 






44 


22 00 


81 50 00 


30 






44 


28 00 


81 27 12 


31 






44 


51 00 


81 36 00 


32 






44 


27 12 


82 20 30 


33 






44 


30 00 


82 50 00 


34 






44 


38 24 


83 13 54 


35 






44 


51 00 


83 15 42 


36 






45 


02 06 


83 22 42 


37 






44 


45 42 


82 47 00 


38 






44 


44 24 


82 03 36 


39 






44 


39 24 


81 22 42 


40 






44 


53 54 


81 26 12 
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Table 3.1. Continued. 



STATION NO. CROSS REFERENCE LATITUDE LONGITUDE 



41 

A2 

43 

44 

45 

94 

46 

47 Master 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 



45' 


"05' 


'00' 


45 


13 


18 


45 


GO 


48 


45 


01 


00 


45 


08 


12 


44 


04 


15 


45 


04 


48 


45 


15 


18 


45 


16 


42 


45 


24 


48 


45 


32 


06 


45 


32 


00 


45 


39 


06 


45 


27 


00 


45 


31 


00 


45 


23 


30 


45 


31 


00 


45 


40 


00 


45 


52 


06 


45 


46 


00 


45 


54 


06 


45 


45 


CO 


45 


40 


30 


45 


42 


12 


45 


48 


48 


45 


50 


42 


45 


51 


48 


45 


56 


06 



81°32' 


18" 


81 49 


12 


82 00 


30 


82 41 


06 


82 59 


00 


83 05 


00 


83 14 


00 


83 20 


48 


82 27 


06 


81 55 


06 


82 02 


42 


82 16 


48 


82 38 


54 


82 54 


54 


83 25 


00 


83 39 


06 


84 05 


00 


83 43 


36 


83 16 


00 


83 01 


42 


83 31 


06 


83 55 


00 


84 11 


12 


84 50 


42 


84 45 


18 


84 34 


00 


84 17 


42 


83 54 


00 


(Cont: 


inued ) 
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Table 3.1. Continued. 



STATION NO. 



CROSS REFERENCE 



LATITUDE 



LONGITUDE 



1980 GEORGIAN BAY SURVEILLANCE STATIONS 



101 


44' 


°43'03 


102 


44 


48 30 


103 


44 


43 30 


104 Master 


44 


38 45 


105 


44 


47 48 


106 


44 


44 12 


107 


44 


53 20 


108 


44 


57 10 


109 


44 


52 18 


110 


45 


03 45 


HI 


45 


55 15 


112 


44 


55 12 


113 


45 


01 36 


114 


45 


08 20 


115 


45 


10 00 


116 


45 


21 13 


117 


45 


14 42 


118 


45 


09 10 


119 


45 


04 00 


120 


45 


13 00 


121 


45 


21 54 


122 


45 


28 50 


123 


45 


33 35 


124 


45 


40 44 


125 


45 


46 40 


126 


45 


50 00 


127 


45 


52 00 


128 


45 


42 12 


129 


45 


35 00 


130 Master 


45 


32 30 


131 


45 


14 18 


132 


45 


16 12 


133 


45 


22 13 


134 


45 


27 10 


135 


45 


31 39 


136 


45 


42 30 


137 


45 


43 00 


138 


45 


53 00 


139 


45 


52 24 


140 


45 


51 52 


141 


45 


56 00 


142 


45 


54 46 


143 


45 


49 52 


144 


45 


58 20 



13 



80' 


°51 


'24" 


80 


52 


18 


80 


37 


00 


80 


10 


00 


80 


14 


36 


80 


26 


06 


80 


17 


50 


80 


08 


06 


79 


58 


05 


80 


11 


28 


80 


36 


21 


80 


52 


30 


80 


52 


36 


80 


31 


24 


80 


17 


48 


80 


29 


12 


80 


52 


30 


81 


04 


03 


81 


15 


14 


81 


13 


36 


81 


11 


24 


80 


50 


15 


80 


36 


38 


80 


50 


20 


80 


45 


15 


80 


54 


00 


81 


00 


00 


81 


05 


24 


81 


05 


00 


81 


22 


00 


81 


26 


24 


81 


35 


00 


81 


35 


06 


81 


43 


46 


81 


40 


10 


81 


37 


12 


81 


22 


30 


81 


06 


30 


81 


15 


30 


81 


32 


08 


81 


31 


04 


81 


35 


42 


81 


47 


19 


81 


41 


55 


(Continued) 



Table 3.1. Concluded. 



STATION NO. 



CROSS REFERENCE 



LATITUDE 



LONGITUDE 



1980 NORTH CHANNEL SURVEILLAJNCE STATIONS 



68 




46' 


^02 '30 


69 




46 


04 42 


70 




46 


08 12 


71 


Master 


46 


14 00 


72 




46 


13 36 


73 




46 


11 12 


74 




46 


08 54 


75 




46 


05 00 


76 




46 


00 00 


77 




45 


58 12 


78 




46 


02 06 


79 




46 


07 24 


80 




46 


00 00 


81 




46 


04 42 


82 




45 


56 18 


83 




46 


00 00 


84 


Master 


46 


05 30 


85 




46 


06 00 


86 




46 


00 18 


87 




46 


03 40 


88 




46 


03 20 


89 




45 


55 00 



83' 


'51' 


'12" 


84 


01 


42 


83 


40 


18 


83 


44 


48 


83 


35 


24 


83 


21 


18 


83 


12 


04 


83 


25 


00 


83 


26 


00 


83 


11 


54 


83 


00 


CO 


82 


53 


09 


82 


51 


21 


82 


44 


36 


82 


45 


30 


82 


33 


00 


82 


33 


24 


82 


25 


30 


82 


23 


18 


82 


11 


50 


82 


00 


00 


82 


09 


40 
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Table 3.2. Cruise dates for Lake Huron study. 



Cruise Number 



Dates 1980 



Average Julian 
Day 



Cruise One 

Lake Huron 

Southern Basin 
Northern Basin 

North Channel 

Georgian Bay 



Cruise Two 

Lake Huron 

Southern Basin 
Northern Basin 

North Channel 

Georgian Bay 



Cruise Three 

Lake Huron 

Southern Basin 
Northern Basin 

North Channel 

Georgian Bay 



Cruise Four 

Lake Huron 

Southern Basin 
Northern Basin 

North Channel 

Georgian Bay 



13-15 April 
15-20 & 23 April 
20-21 April 
24-27 April 



10-11 May 

11-16 May 

16-17 May 

18-21 May 



Cruise Five 

Lake Huron 

Southern Basin 
Northern Basin 

North Channel 

Georgian Bay 



18-19 July 
20-25 July 
25-27 July 
27-30 July 



8-11 Sept. 

11-15 Sept. 

15-16 Sept. 

16-21 Sept. 



104 
108 
110 
116 



130 
133 
136 
140 



28-29 May 
29 May-4 June 
4-5 June 
5-7 June 



148 
152 
156 
157 



200 
203 
207 
210 



252 
256 
258 
261 



(Continued) 
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Table 3,2. Concluded. 



Average Julian 
Cruise Number Dates 1980 Day 



Cruise Six 

Lake Huron 

Southern Basin 22-27 Oct. 298 

Northern Basin 27-31 Oct. 302 

North Channel 31 Oct.-l Nov. 304 

Georgian Bay 2-4 Nov. 307 



Dates 1981 



Winter Cruise One 
Lake Huron 

Southern Basin 20-21 Jan. 20 



Winter Cruise Two 
Lake Huron 

Southern Basin 24-25 Feb. 56 
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Table 3,3. Depth in meters of permanent stations sampled in Lake Huron 1980 
surveillance study. 





Station 


Cruise 


Cruise 


Cruise 


Cruise 


Cruise 


Cruise 


Number 


One 


Two 


Three 


Four 


Five 


Six 


01 


11.5 


11.5 


10.5 


12.0 


11.0 


11.0 


02 


23.5 


23.75 


23.0 


24.0 


24.0 


24.0 


03 


19.0 


15.5 


14.0 


16.0 


14.0 


13.0 


04 


11.0 


10.25 


10.0 


11.0 


9.5 


10.0 


05 


12.5 


12.5 


12.0 


13.0 


11.0 


11.5 


06 


51.0 


52.0 


51.5 


53.0 


52.0 


51.5 


07 


10.0 


9.5 


11.0 


11.0 


12.0 


12.5 


08 


25.5 


27.5 


30.0 


27.0 


30.5 


30.0 


09 


54.0 


60.0 


59.0 


63.0 


61.0 


59.0 


10 


12.0 


11.5 


12.0 


11.0 


12.5 


13.0 


11 


10.0 


11.0 


10.5 


8.0 


11.0 


10.5 


12 


84.0 


88.0 


87.0 


85.0 


88.0 


87.5 


13 


11.5 


15.0 


17,5 


25.0 


18.0 


17.0 


14 


12.2 


14.0 


13.5 


21.0 


14.0 


13.0 


15 


64.5 


67.0 


66.5 


68.0 


68.0 


67.0 


16 


46.0 


54.0 


47.0 


51.0 


48.0 


48.0 


17 


15.0 


25.0 


15.0 


17.0 


15.0 


15.5 


18 


26. 


29.5 


28.5 


32.0 


29.0 


29.0 


19 


31.0 


41.0 


35.0 


35.0 


36.0 


35.0 


20 


45.0 


51.0 


46.0 


51.0 


46.0 


46.0 


21 


40.0 


45.0 


41.0 


44.0 


42.0 


41.0 


22 


18.0 


20.0 


20.5 


22.0 


21.0 


21.0 


23 


13.0 


14.75 


10.5 


11.0 


12.5 


12.0 


24 


64.0 


65.5 


62.0 


63.0 


63.0 


62.0 


25 


10.5 


11.75 


12.0 


14.0 


14.0 


13.0 


26 


57.0 


45.0 


58.0 


52.0 


57.0 


58.0 


27 


61.0 


57.0 


56.0 


53.0 


56.0 


55.0 


28 


— 


28.0 


17.0 


19.0 


18.5 


19.0 


29 


123.0 


124.0 


119.0 


113.0 


125.0 


122.0 


30 


30.0 


20.0 


20.0 


22.0 


23.5 


21.0 


31 


— 


89.5 


86.0 


100.0 


88.0 


86.0 


32 


79.0 


80.0 


80.0 


82.0 


81.0 


80.0 


33 


— 


72.0 


55.0 


73.0 


73.0 


73.0 


34 


11.0 


10.0 


11.0 


12.0 


10.5 


10.0 


35 


15.0 


15.0 


18.0 


16.0 


20.5 


16.0 


36 


13.0 


10.5 


10.0 


9.0 


11.0 


10.5 


37 


76.0 


76.0 


74.0 


78-0 


75.0 


72.0 


38 


— 


138.0 


134.0 


144.0 


136.0 


134.0 


39 


28.0 


30.0 


27.2 


33.0 


29.0 


30.0 


40 


~ 


21.5 


20.0 


29.0 


22.5 


20.5 



(Continued) 
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Table 3.3. Continued. 





station 


Cruise 


Cruise 


Cruise 


Cruise 


Cruise 


Cruise 


Number 


One 


Two 


Three 


Four 


Five 


Six 


41 


17.0 


29.5 


20.0 


22.0 


21.0 


23.0 


42 


50.0 


34.0 


42.5 


45.0 


38.0 


40.0 


43 


144.0 


170.0 


176.0 


172.0 


182.0 


182.0 


44 


161.0 


166.0 


164.0 


166.0 


164.0 


164.0 


45 


74.0 


74.0 


101.0 


78.0 


102.0 


111.0 


46 


12.0 


14.0 


13.0 


15.0 


13.0 


13.0 


47 


15.0 


21.5 


23.5 


12.0 


20.0 


23.0 


48 


109.0 


104.0 


111.0 


110.0 


111.0 


110.0 


49 


45.0 


47.5 


49.0 


55.0 


45.5 


48.0 


50 


30.0 


44.0 


29.0 


31.0 


29.5 


29.0 


51 


63.0 


73.0 


63.0 


65.0 


65.0 


64.0 


52 


64.0 


75.0 


65.0 


65.0 


67.0 


66.0 


53 


93.0 


83.0 


90.0 


86.0 


92.0 


91.0 


54 


129.0 


122.0 


118.0 


125.0 


110.0 


115.0 


55 


15.0 


23.0 


15.0 


10.0 


21.0 


22.5 


56 


19.0 


18.5 


24.0 


27.0 


22.5 


22.0 


57 


108.0 


122.0 


132.0 


122.0 


126.0 


130.0 


58 


33.0 


56.5 


13.0 


38.0 


15.5 


56.5 


59 


22.0 


36.0 


18.0 . 


24.0 


20.5 


19.0 


60 


31.0 


34.0 


30.0 


28.0 


24.5 


38.5 


61 


112.0 


117.0 


118.0 


116.0 


117.0 


116.0 


62 


45.0 


46.0 


46.0 


47.0 


45.0 


45. C 


63 


13.0 


13.0 


16.5 


14.0 


16.0 


16.0 


64 


15.0 


11.0 


55.5 


39.0 


48.0 


52.0 


65 


48.0 


50.0 


59.0 


51.0 


61.0 


60.0 


66 


66.0 


70.0 


70.0 


71.0 


70.0 


69.0 


67 


43.0 


38.0 


29.0 


32.0 


33.5 


33.5 


68 


10.0 


11.5 


17.0 


18.0 


16.0 


17.0 


69 


15.0 


7.5 


11.7 


15.0 


10.5 


11.0 


70 


20.0 


22.0 


22.5 


23.0 


22.5 


22.0 


71 


33.0 


32.0 


36.0 


35.0 


36.0 


34.0 


72 


37.0 


39.0 


23.5 


27.0 


25.5 


28.0 


73 


23.0 


22.0 


19.5 


19.0 


19.0 


20.0 


74 


13.0 


14.0 


13.0 


14.0 


14.5 


13.0 


75 


49.5 


45.0 


45.7 


47.0 


50.0 


50.0 


76 


52.0 


56.5 


58.0 


59.0 


60.0 


60.0 


77 


35.0 


77.0 


78.0 


75.0 


78.0 


79.0 


78 


— 


46.0 


45.0 


47.0 


47.0 


46.0 


79 


— 


25.0 


28.0 


28.0 


25.0 


27.0 


80 


~— 


40.0 


40.0 


41.0 


41.5 


39.5 



(Continued) 



18 



Table 3.3. Continued. 





Station 


Cruise 


Cruise 


Cruise 


Cruise 


Cruise 


Cruise 


Number 


One 


Two 


Three 


Four 


Five 


Six 


81 


«.— 


32.5 


32.5 


34.0 


33.0 


33.0 


82 


— 


25.0 


30.0 


29.0 


30.0 


25.5 


83 


— 


32.5 


31.0 


31.0 


31.0 


32.0 


84 


— 


36.0 


36.0 


37.0 


37.0 


36.0 


85 


— 


20.0 


29.5 


30.0 


32.0 


27.0 


86 


— 


30.0 


29.5 


18.0 


30.0 


30.0 


87 


— 


31.5 


46.0 


36.0 


48.0 


33.0 


88 


— 


31.0 


25.0 


33.0 


33.0 


32.5 


89 


— 


31.0 


40.0 


45.0 


42.0 


40.0 


90 


33.5 


42.5 


42.5 


43.0 


43.0 


42.5 


91 


71.0 


78.5 


78.5 


90.0 


82.0 


78.0 


92 


59.0 


62.5 


62.0 


65.0 


63.5 


62.5 


93 


83.0 


89.5 


88.5 


88.0 


89.0 


89.0 


94 


13.2 


20.0 


14.5 


20.0 


15.0 


14.5 


101 


90.0 


90.0 


90.0 


89.0 


91.0 


90.0 


102 


96.0 


95.0 


91.5 


95.0 


95. C 


95.0 


103 


31.0 


31.0 


34.0 


26.0 


33.0 


34.0 


104 


59.0 


58.0 


58.0 


59.0 


57.0 


57.0 


105 


53.0 


59.0 


60.0 


60.0 


61.0 


60.0 


106 


79.0 


89.5 


88.5 


88.0 


88.0 


85.0 


107 


68.0 


65.0 


67.0 


65.0 


67.0 


67.0 


108 


46.0 


42.0 


53.0 


52.0 


52.0 


53.0 


109 


— 


41.0 


29.0 


43.0 


32.5 


31.0 


110 


45.0 


36.0 


51.0 


38.0 


58.5 


52.0 


111 


48.0 


57.0 


62.0 


57.0 


66.5 


58.5 


112 


86.0 


89.0 


87.0 


90.0 


91.0 


89.0 


113 


70.0 


75.0 


74.0 


75.0 


74.0 


74.0 


114 


58.0 


69.0 


70.0 


69.0 


69.0 


69.0 


115 


37.0 


48.0 


32.6 


30.0 


36.0 


40.5 


116 


26.0 


26.5 


24.0 


27.0 


33.0 


24.0 


117 


80.0 


79.0 


78.0 


79.0 


79.0 


79.0 


118 


97.0 


95.0 


96.0 


86.0 


97.0 


96.0 


119 


47.0 


49.5 


54.0 


55.0 


56.0 


56.0 


120 


139.0 


138.0 


136.0 


136.0 


138.0 


136.0 


121 


88.0 


87.0 


87.0 


91.0 


87.5 


86.0 


122 


44.0 


46.0 


46.0 


47.0 


46.0 


45.0 


123 


16.0 


17.5 


17.0 


17.0 


20.0 


15.5 


124 


36.0 


10.0 


40.0 


41.0 


41.0 


41.0 


125 


15.0 


15.0 


14.7 


18.0 


17.0 


14.0 


126 


20.0 


22.0 


25.5 


19.0 


15.0 


26.0 
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Table 3.3. Concluded. 





Station 


Cruise 


Cruise 


Cruise 


Cruise 


Cruise 


Cruise 


Number 


One 


Two 


Three 


Four 


Five 


Six 


127 


17.0 


19.0 


19.0 


19.0 


21.0 


19.0 


128 


42.0 


39.0 


52.5 


39.0 


48.0 


47.0 


129 


44.0 


42.5 


42.5 


44.0 


44.0 


43.5 


130 


63.0 


64.0 


59.0 


65.0 


62.0 


61.0 


131 


122.0 


116.0 


65.0 


136.0 


62.0 


57.0 


132 


107.0 


105.0 


112.0 


105.0 


110.0 


104.0 


133 
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Dissolved nutrient samples were prepared by vacuum filtration within an 
hour of sample collection. A 47-mm diameter 0.45 y Millipore® filter (HAWP) 
held in a polycarbonate filter holder with a polypropylene filter flask was 
prewashed with 100 to 200 mL of demineralized water or sample water. New 125-mL 
polyethylene sample bottles with linerless closures were rinsed once with 
filtered sample prior to filling with filtrate for subsequent nutrient analyses. 
An aliquot was removed for dissolved orthophosphate and dissolved silica deter- 
minations after which the remainder was preserved with 1 mL/L concentrated 
sulfuric acid for total dissolved phosphorus determinations. 

Sampling methods on the USCGC Bramble were similar except that Niskin® 
bottles were loaded and subsequently retrieved directly from the deck. Sampling 
methods on the Canadian CSS Limnos used a Rosette® sampler to obtain water 
samples. 

Methods used by USEPA's Great Lakes National Program Office (GLNPO) can be 
found in greater detail in Rockwell et al. 1980. Table 3.4 provides a listing 
of the variables measured during cruises one, two, and four. A brief overview 
of the methods is presented below. CCIW methods are listed in Table 3.5 and 
compatible USE PA methods are summarized. 

Aesthetic . Reports of any unusual visual conditions were made. 

Air Temperature was measured with a dial scale bimetallic helix thermometer. 

Wind Speed and Direction . Readings from a permanently mounted Danf orth Marine 
Wind Direction and Speed Indicator were taken and recorded to the nearest degree 
and mile/hour while the vessel was stopped. 
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Table 3.5, The Canadian CSS Limnos was used on cruises three, five, and six. 
During these cruises samples for the following parameters were 
collected and are code referenced to the STAR Dictionary of Ana- 
lytical Methods and Units. On cruises one, two, and four where 
the R/V Simons was used and methods were compatible, the STAR 
Code is also referenced in the Table. 



Code 



Parameter Description 



Cruises Used 



1 Depth 

30 Secchi Disc Depth 

31 Forel ULE Secchi Disc Color 
100 Temperature of Water 

104 Temperature EBT 

124 Transparency (Percent) 

160 Specific Conductance at 25°C 

213 pH 

218 pH Temperature 

223 Total Filtered Alkalinity 

225 Integrated Particulate Organic Carbon 

238 Integrated Dissolved Organic Carbon 

245 Dissolved Oxygen Concentation 

247 % Saturation of Dissolved Oxygen 

254 Total Kjeldahl Nitrogen (Not Filtered) 

260 Total Phosphorus 

263 Dissolved Orthophosphorus 

264 Total Filtered Phosphorus 

267 Integrated Total Particulate Nitrogen 

270 Soluble Ammonia - N. 

276 Soluble Nitrate + Nitrite - N. 

295 Soluble Reactive Silica 

613 Integrated Chlorophyll a^ Uncorrected 

614 Integrated Chlorophyll a^ Corrected 
703 Coliform Fecal MF 

706 Streptococci Fecal MF 

709 Pseudomonas aeruginosa 

720 Aeorobic Heterotrophs - Standard Plak Count 

722 Total Micro Colonies Flores Count 



1. 2, 
1, 2, 



1. 2, 



1, 2, 



1. 2, 

1. 2, 
1, 2, 
1, 2, 
1, 2, 
1, 2, 
1, 2, 
1, 2, 
1, 2, 
1, 2, 
1. 2, 
1, 2, 
1. 2, 
1, 
1, 
1. 
1, 
1, 



4, 5 

4, 5 

5 

5 

4, 5 

4, 5 

5 

4, 5 

5 

5 

4, 5 

4 

4, 5 
4, 5 
4, 5 
4, 5 
4, 5 
4. 5 
4, 5 
4, 5 
4, 5 
4, 5 
4. 5 
4, 5 



6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 

6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
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Secchl Disc Depth was estimated at each station on all cruises by use of a 
standard 30-cm, all-white Secchi disc. 

Wave Height , Average wave height (valley to crest distance) was estimated at 
each station by the senior crew member on the bridge. 

Surface Water Temperature was determined by use of a mercury thermometer with 
0.1°C divisions. 

Temperature and Light Transmission Profiles , Vertical profiles were determined 
at each station from surface to bottom with a Martek® Model EBT/XMS electronic 
bathythermograph/transmissometer with a 1-meter folded light path. 

Turbidity , Turbidity was measured with a Turner Designs® digital nephelometer 
within 2 hours of sample collection. 

Specific Conductance , Specific conductance was measured at 20°C within 2 hours 
of sampling using a Barnstead® Model PM70CB conductivity bridge and a micro 
conductivity cell (VSI 3403, K=1.0), 

Dissolved Oxygen , Dissolved oxygen was measured on water samples from selected 
stations on some cruises. Analyses were made by the azid modification of the 
Winkler test (EPA 1974) or by a probe which was calibrated daily against the 
Winkler method, 

pH , pH analyses were made by electrometric measurement within 15 minutes 

of sample collection. Readings were recorded to the nearest 0.01 pH unit from 

an Orion® Model 701 pH meter equipped with an automatic temperature compensation 

probe and combination glass membrane with a silver/silver chloride internal 

electrode element. 
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Total Alkalinity as CaCO-^ ^ Total alkalinity was determined within 2 hours 
of sampling by titration of a 100-mL aliquot to pH 4.5 with 0.02 N H2SO4. 

Total Ammonia Nitrogen . Total ammonia nitrogen analyses were performed with a 
Technicon Autoanalyzer System II using a modification of Technicon's Industrial 
Method 154-7 IW/Tentative (Hiller and Van Slyke 1933). 

Total Nitrate and Nitrite Nitrogen . A Technicon® Autoanalyzer was used with 
Technicon's Industrial Method No. 158-71W (Armstrong et al. 1967). 

Dissolved Orthophosphate . Samples were analyzed for orthophosphate using a 
Technicon® Autoanalyzer System II and Technicon 's® Industrial Method 155-7 IW 
(Murphy and Riley 1962). 

Total Phosphorus and Total Dissolved Phosphorus . Conversion of the various 
forms of phosphorus to orthophosphate was by an adaptation of the acid persul- 
fate digestion method (Gales et al. 1966). The resulting orthophosphate was 
determined by the Technicon® Autoanalyzer System II and Technicon' s® Industrial 
Method 155-71W (Murphy and Riley 1962). 

Chloride . A Technicon® Autoanalyzer System II was used with Technicon 's® 
Industrial Method No. 99-70W (Zall et al . 1956, O'Brien 1962). 

Sulfate. Samples were analyzed for sulfate with a Technicon® Autoanalyzer using 
Technicon's® Industrial Method 118-71W (Lazrus _et al . 1965). 

Dissolved (Reactive) Silica . A Technicon® Autoanalyzer System II was used with 
Technicon' s® Industrial Method No. 186-72W/Tentative. 
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Chlorophyll "a" and Phaeophytin > Samples for chlorophyll (100 mL to 500 mL) 
were analyzed using a Turner® Model 1430 dual mono chroma tor spectrof luorometer 
(Strickland and Parsons 1972). 

Quality Assurance Used by GLNPO . Data quality assurance, evaluation, and con- 
trol were achieved by techniques given in detail in Rockwell et^ al. 1980. 
In brief, quality assurance consisted of check standards, reagent blanks, dupli- 
cate samples, split samples, and performance evaluation samples (unknowns). 
Maximum permissible shelf life was indicated for each analysis, and no data were 
taken from samples whose shelf life exceeded this value. The results of field 
split samples are given in Table 3.6. 

Data Base . Lake Huron was sampled extensively during 1980 by the U. S. Environ- 
mental Protection Agency (USEPA) and the Canada Centre for Inland Waters (CCIW). 
The data from the 138-station network are available from STORET, the CCIW data 
management group, or The University of Michigan. All nutrient data (TP, TDP, 
DRP, soluble NH4 , soluble NO2 + NO3, soluble SIO2) were performed by chemists 
from CCIW on all cruises. 
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Table 3,6« Differences between split sample analyses. 



Parameter 



Number of Mean Absolute Standard Deviation 
Splits Value of Differences of Differences 



Temperature (°C) 

Turbidity (HTU) 

Specific Conductance 
(umbo/ cm) 

pH (SU) 

Total Alkalinity (mg/L) 

Dissolved Ammonia (yg/L) 

Total Ammonia (yg/L) 

Total Nitrate + Nitrite (mg/L) 

Dissolved Nitrate + 
Nitrite (mg/L) 



176 


0.0057 


0.041 


180 


0.100 


0.348 


186 






182 


0.0437 


0.223 


184 






76 


0.39 


0.69 


118 


0.63 


2.09 


50 


0.0035 


0.0061 



88 



Dissolved Reactive 
Silica (rag/L) 



0.0031 



Dissolved Reactive 






Phosphorous (yg/L) 


266 


0.18 


Chloride (mg/L) 


52 


0.10 


Total Sulfate (mg/L) 


54 


0.17 



92 



0.012 



0.0051 

0.65 
0.15 
0.31 

0.014 
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CHAPTER FOUR 

DESCRIPTIVE LIMNOLOGY OF LAKE HURON 
by Russell Moll 

The 1980 study of Lake Huron was designed to provide extensive temporal and 
areal coverage of nutrient conditions throughout the lake. The sampling grid 
included 138 stations positioned in a low density pattern of approximately one 
station per 100 km^ (see Fig. 3.1). This grid was arranged such that in the 
open lake regions stations were located approximately 10 km apart. In nearshore 
regions, the North Channel, parts of Georgian Bay, and the mouth of Saginaw Bay 
the density of stations was twice as high as offshore, and the station pattern 
was somewhat irregular to accommodate the complex shoreline. Stations located 
in the mouth of Saginaw Bay were also at twice the density (one station/50 km^) 
of the rest of Lake Huron. This high density was designed to monitor the 
exchanges of material between the bay and the open lake. A separate sampling 
program was conducted in Saginaw Bay and is reported by Smith et al. (1983). 
Several stations were situated to monitor exchanges and/or inputs of materials 
into Lake Huron. Some of these locations were: the Detour Passage, the False 
Detour Passage, Mississagi Strait, the Straits of Mackinaw region, and the 
mouths of Georgian Bay and Saginaw Bay. 

Lake Huron was sampled during six cruises in 1980: April, May, June, July, 
September, and November. Two winter cruises were conducted in January and 
February 1981. The arrangement of these cruises was designed to provide data 
throughout the entire thermal cycle of one summer, i.e., from spring horaother- 
mous conditions through fall overturn to winter homothermous conditions. 
A major objective of this study was to develop accurate estimates of nutrient 
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conditions throughout Lake Huron just before thermal stratification (either 
therraocline formation or thermal bar formation). Cognizant of this objective, 
the cruise schedule was developed with the majority of the cruises occurring 
during the spring homothermous or thermal bar periods. 

The experimental design of this study was developed to yield the largest 
amount of information about areal and temporal nutrient distributions. The 
vertical distribution of nutrients was somewhat less intensely measured than the 
areal distribution. From four to seven depths were sampled at each station. 
Sampling depths always included a 1-m sample, a mid-water column sample, and two 
near-bottom samples. During periods of thermal stratification, samples were 
collected in the lower epilimnion, mid-metalimnion, and near the top of the 
hypolimnion in addition to the four depths mentioned above. This sampling 
program provided good but not detailed information about the vertical distribu- 
tion of nutrients. Samples for chlorophyll, particulate organic carbon, and 
particulate organic nitrogen were collected with a system which provided one 
integrated sample from the surface to 20-m depth. More detail on the sampling 
program can be found in Chapter Two, while the analytical methods are given in 
Rockwell et al. (1980). 

The sampling program included a complete suite of most dissolved and 
particulate nutrients, conservative ions, and particulate matter (Table 4.1). 
The complete sampling program for Lake Huron in 1980 included measurement of 
many variables in addition to nutrients, such as metals, organics, and plankton. 
This report is confined to the consideration of nutrients and related variables. 
But, because dissolved nutrient distributions are an important criterion in 
determining trophic status of lakes (Likens 1972), the inferences developed 
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Table 4.1. Chemical and physical variables measured during 1980 Lake Huron 
study. 



Alkalinity 

Ammonia 

Chlorophyll 

Dissolved Oxygen 

Kjeldahl Nitrogen 

Particulate Organic Carbon 

Particulate Organic Nitrogen 

pH 

Secchi Depth 

Soluble Nitrate-Nitrogen 

Soluble Reactive Phosphorus 

Soluble Reactive Silica 

Specific Conductance 

Temperature 

Total Filterable Phosphorus 

Total Phosphorus 

Transmissivity 

Turbidity 
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from analysis of the nutrient conditions have important implications for the 
entire ecosystem. 

Thermal Cycle 

The thermal cycle in the Great Lakes has a very large effect on the overall 
distribution of nutrients and biota (Ragotzkie 1978, Csanady 1978). Each ther- 
mal season produces a different mixing regime and hence distribution of dis- 
solved and suspended materials (Boyce 1974, Scavia and Bennett 1980). The 1980 
Lake Huron study covered the three thermal seasons typically observed in the 
Great Lakes; homothermous , thermal bar, and thermal stratification (Fig. 4.1). 

Cruise one (April) was conducted during inverse thermal stratification. 
Surface water temperatures were below 4'^C at almost all locations, with the 
lake-wide mean water temperature 2.2°C (Table 4.2). The surface (1-m) contour 
plots of temperature for the six cruises are shown in Figure 4.1. 

The second cruise (May) was conducted during early thermal bar conditions. 
Figure 4.1 shows large areas of nearshore surface water which had warmed above 
4°C, indicative of the onset of the thermal bar (Boyce 1974). Surface water 
temperatures increased between cruises one and two (Table 4.3), but the majority 
of Lake Huron was homothermous. 

Cruise three was conducted during late May and early June when thermal bar 
conditions were at an advanced stage. Surface temperatures increased consider- 
ably in the three weeks between cruises two and three (Table 4.4). Surface 
temperatures were in excess of 4°C over most of open Lake Huron and all of 
Saginaw Bay, Georgian Bay, and the North Channel (Fig. 4.1). More than 80% 
of the stations sampled during cruise three had surface water temperatures 
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Figure 4.1 One-meter temperature distribution during 1980. 
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Table 4.2. Cruise one descriptive statistics for 1-m samples. 



VARIABLE 
a.TEMPERAT 
9.PH 
1 1 .CONDUCTI 
17.NH3. 
18.N03 
ig.ORTHO PO 
20.S0L 5102 
21 .K NITROG 
26. TOTAL P 
28.CHL0R0PH 
29 . POC 
31 .SECCHI 



N MINIMUM 

122 .70000 

122 7.4400 

122 94.000 

105 0. 

122 . 19500 

106 .30000 -3 
122 .91000 
113 . 10700 
113 .20000 -2 
113 . 80000 
113 .10600 

71 1 . 0000 



MAXIMUM 
8 . 0500 
a . 4500 
245.00 
.44000 -1 
47500 
.35000 -2 
2 . 4300 
.46300 
.24500 -1 
10.000 
. 90400 
1 2 . 500 



MEAN 
2.2247 
8.0148 
193.05 
.31333 -2 
.28608 
.82547 -3 
1 .4533 
. 18848 
.55080 -2 
1 .6088 
.22512 
7.6986 



STD DEV 
.87189 
. 14050 
22.851 
.51255 -2 
.40026 -1 
.41954 -3 
.27340 
.61875 -1 
.27259 -2 
.96682 
. 13153 
2.7303 



Table 4.3. Cruise two descriptive statistics for l-m samples, 









MINIMUM 



















VARIABLE 


N 


MAXIMUM 


MEAN 




STD DEV 








8.TEMPERAT 


138 


2 . 0000 




1 1 -. 20O 


4.8123 




2.2279 








9.PH 


138 


7 . 4400 




8 . 7200 


7.9834 




. 17307 








1 1 .CONDUCTI 


138 


98 . 000 




237.00 


189.38 




25. 167 








16.NH3 


135 


. 10000 


-2 


.35000 -1 


.25259 


-2 


.38765 


-2 






17.N03 


138 


. 18000 




. 70500 


.27989 




.52352 


-1 






18.0RTH0 PO 


138 


. 20000 


-3 


. 39000 -2 


.62319 


-3 


. 36414 


-3 






19. SOL SI02 


138 


. 78000 




2 . 66C0 


1 .4472 




. 34598 








20. K NITROG 


138 


. 68000 


-1 


.49300 


. 16070 




.62538 


- 1 






21 .CHLORIDE 


128 


1 . 5COO 




8 . 9000 


5.0578 




.90594 








23 TOTAL P 


137 


. 33000 


-2 


. 12400 -1 


.52431 


-2 


. 16295 


-2 






25 CHLOROPH 


138 


. 90000 




7 , 8000 


1 7993 




.88338 








26. POC 


138 


. 24000 


-1 


. 72600 


. 19207 




. 1 1152 








28. SECCHI 


54 


3 . 5000 




16.000 


9.2722 




2.9438 







34 



Table 4,4 • Cruise three descriptive statistics for 1-m samples. 





VARIABLE 


N 


MINIMUM 




MAXIMUM 




MEAN 




STD DEV 







8.TEMPERAT 


138 


2 . 7900 




13.500 




8 . 0004 




3.0495 






9.PH 


138 


7 . 9900 




8 . 7400 




8.2785 




. 15346 






1 1 CONDUCT I 


138 


98.000 




234. OO 




190.96 




26. 180 






14.NH3 


138 


. 10000 


-2 


. 38000 


-1 


.23913 


-2 


.35503 -2 






15.N03 


138 


. 16100 




. 76500 




.27654 




.55868 -1 






16.0RTH0 PO 


138 


. 20000 


-3 


. 14000 


-2 


.70652 


-3 


. 18765 -3 






17. SOL SI02 


138 


.51000 




2 . 6800 




1 .3854 




.35510 






18. K NITROG 


138 


. 60000 


-1 


.32200 




. 12445 




.33915 -1 






19. CHLORIDE 


138 


1 .3000 




7 . 2000 




5. 1065 




.95994 






21 . TOTAL P 


138 


. 28000 


-2 


. 73000 


-2 


.45406 


-2 


.87451 -3 






23.CHL0R0PH 


138 


. 90000 




4 . 7000 




1 .9565 




.60145 






24.P0C 


138 


. 11300 




. 60000 




.22330 




.68143 -1 






2S.SECCHI 


73 


2.5000 




1 2 . 000 




6.8425 




2.0427 





Table 4.5. Cruise four descriptive statistics for 1-m samples, 




















--■ 










VARIABLE 


N 


MINIMUM 


MAXIMUM 




MEAN 


STD DEV 








8.TEMPERAT 


138 


1 1 . 900 




22.000 




18.819 




1 . 706 1 








9 .PH 


138 


7 . 7800 




9.2200 




3. 1885 




. 14852 








1 1 .CONDUCT I 


138 


101 .00 




263 . 00 




188.73 




24.208 








14.NH3 


138 


. 10000 


-2 


. 1 1000 


-1 


.31812 


-2 


. 17561 


-2 






15.N03 


138 


. 13800 




. 32900 




.23962 




. 30087 


-1 






16.0RTH0 PO 


136 


. 10000 


-3 


. 17000 


-2 


. 39338 


-3 


.25912 


-3 






17. SOL SI02 


138 


0. 




2 . 1 300 




1 .0475 




. 34200 








18. K NITROG 


138 


.87000 


- 1 


.25700 




. 15324 




.30401 


- 1 






19. CHLORIDE 


106 


1 . 3000 




7 . 0000 




5.048 1 




.94760 








2 1 .TOTAL P 


138 


. 23000 


-2 


. 90000 


-2 


.41572 


-2 


. 12776 


-2 






23.CHL0R0PH 


138 


. 30000 




9 . 4000 




1 .2420 




.31776 








24 POC 


138 


. 11000 


-1 


.34 100 




.20016 




.53643 


-1 






26.SECCHI 


93 


1 . 5000 




1 7 . 600 




9.0935 




2 . 3806 
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in excess of 4*^0. Nearshore thermal warming had progressed to the stage where 
early thermal stratification had developed up to 10 km offshore. Figure 4,1 
shows large nearshore areas during cruise three with surface temperatures in 
excess of 10°C. 

Lake Huron was completely thermally stratified during cruise four (late 
July), Surface temperatures had reached a mean of 18.8°C (Table 4,5) . The 
depth of the thermocline varied across the lake but ranged from 13 m in the 
southern end of the lake to 5 m in the deepest part of the lake near station 48. 
Figure 4.1 shows that there was a large range of surface temperatures across 
Lake Huron from less than IS^'C to greater than 21 °C. 

Thermal stratification continued through cruise five (mid-September), 
although some cooling had occurred between cruises four and five. Table 4.6 
shows that surface temperatures were lower in September (mean = 16.2°C) than 
July. Figure 4.1 shows that the highest surface temperatures were found in the 
southern and central parts of Lake Huron, while the lowest temperatures were 
nearshore in the north channel. The average depth of the thermocline was 23 m 
in September. 

Cruise six was conducted during late October and early November. Although 
surface water temperatures had cooled considerably from the previous cruise in 
September, Lake Huron was still fully stratified. Figure 4.1 shows uniform 
surface temperatures with a mean of 8.2°C (Table 4.7). The thermocline at 85 m 
had deepened considerably from fall cooling and mixing since cruise five. Those 
parts of Lake Huron which have water depths in excess of 85 m could still be 
characterized as completely stratified, but the thermal gradient between the 
surface and bottom was only 4°C. 
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Table 4.6. Cruise five descriptive statistics for l-m samples. 





VARIABLE 


N 


MINIMUM 




MAXIMUM 




MEAN 


-- 


STO OEV 






8.TEMPERAT 


138 


1 1 . 400 




22 . 600 




16.229 




2.6205 






9.PH 


138 


7 . 7000 




8 . 6600 




8.2254 




. 16387 






1 f CONDUCT I 


137 


95.000 




234.00 




182.23 




16.075 






1 4 . NH3 


138 


. 10000 


-2 


. 13000 


-1 


. 15362 


"2 


. 15291 -2 






15.N03 


138 


. 17200 




. 29300 




.23507 




.20620 -1 






IG.ORTHO PO 


138 


. 10000 


-3 


. 1 7000 


-2 


.73116 


-3 


.24903 -3 






17. SOL SI02 


138 


.56000 




2. 1500 




1 . 0209 




.27873 






18. K NITROG 


138 


. 12100 




. 33400 




. 1764 1 




.38985 -1 






19. CHLORIDE 


138 


1 . 4000 




7 . 3000 




5. 1355 




.68122 






2 1 .TOTAL P 


138 


,29000 


-2 


. 95000 


-2 


. 435^22 


-2 


. 10250 -2 






23.CHLOROPH 


138 


. 80000 




2.8000 




1 .4014 




. 34809 






24.P0C 


138 


. 14500 




. 48400 




.22257 




.53869 -1 






26.SECCHI 


68 


1 . 5000 




11 .000 




6.2485 




2.3193 





Table 4.7. Cruise six descriptive statistics for 1-m samples. 





VARIABLE 


N 


MINIMUM 





MAXIMUM 


MEAN 





STD DEV 


-- 




8.TEMPERAT 


138 


5.2100 




1 1 . 600 


8.2165 




1 .4167 






9.PH 


138 


7. 1500 




8 . 5600 


8. 1822 




. 16022 






11 .CONOUCTI 


137 


99 . 000 




267.00 


198.82 




19.363 






14.NH3 


138 


. 10000 


-2 


. 10000 -1 


.268 12 


-2 


. 17299 


-2 




15.N03 


138 


. 16300 




. 3 1 20O 


.265 10 




.251 12 


-1 




16.0RTH0 PO 


138 


. 20O0O 


-3 


.28000 -2 


.5884 1 


-3 


.30401 


-3 




17. SOL SI02 


138 


. 76000 




2.2100 


1 .2951 




.27798 






18. K NITROG 


138 


. 10300 




. 50400 


. 17501 




.53133 


-1 




19. CHLORIDE 


138 


1 . 3000 




8.0000 


5.2522 




.75283 






2 1 .TOTAL P 


138 


. 33000 


-2 


. 1 1 300 - 1 


.52486 


-2 


. 15121 


-2 




23.CHL0R0PH 


138 


. 60000 




4 . 7000 


1 .5348 




.65698 






24.P0C 


138 


. 1 1400 




. 59600 


.21285 




. 88008 


- t 




26.SECCHI 


69 


1 . 0000 




9 . 0000 


4.2754 




1 .8875 





37 



Conductivity and Chloride Distribution 

The distribution of the conservative tracers conductivity and chloride did 
not show a seasonal trend. These two variables should be unaffected by either 
temperature or the biota (Hutchinson 1957, Wetzel 1975). The areal distribution 
for conductivity and chloride was considered a result of two factors: 
(1) sources of high or low conductivity water, and (2) surface mixing of water 
masses. The second factor, surface mixing, is affected by the thermal cycle. 
For example, the thermal bar tends to act as a barrier to horizontal water mass 
exchange (Boyce 1974). But, Scavia and Bennett (1980) have shown that the 
thermal bar enhances vertical mixing. Figures 4.2 and 4.3 show the surface 
distribution of conductivity and chloride respectively for cruises one through 
six (cruise one chloride data are missing). 

Figures 4.2 and 4.3 show a high degree of agreement in the areal distribu- 
tion of conductivity and chloride, which was confirmed by correlation analyses 
(see Chapter Seven). A previous study (Schelske et al. 1980, Moll et_ al. 1980) 
attributed high chloride levels in southern Lake Huron to inputs from Saginaw 
Bay. Likewise, in the Straits of Mackinac, high chloride water was a result 
of inflowing Lake Michigan water (Moll £t al . 1976). Low chloride water in the 
northeastern Straits region was attributed to inflowing Lake Superior water. 

The areal distribution of conductivity during cruise one shows the effect 
of the two high conductivity and one low conductivity sources. High chloride 
water was observed on the southern edge of Saginaw Bay along the north shore 
of Michigan's "thumb" and in the Straits of Mackinaw. Low conductivity water 
dominated the western end of the north channel, the entire northwestern end 
of open Lake Huron, and parts of Georgian Bay (Fig. 4.2). Several small regions 
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Figure 4.2 One-meter conductivity distribution during 1980. 
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Figure 4.3 One-meter chloride distribution during 1980. 
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of high conductivity water were observed in the nearshore regions; these high 
conductivity regions were interpreted as outflow of small rivers which often 
carry high chloride loads in the early spring (Schelske et al« 1980). The low 
conductivity area of Georgian Bay was attributed to a low chloride river dis- 
charge. 

During cruise two the same conductivity pattern was observed as cruise one, 
except that nearshore conductivities were higher, raising the lake-wide, surface 
conductivity mean from cruise one (Fig. 4.2, Tables 4.2 and 4.3). The higher 
nearshore conductivities of cruise two were attributed to high river discharges 
during the spring thaw. The areal distributions of conductivity and chloride 
during cruise two were very similar (Figs. 4.2 and 4.3). 

Cruises three through six produced very similar conductivity and chloride 
distributions (Figs. 4.2 and 4.3). Conductivity and chloride values in southern 
and open Lake Huron were uniform at approximately 205-210 ymho/cm and 5.-5 to 
6.0 mg/L, respectively. Regions of high conductivity and chloride water in Lake 
Huron's southern basin were confined to the edge of Saginaw Bay along Michigan's 
*' thumb." This area appears to be a location where Saginaw Bay water typically 
mixes into open Lake Huron (Schelske et al. 1980). Two regions of low chloride 
and conductivity water were observed during cruises three to six; the North 
Channel and adjacent open lake where Lake Superior water mixes into Lake Huron, 
and the north edge of Georgian Bay. The distribution of low chloride and con- 
ductivity water in the North Channel and adjacent region changed somewhat from 
cruise to cruise and was attributed to wind-induced mixing. For example, the 
cruise four areal distribution shows an almost north-south gradient between the 
low and average chloride waters while cruise six showed an east-west gradient. 
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Despite changes in distribution, conductivity and chloride concentrations 
changed very little from April through November, 

Seasonal Nutrient Cycles 

Soluble Silica — Research over the past twenty years has shown soluble 
silica undergoes a consistent annual cycle in the Great Lakes (Schelske and 
Stoermer 1971, Kilham 1971, Rousar 1973). This annual silica cycle has been 
attributed to increased phosphorus levels which promote diatom growth and hence 
silica uptake (Schelske and Stoermer 1971). The annual soluble silica cycle 
in surface waters begins with high soluble levels in the spring, goes through 
a severe depletion in the summer, and slow replenishment throughout the fall 
and winter (Rousar 1973, Rousar and Beeton 1973). 

This annual pattern of soluble silica concentrations was observed in Lake 
Huron, but to a much lesser extent than in Lake Michigan (Schelske and Roth 
1973). Lake Huron water is composed primarily of a mixture of Lake Michigan 
and Lake Superior outflows. Lake Superior waters have relatively high and 
consistent silica levels, while Lake Michigan waters readily display an annual 
cycle as described above (Rockwell et al. 1980). Because the Lake Superior 
inflow vastly exceeds the Lake Michigan inflow, soluble silica cycles in north- 
western Lake Huron showed very little seasonal change (Fig. 4.4). 

Soluble silica concentrations in central Lake Huron, southern Lake Huron, 
and Georgian Bay somewhat followed the annual cycle described above in conjunc- 
tion with the thermal cycle. Figure 4.4 shows high surface silica levels for 
cruise one with a mean of 1.5 mg/L Si02 (Table 4.2). The areal distribution 
of silica was the same throughout all cruises: low soluble silica levels near 
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Figure 4.4 One-meter soluble silica distribution during 1980. 
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the Straits of Mackinaw (Lake Michigan water), and high levels in the North 
Channel and adjacent Lake Huron (Lake Superior water), 

A slight decline in surface silica levels was observed between the first 
and third cruises as the mean dropped to 1.4 mg/L« Thermal bar conditions were 
present during both cruises two and three. Algal utilization of silica was 
indicated between the thermal bar and shore by a band of low silica water along 
most of the shoreline. Open Lake Huron soluble silica levels declined only very 
slightly (0.1 mg/L) from spring concentrations. 

Soluble silica concentrations continued to decline through cruises four and 
five which were conducted during thermal stratification. The surface means of 
soluble silica were 1.0 mg/L for cruises four and five. Figure 4.4 shows a 
continuation of the pattern of lowest silica levels nearshore during cruise 
four. Cruise five, on the other hand, had fairly uniform areal silica concen- 
trations. 

Some replenishment of soluble silica levels from fall mixing was evident 
during cruise six. Silica levels increased to a mean surface value of 1.3 mg/L. 
In summary, soluble silica displayed an annual cycle similar to the one observed 
in Lake Michigan, but the level of depletion was much lower than has been 
observed in Lake Michigan (Rockwell et al. 1980). 

Orthophosphate — Orthophosphate has been considered the limiting nutrient 
for algal growth in most of the Great Lakes (Beeton 1969). A consequence of the 
high algal demand for orthophosphate is that concentrations tend to be very low, 
especially in the epilimnion. Figure 4.5 shows the surface concentration of 
orthophosphate was low throughout all six cruises. The mean surface orthophos- 
phate concentration ranged from 0.8 yg/L for cruise one to 0.4 ug/L for cruise 
four. 
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Figure 4.5 One-meter soluble reactive phosphorus distribution during 1980. 
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Areal patterns of orthophosphate were not evident during any of the 
cruises. In a few instances, high levels of orthophosphate were observed near- 
shore, probably from river discharges. High orthophosphate water was not 
observed in the mouth of Saginaw Bay. Studies by Moll et al. (1980) showed that 
large algal populations in inner Saginaw Bay reduced most soluble nutrient 
levels below open lake ambient concentrations. Areal patterns of orthophosphate 
could be characterized as uniformly low throughout the entire lake. Similarly, 
seasonal patterns could be called uniformly low throughout all six cruises. 

Nitrate-Nitrogen — In predominantly phosphorus-limited oligotrophic lakes 
such as Lake Huron, nitrate-nitrogen is not considered a limiting nutrient 
(Wetzel 1975). Nonetheless, areal and seasonal patterns of nitrate-nitrogen 
were observed in Lake Huron (Fig. 4.6). Surface concentrations of soluble 
nitrate-nitrogen ranged from a mean of 0.28 mg/L for cruises one through three 
to 0.24 mg/L for cruise four. This small range was indicative of a relatively 
low level of nitrogen utilization by phytoplankton. 

Areal patterns of nitrate distribution were primarily determined by the 
inflow of high or low nitrate water into Lake Huron. Davis et al. (1980) 
measured extremely high nitrate-nitrogen concentrations along the Province of 
Ontario shore during thermal bar conditions. They attributed the high concen- 
trations to high nitrate river outflows during spring runoff. The highest 
nitrate concentrations of all six cruises were found in the same region, the 
Ontario shore of southern Lake Huron (Fig. 4.6). The high nearshore nitrate 
levels were observed during and just prior to thermal bar conditions of cruises 
one to three. During cruises five and six, surface nitrate-nitrogen concentra- 
tions were uniform across the entire lake. Thus, the nitrate levels along the 
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Figure 4.6 One-meter nitrate plus nitrite distribution during 1980. 
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Ontario shore were only slightly elevated above ambient open lake levels. 
Nitrate-nitrogen annual cycles and areal distribution appear to be controlled 
primarily by exogenous nitrogen sources and not algal utilization. 

Ammonia-Nitrogen — The nutrient dynamics of ammonia tend to fall between 
orthophosphate and nitrate (Fogg 1975). While ammonia is not considered a 
limiting nutrient, it is a highly available form of nitrogen for algal uptake 
(Eppley et al. 1969). As a result, ammonia concentrations remain at constant 
low levels. Ammonia originates from one of two sources, either aquatic animal 
excretion (zooplankton and fish excretion) or local domestic sewage outfalls. 
Large sewage outfalls are not abundant in Lake Huron. The annual range in 
ammonia concentration was small, with a mean of 3.1 yg N/L for cruises one and 
four to 1.5 yg N/L for cruise five (Tables 4.2 to 4.7). Figure 4.7 shows that 
ammonia levels were consistent and low across the open lake for all cruises. 
Occasionally a nearshore area would have higher ammonia levels than the open 
lake, especially in southern Lake Huron. Ammonia areal distribution was analo- 
gous to orthophosphate distribution; concentrations were consistently low and 
uniform and probably kept at low levels by algal utilization. 

Total Nitrogen and Phosphorus — The areal distribution of total (Kjeldahl) 
nitrogen and total phosphorus for the six research cruises are shown in Fig- 
ures 4.8 and 4.9, respectively. The seasonal cycle and areal distribution 
of these two variables are tied closely to phytoplankton biomass and produc- 
tivity (Paerl et al. 1975). Because Lake Huron is primarily oligotrophic, it 
supports only low levels of algal growth (Schelske £t al . 1980). Therefore, 
throughout most of the year and across much of the open lake total nitrogen and 
phosphorus levels were uniformly low. But, high levels of total nitrogen and 
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Figure 4.7 One-meter airanonium ion distribution during 1980. 
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Figure 4.8 One-meter Kjeldahl nitrogen distribution during 1980. 
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Figure 4.9 One-meter total phosphorus distribution during 1980. 
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phosphorus were occasionally found in the nearshore zone where soluble nutrient 
levels were elevated. 

Cellular algal phosphorus and nitrogen levels are determined, in part, 
by external nutrient sources (Fogg 1975). As mentioned above, orthophosphate 
concentrations were uniformly low across Lake Huron throughout all six cruises. 
Soluble nitrogen levels were uniform across the six seasons, but occasionally 
elevated in nearshore areas. Total nitrogen and phosphorus reflect the same 
pattern; consistent across all six cruises (Tables 4.2-4.7), and uniform across 
most of the lake except in some nearshore areas. 

Seasonal Biological Cycles 

Particulate Organic Carbon — In oligotrophic lakes such as Lake Huron, 
sources of particulate organic carbon originate primarily from biological pro- 
ductivity (Wetzel 1975). Thus algal cells and zooplankton contribute most, 
if not all, of the particulate organic matter to the open lake. Particulate 
matter from river runoff creates high levels of particulate organic carbon in 
the nearshore zone. 

The areal distribution of particulate organic carbon (Figure 4.10) was 
fairly consistent in the open lake and elevated concentrations near shore. 
The open lake concentrations , which are dominated by phytoplankton and zooplank- 
ton, showed some annual pattern; highest particulate carbon concentrations were 
observed in the late spring and early summer, the time of the spring algal 
bloom. But the total annual range in particulate organic carbon in the open 
lake was small. 
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Figure 4.10 One-meter particulate organic carbon distribution during 1980. 
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Particulate organic carbon distribution in the nearshore areas apparently 
was affected by river runoff and locally high concentrations of phytoplankton. 
Particulate organic carbon levels were consistently high in the following near- 
shore areas: the Ontario shore of southern Lake Huron, the southern edge of 
outer Saginaw Bay, eastern Georgian Bay, and the southern shore of the Straits 
of Mackinaw. The areal extent of these elevated particulate organic carbon 
regions suggested that river runoff alone was unlikely to produce the elevated 
concentrations. Rather, the high nutrient loads from river runoff were probably 
stimulating algal productivity and increasing particulate carbon levels. Davis 
et al. (1980) found this mechanism to operate along the Ontario shore of south- 
ern Lake Huron in 1974. Moll et al. (1980) found a similar series of events in 
outer Saginaw Bay in 1975. Comparing the particulate organic carbon distribu- 
tion to the chlorophyll distribution indicates a large portion of organic carbon 
contains chlorophyll (i.e., phytoplankton). 

Chlorophyll — The distribution of chlorophyll bioraass for the six cruises 
is shown in Figure 4.11. Chlorophyll biomass is closely tied to phytoplankton 
concentration, especially in surface samples. Furthermore, chlorophyll biomass 
has been used as a term in several models to predict primary productivity and 
as a response variable to changing nutrient concentrations (Bierman and Dolan 
1981). Because of the relationships among chlorophyll, nutrients, and primary 
productivity, chlorophyll distributions have been thoroughly analyzed on both 
temporal and spatial scales. 

Figure 4.11 presents the areal distribution of surface (1-20 m) chlorophyll 
levels for the six 1980 cruises. An annual cycle of surface chlorophyll values 
has typically been observed throughout the Great Lakes as follows: a spring 
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Figure 4.11 One-meter chlorophyll distribution during 1980. 
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bloom of phytoplankton follows the annual minimum values during the winter; 
relatively low surface chlorophyll levels are typical of mid-summer, followed by 
a small fall algal bloom (Glooschenko and Moore 1973, Fee 1976, Munawar and 
Burns 1976, Vollenweider et al, 1974). The relatively low and invariant surface 
chlorophyll concentrations found in Lake Huron were indicative of oligotrophic 
lakes. Tables 4.2-4.7 show the mean surface chlorophyll values ranged from 
1.24 ug/L during cruise four to 2.00 yg/L during cruise three. The annual 
pattern of chlorophyll in Lake Huron was consistent with the cycle described 
above, but the difference between spring bloom conditions and summer minimum was 
only 0.76 yg/L. The hypothesis presented by Moll and Stoermer (1982) indicates 
that the difference between the size of the spring bloom and annual minimum is 
just as important as the total biomass. Lake Huron's uniform annual chlorophyll 
levels further suggest oligotrophic conditions. 

The areal distribution of chlorophyll is often used as an indication of 
areas of high algal growth due to nutrient loading (Holland and Beeton 1972, 
Robertson et^ ad. 1971). Considered in this fashion, the chlorophyll distribu- 
tions shown in Figure 4.11 again reflect the oligotrophic character of Lake 
Huron. During cruises one, two, three, and six, elevated chlorophyll concentra- 
tions were observed along the Ontario shore of southern Lake Huron. Other near- 
shore areas which exhibited elevated chlorophyll levels were: southern Saginaw 
Bay along Michigan's "thumb," extreme southern Lake Huron, eastern Georgian Bay, 
and the Straits of Mackinaw. During cruises one, two, and four, there were one 
or two stations with high chlorophyll levels (>7.0 yg/L). But, with few excep- 
tions, chlorophyll concentrations in the nearshore regions mentioned above 
rarely exceeded lake-wide means by 3.0 yg/L. 
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One nearshore area of Lake Huron that deserved particular attention was the 
tip of Michigan's "thumb." In 1974 and 1975, this region was characterized by 
very high chlorophyll levels as a consequence of eutrophication in Saginaw Bay 
(Moll et al. 1980). Since that study, remedial measures have been used in an 
attempt to control Saginaw Bay pollution. The results shown in Figure 4.11 show 
very little increase in chlorophyll levels in this region above the adjacent 
open lake conditions. 

Secchi Disc — Water transparency as measured by the Secchi disc technique 
cannot be plotted for the 1980 study. Large blocks of the Secchi disc readings 
are missing because readings could not be made at night when the ships were 
operating 24 hours a day. Nonetheless, the statistics given in Tables 4.2-4.7 
allow some comparison between 1980 Secchic depths and previous studies. 

Secchi disc transparency typically follows the annual cycle of chlorophyll 
concentrations (Ladewski and Stoermer 1973). This close correspondence between 
chlorophyll and Secchi depth readings has led many researchers to suggest Secchi 
readings can determine trophic status, chlorophyll concentrations, nutrient 
concentrations, etc. (Carlson 1977). Although Secchi depths can indicate water 
clarity with high reproducibility, they cannot reliably estimate trophic status, 
chlorophyll biomass , etc. In Lake Huron, Secchi depths ranged from 1.0 to 16.0 
meters (Tables 4.2-4.7) during 1980. Smallest Secchi depth readings were found 
during the third cruise (June), which occurred during thermal bar/spring bloom 
conditions. Largest Secchi depths were observed in the open lake during the 
fourth (July) cruise. It should be noted that the pattern of missing data has 
a large influence on the Secchi results. Many missing observations from the 
middle of Lake Huron would bias the mean depth toward smaller depths measured 
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in the nearshore and vice versa. As a consequence of the possible bias created 
by the missing observations, the Secchi depth data provide few inferences beyond 
the relationship with chlorophyll bioraass. 

pH — Fluctuations in pH values have been used as an estimate of free-water 
primary productivity (Hall and Moll 1975 )• This technique hypothesizes that 
falling pH reflects an increase in dissolved CO2 due to plant and animal respi- 
ration; increasing pH indicates increasing levels of dissolved oxygen from 
primary productivity. However, the pH method of estimating primary productivity 
lacks sufficient accuracy for use in oligotrophic environments such as Lake 
Huron. Furthermore, for the pH method to be useful, one station must be 
observed throughout one 24-hour cycle. Finally, pH measurements from different 
stations at different times of the day are not comparable. 

In a well-buffered, oligotrophic system such as Lake Huron, pH changes were 
expected to be minimal. Tables 4.2-4.7 and Figure 4.12 confirm this expecta- 
tion. pH changes throughout the year were small: mean surface readings of 7.98 
to 8.22. Likewise, areal distributions showed uniformity across the entire lake 
in all six cruises. These small changes coupled with the problem of readings 
made at different locations at different times of day prevented any ecological 
inferences being drawn from the pH data. More meaningful analysis of the pH 
data came from consideration of pH readings with depth (see below). 

Vertical Distribution of Temperature, 
Nitrogen, Phosphorus, and Silicon 

One station representative of each of the four major basins of Lake Huron 
was used to document the vertical distribution of temperature and nutrients 
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Figure 4.12 One-meter pH distribution during 1980, 
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between April and November 1980. The four stations plotted in Figures 4.13 to 
4.16 are stations 15, 44, 78, and 117 from the southern basin, central basin, 
North Channel, and Georgian Bay, respectively. Results are plotted for tempera- 
ture, dissolved silica, nitrogen (ammonium ion, Kjeldahl nitrogen, nitrate + 
nitrite), and phosphorus (total phosphorus, total dissolved phosphorus, ortho- 
phosphorus). Data from station 78 (North Channel) are not available for the 
first cruise (April) because the station was ice covered and inaccessible to the 
research vessel. Additional data are available from station 15 (southern basin) 
during January and February of 1981. The results presented below are for one 
variable at a time rather than by geographic regions; this presentation is used 
because the vertical profiles for each variable are similar among all four 
geographic regions . 

Temperature — Thermal conditions in Lake Huron during the April cruise 
could be described as either late winter or early spring. All locations had 
inverse thermal stratification with warmer, more dense 4°C water at the bottom 
of the lake. This inverse stratification was most pronounced at station 44 
(central basin), but was present in the other geographic regions as well. 
Homothermous conditions were present during cruise two (May) at all locations 
and during cruise three (late May-early June) at all locations except the North 
Channel. At the latter location, signs of early thermal stratification were 
present. These signs were an extremely shallow thermocline (5-6 m) and cold 
surface temperatures (10-11°C). 

Thermal stratification was present throughout Lake Huron during cruises 
four through six (July to November), with the exception of station 78 (North 
Channel) in November. The thermal cycle was similar in all four basins and was 
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as follows: intense stratification was present during cruise four (July), but 
the epilimnion was relatively thin (5-10 m) . Stratification continued through 
September, but the epilimnion was much thicker (15-30 m) by cruise five. 
The epilimnion became very thick in November (50-80 m) , although the temperature 
difference between the epilimnion and hypolimnion was relatively small (5C°). 
Station 78 in the North Channel had warm (9C**) homothermous conditions in 
November, The two winter cruises (January and February) in the southern basin 
showed a return to inverse stratification. 

Dissolved Silica — The vertical distribution of dissolved reactive silica 
followed a similar annual pattern in all four basins of Lake Huron, During 
periods of inverse thermal stratification, dissolved silica concentrations were 
higher near the bottom than at the surface. As the lake became homothermous and 
was mixed from surface to bottom, dissolved silica became evenly distributed 
throughout the water column. The onset of thermal stratification created a 
large difference between surface and bottom dissolved silica concentrations, 
with surface waters having much lower levels. The difference in dissolved 
silica concentrations between surface and bottom waters increased from July 
to September. As Lake Huron cooled during the fall and early winter, dissolved 
silica was somewhat replenished in the surface waters. But a difference between 
the lower surface concentrations and higher bottom values persisted through 
February in the southern basin. 

Dissolved silica profiles at station 78 in the North Channel differed 
slightly from the general trend described above. During cruise three, higher 
concentrations of dissolved silica were observed at the surface rather than 
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the bottom. Furthermore, dissolved silica concentrations were uniform from 
the surface to the bottom during cruise six. 

The annual cycle of vertical profiles of dissolved silica has been observed 
in other Great Lakes, most notably Lake Michigan (Schelske and Stoermer 1971, 
Rockwell et al. 1980). The vertical patterns of dissolved silica in Lake Michi- 
gan have been attributed to intense silica utilization by diatoms and silico- 
flagellates (Schelske and Stoermer 1971). The onset of thermal stratification 
is accompanied by the spring diatom bloom which rapidly utilizes most of the 
available surface dissolved silica. Algal cells are then transported downward 
by either sinking or in fecal pellets, increasing near-bottom silica concentra- 
tions. This downward flux of silica continues throughout the summer until 
differences in surface to bottom concentrations are large. These same processes 
appeared to occur throughout Lake Huron in 1980 because the annual cycle of 
dissolved silica was similar to that of Lake Michigan. One interesting aspect 
of the Lake Huron dissolved silica cycle is that a uniform vertical profile was 
not apparent until just before the onset of summer thermal stratification 
(Figure 4.13 - Cruise 1 compared with winter cruises 1 and 2). 

Dissolved Nitrogen — Figures 4.13 to 4.16 show that the differences in 
dissolved nitrogen concentrations (NO3 + NO2 and NH3) between surface and bottom 
waters are relatively small. Furthermore, dissolved nitrogen was uniformly 
distributed throughout the water column by cruise one (April) and remained 
uniformly distributed until cruise four (July). Dissolved nitrogen follows 
approximately the same annual patterns as dissolved silica for similar reasons; 
summer phytoplankton growth reduces surface dissolved nitrogen concentrations, 
while bottom concentrations increase from accumulation of decaying material 
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(Wetzel 1975). Ammonium concentrations are somewhat more controlled by animal 
(mostly zooplankton) excretion. Therefore, the vertical distribution of ammo- 
nium usually differs from the vertical distribution of nitrate + nitrite during 
periods of active zooplankton growth. July (cruise 4) was a period of high 
zooplankton abundances in Lake Huron (Evans 1983). Because nitrogen is usually 
not a limiting nutrient for algal growth in the Great Lakes (Lin and Schelske 
1981), surface nitrogen depletion is rarely as intense as dissolved silica. 
Nitrogen would not be expected to limit algal growth in southern Lake Huron 
because of large loadings during the spring runoff (Davis ^ al . 1980). 

Particulate Nitrogen and Particulate Phosphorus — Total Kjeldahl nitrogen 
(particulate) had a more irregular pattern of vertical profiles than did dis- 
solved nitrogen. Particulate forms of nutrients tend to have highly variable 
vertical profiles during periods of thermal stratification (including inverse 
thermal stratification). Primary production (algal photosynthesis) is the major 
process that converts dissolved nutrient pools into particulate pools (Wetzel 
1975). Once the nutrient has gone through this phase change, a variety of 
processes affect its distribution as particulate nitrogen and phosphorus. These 
processes include sinking, advection, grazing, and regeneration. Because each 
of these processes can dominate particle dynamics at different locations and 
times, the vertical distributions of particulate nitrogen and phosphorus in Lake 
Huron are usually non-uniform and highly variable. In general, processes 
affecting particulate phosphorus and nitrogen pools are important only during 
periods of high primary productivity (spring and summer). As a result, vertical 
profiles of particulate nitrogen and phosphorus showed no change among the early 
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spring cruises (one to three) and large differences among the late spring and 
summer cruises (four and five). 

Phosphorus — The vertical distributions of dissolved orthophosphate and 
total dissolved phosphorus followed a pattern similar to dissolved nitrogen. 
The distribution of soluble phosphorus was uniform during homothermous condi- 
tions (May and June). During stratified conditions, phosphorus concentrations 
were always lower in surface waters than in bottom waters; this surface to 
bottom difference was also present during periods of inverse thermal stratifi- 
cation. 

Dissolved phosphorus is usually regarded as the limiting factor for algal 
growth in Lake Huron (Lin and Schelske 1981). As a result, this nutrient is 
almost always at low concentrations and rapidly recycled throughout the epilim- 
nion. Available phosphorus pools are small and rapidly utilized. These proc- 
esses keep phosphorus concentrations uniformly low throughout the water column, 
and reduce the magnitude of the surface to bottom differences in dissolved 
phosphorus concentrations. The vertical profiles of dissolved and total dis- 
solved phosphorus in Lake Huron (Figs. 4.13 to 4.16) show relatively uniform 
distributions but low concentrations. 

Conservation of Nutrient Mass — The processes of advection and primary 
production can serve to cause overall seasonal depletion of a nutrient in one 
region of Lake Huron. The mechanisms work as follows: Photosynthesis causes 
a major shift from the dissolved to particulate phase of most nutrients in the 
summer. Currents could then move the particulate nutrients away from the 
regions of high photosynthesis. The net result is a local, seasonal depletion 
of nutrients. This process can be partially evaluated by comparing the water 
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column total for each nutrient among the six cruises. Figures 4.13 to 4.16 show 
that the only region of net loss or gain was for dissolved silicon at station 78 
(North Channel). This station shows a net loss of dissolved silica between 
cruises two and six. 

Summary 

The limnology of Lake Huron during 1980 indicates two major properties 
of the ecosystem: this is an oligotrophic ecosystem, and the lake has several 
important external inputs of nutrients. The oligotrophic characteristics of the 
ecosystem are reflected in two ways. First, nutrient concentrations were low in 
1980, compared to most other lakes, indicating low levels of internal and exter- 
nal loadings. Second, the annual ranges of most nutrients were small. These 
small ranges suggest small seasonal loadings which in turn stimulate small algal 
crops. 

The external inputs to Lake Huron take several forms based on location 
and season. Water from Lakes Michigan and Superior flow into Lake Huron and 
mix in the vicinity of the North Channel and Straits of Mackinac. These waters 
tend to have contrasting properties, and Lake Huron ultimately takes on charac- 
teristics as a mixture of the two water masses. Georgian Bay is often dominated 
by runoff from rivers flowing across the Canadian Shield. Southern Lake Huron 
is a mixture of water from outer Saginaw Bay, central Lake Huron, and runoff 
from the Ontario shore. Each of these three inputs to southern Lake Huron tends 
to dominate the southern basin in a different season. 
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CHAPTER FIVE 

SPATIAL DISTRIBUTION OF DISSOLVED NUTRIENTS: 

A COMPARISON BETWEEN 1980 AND 1974 

by David C. Rockwell 

Introduction 

Water quality in Lake Huron, Georgian Bay, and the North Channel has a 
large amount of spatial variability which is due, in part, to nearshore nutrient 
loadings and discharges from Lakes Michigan and Superior, Lake-wide average 
concentrations of most nutrients usually have large variances because of this 
spatial heterogeneity. Fixed geographical areas are compared in this chapter 
because water quality is evaluated in this manner by the public. Identification 
of environmental changes within fixed areas is a conservative measure of change 
because small magnitude changes may be missed due to the inherently larger vari- 
ances associated with fixed geographical area analysis relative to the homogene- 
ous water mass approach (see Chapters 6 and 8). By dividing Lake Huron, Geor- 
gian Bay, and the North Channel into 19 open lake and 8 nearshore areas, differ- 
ent areas in 1980 were compared with one another; and the 1980 areas were com- 
pared with their 1974 counterparts to detect gross changes between the two 
years. A statistical comparison of the surface layer or epilimnion (S/E layer) 
with the lower layer or hypolimnion (L/H layer) samples for the major nutrients, 
specific conductivity, pH, and alkalinity is contained in Appendix A. Statis- 
tically significant differences between S/E and L/H means were determined when 
the limnion 95% confidence intervals did not overlap. Appendix B contains 
statistical summaries of the nearshore areas and the open lake areas for the 

1974 and 1980 seasons. 
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1974 Segmentation of Lake Huron 

During the analysis of the 1974 Lake Huron surveillance data, a scheme 
of segmenting the lake was devised in an effort to reduce variance components 
associated with average nutrient concentrations. In this scheme, Lake Huron, 
Georgian Bay, and the North Channel were divided into 19 segments, each of which 
was considered a relatively "homogeneous water mass" (IJC 1976, 1977). Each 
segment was a compact geographic region with similar water quality (Fig. 5.1). 
Nearshore waters of Lake Huron, the North Channel, and Georgian Bay were divided 
into nine segments A to I (Table 5.1). Saginaw Bay segment H was not considered 
in the analysis of the 1980 data. Nearshore waters were defined in 1974 as a 
water column with a depth of 15 m or less and/or within 3 km from the nearest 
shore. Open waters were geographically segmented into 19 areas, with area 3 
covering parts of the North Channel and Georgian Bay (Table 5.2). Open waters 
were defined as greater than 3 km from shore and/or a water column deeper than 
15 m in depth. Seasons for 1974 were defined as spring (April, May, June); 
summer (July, August, September); and fall (October, November, December). 
Cruises 1 through 3 were considered spring, cruises 4 and 5 were summer, and 
cruise 6 was fall. Epilimnetic waters were defined as the upper 10 m of the 
water column in 1974. This definition was convenient in the nearshore areas. 
In deep water areas, the vertical sampling pattern seldom placed more than the 
1-m sample in this layer. 

1980 Segmentation of Lake Huron 

Comparison of 1980 results with 1974 results was conducted by dividing the 
nearshore and open water areas of Lake Huron into segments similar to the 1974 
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Figure 5.1 Segmentation for Lake Huron, the North Channel, and Georgian Bay 
(from IJC 1976 and 1977). 
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Table 5.1. 1974 near shore area acronyms, 1980 nearshore area acronyms, 
and 1980 monitoring sites assigned to area. 



1974 


1980 


Maximum 


_ , . , — „ 





Acronym 


Acronym 


Depth 


1980 Stations Assigned 


to Area 


Nearshore Area A 




17 


LH 1, 3 




Nearshore Area B 




11 


LH 4, 5, 10, 11 




Nearshore Area c'^ 




28 


LH 28, 30, 41 




Nearshore Area D 




134 


GB 101, 102, 103, 119, 


131, 132 


Nearshore Area E 




56.5 


GB 110, 115, 116, 123, 


125, 



126, 127, 138, 139, 141, 142 

Nearshore Area F 46 LH 71 , 72, 73, 74, 79, 81 , 84, 

85, 87, 88 

LH 34, 35, 36, 46, 47, 55 

LH 7, 13, 14, 17, 23, 25, 94 

LH 69 



Nearshore Area G 


21.5 


Nearshore Area I 


NS-I 23 


River 


river 13^ 



^Comparison of 1974 data with 1980 data required the exclusion of stations 
LH39 and 40 which appear to be included within nearshore area C of Figure 5-1. 
These stations were not included because they were more than 3 km from shore 
and had maximum depths greater than 15 m. 
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Table 5.2. 1974 open lake area acronyms, 1980 open lake area acronyms, 
and 1980 monitoring sites assigned to area. 



1974 


1980 


1980 Maximum 




Acronym 


Acronym 


Depth 


1980 Stations Assigned to Area 






North Channel 




Area 1 


WNC 


58 


68, 70-72, 75-76 


Area 2 


CNC 


76 


73-74, 77-86 


Area 3a 


ENC 


46 
Lake Huron 


87-89 


Area 4 


A4LH 


68 


60, 65-67 


Area 5 


A5LH 


53.5 


56, 62-64 


Area 6^ 


NBLH 


180 


27, 29, 31-33, 37-38, 43-45 
48, 51-54, 57, 61 


Area 7^ 


SBM 


62 


16, 18-21, 24, 26 


Area 8^ 


SBLH 


87 


6, 9, 12, 15, 90-93 


Area 9 




52 
Georgian Bay 


42, 49 


Area 10 




53 


136 


Area 11 




134 


130-135 


Area 12 




134 


118-120 


Area 13 




96 


101-102, 112 


Area 14 




86 


103-104, 106 


Area 15 




66 


105, 107-109 


Area 16 




56.5 


110, 115-116, 123 


Area 17 


CBGB 


95 


111, 113-114, 117-118, 121-122, 
124, 128-130, 137 


Area 18 




27 


125-127, 138-139 


Area 3 




41 


140-144 



^Stations 39-40 (Area 6), station 22 (Area 7), and station 8 (Area 8) are in 
depths greater than 15 m or more than 3 km from nearest shore but do not 
cluster with the associated areas. 
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scheme (Tables 5,1 and 5.2, Fig, 5,1). Ten of the nineteen segments from the 
1974 scheme were selected to compare with the 1980 surveillance results (Tables 
5.1 and 5.2). When epilimnetic and hypolimnetic chemical concentrations are 
presented for 1980 data in tables or graphs, epilimnetic refers to surface 
samples during unstratified periods and hypolimnetic refers to the lower layer 
during unstratified or isothermal periods. The lower layer is defined for each 
area and depends on the sample pattern. The epilimnion is approximated by 
taking all samples with a water temperature above a given temperature depending 
on the area and cruise involved. Similarly, the hypolimnion is approximated by 
taking all samples below another given temperature. These two temperatures were 
determined from the EBT temperature traces for all stations within an area for 
each cruise. S/E is used for surface/epilimnion data and L/H is used for lower 
layer /hypolimnion data. 

Nine of the segments selected for comparison were open water areas. 
The tenth area is nearshore area 1 where Saginaw Bay was expected to affect 
Lake Huron. Five of the open water areas were in Lake Huron. The northern and 
central basins of Lake Huron (NBLH) were combined and represented the largest 
segment in the open waters of Lake Huron. It was least directly affected by 
external inputs such as tributaries, flow from the other water bodies, and 
anthropogenic activities. The southern basin of Lake Huron (SBLH) was the 
southernmost extension of Lake Huron and was the area most likely affected by 
Saginaw Bay inputs. This basin was principally affected by the addition of 
dissolved materials derived from the mixing of Saginaw Bay waters containing 
anthropogenic inputs (Schelske and Roth 1973) with the northern and central 
basin waters. Saginaw Bay mouth (SBM) was represented by stations placed near 
and outside of the bay mouth to monitor effects of the bay on Lake Huron. 
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The nearshore zone - I (NS-I) represented the western Saginaw Bay nearshore area 
(stations LH 23, 25) and eastern Saginaw Bay nearshore (stations LH 17, 94) 
and the western nearshore along southern Lake Huron (stations LH 7, 13, 14). 
This segment was chosen to represent nearshore areas and to focus on the impact 
of Saginaw Bay upon Lake Huron. Two northern basin segments, A4LH and A5LH, 
represented areas where Lake Superior and Lake Michigan waters intermix with 
Lake Huron water. 

The central basin of Georgian Bay (CBGB) was the largest Georgian Bay seg- 
ment, both in surface area and in volume. Because the boundaries were set more 
than 3 km from shore, it was anticipated to be the least affected by tributaries 
and interactions with other water bodies. This segment was considered represen- 
tative of the chemistry and biology of Georgian Bay. 

The North Channel was divided into three open lake segments corresponding 
to the 1974 scheme, west North Channel (WNC), central North Channel (CNC), and 
east North Channel (ENC) (Table 5.1). The data for all three segments have been 
included because of the high variability of this water body. There was a large 
east to west gradient in the water chemistry for many of the 1980 cruises. 
The chemical composition of the North Channel waters was primarily affected by 
the outflow from Lake Superior. The chemical gradients from west to east were 
very pronounced; and, as a consequence, no single area in the North Channel was 
representative of the area as a whole (Warry 1978a). 

Northern Lake Huron waters were affected by the outflow of Lake Michigan 
through the Straits of Mackinaw and of Lake Superior through the St. Marys 
River. The northern areas of Lake Huron have concentrations of conservative 
elements which were approximately a mixture of 40% Lake Michigan and 60% Lake 
Superior water, the relative proportions of the two largest inputs to the 
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northeastern part of the lake (Schelske and Roth 1973). An additional source 
was the dissolved materials input by Canadian and United States tributaries. 
Canadian tributaries pass through silica-rich granites and gneisses of the 
Canadian Shield and carbonate-rich limestones and shales of the St. Lawrence 
platform (Warry 1978b). United States tributaries pass through carbonates, 
dolomite, shale, sandstone, limestone, siltstone, and unconsolidated sediments 
consisting of sand and gravel beds in glacial drift and post-glacial alluvium 
(US Geological Survey 1975). 

Results and Discussion 

Phosphorus — The three forms of phosphorus routinely measured during 
the surveillance program were total phosphorus (TP), total dissolved phosphorus 
(TDP), and soluble reactive phosphorus (SRP). From these analyses, particulate 
phosphorus and soluble organic phosphorus were computed. Annual mean total 
phosphorus (TP), total dissolved phosphorus (TDP), and soluble reactive phos- 
phorus (SRP) concentrations were higher in the L/H layer (Tables 5.3-5.24). 
In shallow water areas, vertical mixing minimized differences between the S/E 
and L/H layers. L/H concentrations were frequently greater in areas isolated 
from shore influences and in the open lake (Figs. 5.2-5.4). The greater concen- 
tration of TP in the L/H layer was probably due to algal utilization of phos- 
phorus in the S/E layer and the subsequent settling of the particulate forms 
(Hutchinson 1957, Paerl et al. 1975). Thus, statistically significant TP, TDP, 
and SRP concentration differences between S/E and L/H layers were observed 
primarily during the stratified period in the open lake (see Chapter 4). 
Thermal stratification was completed by cruise 4 and continued through cruise 6, 
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The seasonal ranges of TP, TDP, and SRP concentrations for Lake Huron and North 
Channel were 2.5-58.1, 1.4-7.6, and 0.1-5.2 yg/L, respectively. For Georgian 
Bay they were 2.3-15.3, 1.4-6.8, and 0.1-4.0 yg/L, respectively. 

The highest observed TP, TDP, and SRP occurred during the first two spring 
cruises when the water column was homo the rmous and thoroughly mixed by storms. 
Although low values for these forms of phosphorus were found throughout the 
season, SRP values were usually at a minimum during cruise 4 when water tempera- 
tures were at their maximum. This would be expected due to orthophosphate being 
considered the limiting nutrient (Beeton 1969) and to dissolved phosphorus 
limiting algal growth during the stratified period (Lin and Schelske 1981). 

Nearshore area B had the highest annual mean cruise TP concentrations 
(6.7 yg/L) of all areas with several cruise means exceeding 8 yg/L (8.3 yg/L 
spring - 8.5 yg/L fall). The St. Marys River station TP concentrations ranged 
from 6.6 to 10 yg/L with an annual cruise mean of 8.4 yg/L in the S/E samples 
(Table 5.23). 

TP mean cruise concentrations in the S/E layer during the six cruises 
ranged from 3.5 to 5.1 yg/L in the open lake areas NBLH, SBLH, CBGB, and SBM. 
Eighteen of the twenty four TP cruise means were between 4 and 5 yg/L and five 
were between 3.5 and 4.0 yg/L. Annual cruise means for TP concentrations were 
4.1, 4.3, 4.5, and 4.7 yg/L for CBGB (Table 5.5), SBLH (Table 5.21), SBM 
(Table 5.13), and NBLH (Table 5.3), respectively. Statistically significant 
differences between layers occurred for TP concentrations in 10 of 29 areas 
(Appendix A, Table A-1). The majority of these events occurred during cruise 4. 
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In the epiliranion, TP levels in the North Channel ranged from 5.7 to 
7.2 ug/L in WNC, 4.0 to 6.6 yg/L in CNC, and 4.7 to 7.5 yg/L in ENC (Tables 5.7, 
5.9, and 5.11, respectively). These ranges were similar to those found in 1974 
(Warry 1978b). 

TP levels in CBGB ranged from 3.1 to 8.1 yg/L in the epilimnion in 1974 
(Warry 1978a). The 1980 TP levels in CBGB ranged from 3.5 to 4.5 yg/L 
(Table 5.5). 

TDP and SRP concentrations in the North Channel and Georgian Bay during 
1980 had ranges similar to 1974 observations (Warry 1978a,b). Cruise mean TDP 
concentrations for 1980 were frequently between 2 to 3 yg/L in all areas of the 
lake and for both the S/E and L/H layers (Tables 5.3-5.24). There were only 
four occurrences when TDP concentration differences between these layers were 
statistically significant. These cases occurred in the northern and southern 
basins of Lake Huron during the stratified period (cruises 4 through 6) (Appen- 
dix A, Table A-2). 

Seasonal fluctuations were expected for TDP and SRP due to algal utiliza- 
tion (Hutchinson 1957, Beeton 1969, Lin and Schelske 1981). The changes which 
occurred in TDP and SRP concentrations were not similar to the changes observed 
for silica and nitrogen. Silica and nitrogen had a single annual minimum. 
Major reductions in TDP and SRP occurred every other cruise in some areas. 
The mean TDP maximum depletion was 30% with a 95% confidence interval about 
the mean from 25 to 35% (Table 5.25). The mean SRP maximum depletion was 56% 
with the 95% confidence interval about the mean from 51 to 61% (Table 5.25). 

SRP mean cruise concentrations seldom exceeded 1 yg/L in Lake Huron, and 
the majority of these occurrences were during cruises 1 and 5 (Tables 5.3-5.24). 
The SRP mean layer concentration differences were statistically significant in 
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areas A5LH, NBLH, SBM, SBLH, and CBGB, as well as area 9, area 11, and area 13 
during cruises 5 and 6 (Appendix A, Tables A-15 - A-30), 

The majority of maximum SRP and TDP mean cruise concentrations were 
observed during cruise 1, These levels reflect the nutrient concentrations 
prior to the early spring phytoplankton bloom. TDP mean cruise concentrations 
were highest in A4LH (3.8 yg/L), A5LH (4.2 yg/L), WNC (3.9 yg/L), and nearshore 
area F (4.4 yg/L). These areas were sampled shortly after ice-out conditions 
during cruise 1. Numerous small icebergs and floes were encountered in A4LH, 
A5LH, WNC, and nearshore area F. 

Dissolved Reactive Silica — Dissolved reactive silica (DRS) concentrations 
reflect seasonal uptake of this nutrient which is important for diatoms and 
silica flagellates to form their cell walls (Schelske and Stoermer 1971). 

Spring areal mean cruise DRS concentrations are summarized for the 
St. Marys River, North Channel, Georgian Bay, and Lake Huron. In the St. Marys 
River, station 69, the spring S/E layer's DRS concentration was 2.23 mg/L 
reflecting high DRS concentration of Lake Superior waters (Schelske and Roth 
1973). In 1980, in the North Channel nearshore, area F, the maximum spring 
areal mean cruise DRS concentration was 2.2 mg/L (Table 5.26). In 1974, high 
epilimnetic DRS concentrations were also found in the North Channel and ranged 
from 2.04 to 2.46 mg/L (Warry 1978b). Georgian Bay spring DRS areal mean cruise 
concentrations ranged from 1.2 to 1.6 mg/L with the highest levels found along 
the northern shoreline (area 18) and the lowest along the southern shoreline 
(area D). Typical Georgian Bay open lake areas had spring DRS areal mean cruise 
concentrations around 1.2 mg/L (areas 13 and 16). CBGB had spring DRS areal 
mean cruise concentrations of 1.2 mg/L in 1980 (Table 5.5). In 1974, Warry 
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(1978a) observed 1.4 mg/L spring DRS concentrations in CBGB in 1974. In 1980, 
Lake Huron spring DRS areal mean cruise concentrations ranged from 1.5 to 
1.6 mg/L in NBLH and SBLH. Nearshore spring DRS areal mean cruise concentra- 
tions in Lake Huron were similar in magnitude to the concentrations of adjacent 
open lake areas. 

Seasonal variation in DRS is often taken as a measure of nutrient utiliza- 
tion (Schelske and Roth 1973). The amount of depletion of DRS was estimated by 
comparing the maximum observed spring areal cruise mean of the S/E samples with 
the minimum observed summer areal cruise mean (Table 5.26). This measure of 
depletion is valid only in those open lake areas which are unaffected by upwell- 
ing and located beyond the influence of DRS point sources such as tributaries 
and other large lakes. The resulting average depletion in all areas was 35 + 2% 
with an average absolute depletion of 0.53 + 0.04 mg Si02/L. 

The general pattern in many of the geographical areas was DRS concentration 
depletion in the epilimnion and DRS concentration increases in the hypolimnion 
(Fig. 5.5). The DRS concentrations in the S/E and L/H layers are illustrated in 
Chapter 4, Figure 4.14. Six vertical (NBLH, SBLH, SBM, CBCG, A4LH, and A5LH) 
DRS traces show decreasing S/E concentrations during the early cruises followed 
by gradually increasing S/E concentrations during cruises 5 and/or 6. This 
resulted from a deepening of the epilimnion which mixed with hypolimnetic 
waters. The difference in concentrations between these layers became statisti- 
cally significant during cruises 4 and 5 in almost all areas (Appendix A, 
Table A-4). For most areas, statistically significant differences between S/E 
and L/H occurred during thermal stratification except for two cases. 



135 




as 



sn on 78 

1/2QIS OU 



0^ 



2 

2 V ~ 



2 * Si 

^ § 11 






■as 



i 




2? 


. 




O 


^ 


00 


8| 


00 




c 




•H 




J-i 


91 


3 




T3 




d 




CJ 



1/2CIS au 






ii u,« 



g^ ii 

°° Q if 



is( £ 



sS 




I 

Ui 

>- 

3 

'a 



is 

£ S 55 



"»> 



9 
i 

>- 

3 
fas 



SI on s-| 

/20IS OU 



onr 



1/ZQIS 3? 



a 

0) 

> 
1— ! 
O 
CO 
Cfi 
•H 

M-l 
O 

c 
o 






in 
in 

WD 



136 




ji 
i 

■ai 

t— 

4 

s 




zi 



>- 

S 
.1 



1/20IS 3U 



mf 



SI «n sn 

1/20IS DU 



o«r 



a „ 



Ij m ^ 



E sis 






8S 

ii 




■ZS 



>- 

s 

1.^ 



LU "" S _i 

S« =.= 

£ * Si 

GC o 5 2 




1/30IS 9U 



S!S 


/— N 


CJ 


0) 


d 


a 


'» 


c 


•H 


>- 


o 


2i 

•-t 




4 




i| 


LO 




LO 




0) 


s 


^ 




D 




CxO 



sn on ST 

1/20IS OU 



137 



|5 

is Sll 

S 2 11 

Vi it 



- 5 



So 



8§ 
11 



1/20IS OU 




•^1 



>- 

fas 



ai 



c 
>- 

CL 

d 

»5 



■s 



13 

^2 



^S Sll 



S8 
11 



IMI 



1/20IS Ou' 



8 

u 

la" 

>. 

3 

i 



'T3 
0) 

u 
o 



(D 
50 



1/20IS OU 



138 



The St. Marys River, station 69, had no observed depletion of DRS in the 
S/E waters. The WNC area showed the second smallest depletion (19%) of all the 
areas (Table 5.26). The WNC, like the St. Marys River station, received water 
from Lake Superior and did not show a large seasonal change. There was only one 
cruise where the DRS concentration differences between the epilimnion and the 
hypolimnion exceeded .3 mg/L even though the water column is relatively deep 
(58 m maximum depth). 

In the areas where statistically significant differences in DRS concentra- 
tions between the limnions were not noted, seasonal DRS depletion occurred. 
These depletions were as large as .6 mg/L in absolute magnitude, or a 55% deple- 
tion. Several reasons for the lack of statistical significance, in spite of the 
magnitude of the DRS differences, were the high variability of the observed DRS 
concentrations (areas 3 and 18), the shallow water column (nearshore areas A, B, 
and G), or a DRS concentration depletion in the entire water column between 
cruises (WNC). 

Station 78 showed complete depletion of DRS in the entire water column. 
This station represented most of the North Channel in 1980 (WNC and ENC). 
At this station, summer DRS concentrations decreased and did not return to the 
spring DRS concentrations by the sixth cruise (beginning of the fall overturn 
conditions). The DRS concentration at station 78 declined from 2.16 mg/L 
(cruise 1) to 1.61 mg/L (cruise 6). This complete depletion pattern was seen 
in areas of WNC, CNC, and nearshore areas C and I. In these areas, spring DRS 
concentrations cannot be restored by simply mixing in the hypolimnetic waters as 
in the deep water lake areas. Only by a resupply of DRS from the St. Marys 
River can this area be restored to spring DRS levels. Nearshore areas C and I 
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can be resupplied when the offshore waters are mixed into the nearshore zones 
during the winter period. 

Nitrate plus Nitrite — Dissolved nitrate + nitrite (DNN) concentrations 
decreased in the S/E layer as the water temperature increased (Fig. 5.6) • 
The difference between the observed spring maximum concentration and the 
observed summ**'" minimum concentration in the S/E layer is an indication of the 
utilization of the DNN, or nutrient depletion (Schelske 1975). Differences 
between S/E and L/H DNN concentrations were almost always statistically signifi- 
cant during stratified conditions in all areas of Lake Huron, Georgian Bay, and 
the North Channel (Appendix A, Table A-5). The L/H layer DNN concentrations 
were generally greater than those found in the S/E layer. This can be due to 
the algal uptake of DNN in the epilimnion where algal activity is enhanced in 
the photic zone (Hutchinson 1957, Wetzel 1975). 

Maximum DNN depletions observed in 1980 are shown in Table 5.27 for all 
areas. The comparison was made between the highest observed DNN concentration 
in the nearshore of 283 + 17 g/L and open lake of 276 + 10 g/L (usually 
cruise 1) with the lowest observed DNN concentration in all areas of 226 + 
4 g/L (which most frequently was observed on cruise 4). Maximum percent deple- 
tion in nearshore areas of DNN averaged 19.6% (95% C.I. = +4.2%) with areas 
A and B removed. There was an apparent source of DNN in nearshore areas 
A and B. These two areas had the highest spring concentrations in 1980. 
Elevated levels persisted from cruise 1 to cruise 3. In these areas, maximum 
DNN concentrations were found (station 11 - cruise 1 - 485 g/L, station 4 - 
cruise 2 - 735 g/L, and station 3 - cruise 3 - 765 g/L). Maximum percent 
depletion in open lake areas averaged 18.9% (95% C.I. = 4-2.7%). 
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In 1974, Warry (1978a) reported a CBGB areal epilimnetic raaximuni DNN con- 
centration of 260 yg/L. This compares to a 277 yg/L DNN concentration in 1980. 
For a North Channel epilimnetic spring maximum, Warry (1978b) reported DNN 
concentrations of 300 yg/L (WNC), 270 yg/L (CNC), and 280 yg/L (ENC). These 
compare with 296 yg/L (WNC), 300 yg/L (CNC), and 229 yg/L (ENC) in 1980. 

Ammonia — Relatively stable and homogeneous dissolved ammonia (DA) con- 
centrations were observed throughout the water column in NBLH and SBLH during 
the first and second cruises (Tables 5.3-5.24, Fig. 5.7). Samples from these 
two areas ranged from 1.6 to 2.2 yg/L for cruises 1 and 2. Differences in DA 
concentrations between S/E and L/H layers became statistically significant 
during cruise 4 (Appendix A, Table A-6) in 14 of the 27 areas. The concentra- 
tions observed in areal hypolimnetic DA were two to four times the observed 
epilimnion DA concentrations found in each area during cruise 3 and/or 4. 
Maximum DA concentrations in the hypolimnetic layer ranged from 6.6 yg/L CR3-CNC 
to 11.4 yg/L CR4-SBLH, with six areas between 9-12 yg/L (ENC, SBM, SBLH, area 
12, area 15, and CBGB). Warry (1978a) reported higher hypolimnetic DA concen- 
trations in the North Channel, ranging from 15-16 yg/L in CNC and ENC in June. 
In these areas, the hypolimnetic DA concentrations found in June were two to 
five times those of other cruises (Fig. 5.7). 

Higher DA concentrations were observed for cruises 3 and 4 in most of the 
hypolimnetic traces for these areas. North Channel DA concentrations in WNC 
were high during cruise 1 (>12 yg/L). Warry (1978a) reported levels as high as 
24.9 to 32.3 yg/L during May in WNC. Near shore area A was the only place where 
DA concentrations exceeded those found in WNC (15.5 yg/L - cruise 4 hypolim- 
nion) . 
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In CNC (cruise 4), SBM (cruise 5), SBLH (cruise 5), and open lake area 13 
(cruise 5) the areal DA concentrations for the metalimnion showed increased 
levels when compared to the epilimnion and hypolimnion. This DA pattern was 
also found in Lake Michigan's southern basin in 1976 and 1977 (Rockwell et al, 
1980). Increased DA concentrations are often attributed to increased zooplank- 
ton excretion or dieoff and decay of algal forms (Evans 1983) • 

Specific Conductivity — In 1974, specific conductances ranging from 133 
to 191 ymho/cm in WNC, 155 to 198 ymho/cm in CNC, and 164 to 181 ymho/cm in ENC 
were reported by Warry (1978a), In 1980 the corresponding ranges were 125.9 
to 178 ymho/cm in WNC, 155 to 194,8 ymho/cm in CNC, and 157.7 to 172.6 ymho/cm 
in ENC (Tables 5.3-5.24). 

In 1974, specific conductances ranging from 182 to 194 ymho/cm in CBGB 
epilimnetic layer were reported by Warry (1978b). In 1980 the corresponding 
ranges in CBGB were 178.8 to 193.7 ymho/cm (Tables 5.3-5.24). 

Statistically significant conductivity differences between S/E and L/H 
layers developed during most of the stratified cruises even though this variable 
is not expected to be influenced by temperature or biota (Fig. 5.8) (Hutchinson 
1957, Wetzel 1975). In addition, no observations of calcium carbonate precipi- 
tation, "milky waters," were reported. Cruise 4 had sixteen of the twenty seven 
areas showing statistically significant differences (Appendix A, Table A-7). 

Chloride and Sulfate — During 1980, chloride and sulfate concentrations in 
Lake Huron had spatial distributions which were similar to the specific con- 
ductance distribution. In the St. Marys River, chloride concentrations ranged 
from 1.4 to 6.3 mg/L and sulfate concentrations ranged from 3.2 to 16.7 mg/L. 
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Annual average A5LH chloride and sulfate concentrations were 6.3 mg/L and 
15.7 mg/L in the S/E layer and 5.8 mg/L and 15.7 rag/L in the L/H layer, respec- 
tively (Tables 5.19-5.20). These values represent a mixing of higher concentra- 
tion Lake Michigan waters (7.7 + .1 mg/L and 21.1 + .4 mg/L chloride and sul- 
fate, respectively) (Rockwell et al. 1980) with lower concentration Lake Huron 
waters where NBLH area concentrations of 5.4 mg/L chloride and 15.3-15.5 mg/L 
sulfate were found in 1980 (Tables 5.3-5.4). 

In 1974, Warry (1978a) found North Channel chloride and sulfate concentra- 
tions to increase generally from west to east with low levels of 2.7-5.4 mg/L 
chloride and 8.0-15.1 mg/L sulfate in the WNC, intermediate levels of 4.4- 
5.5 mg/L chloride and 13.7-15.5 mg/L sulfate in CNC, and higher levels of 4.4- 
4.8 mg/L chloride and 14.7-15.6 mg/L sulfate in the ENC. In 1980 similar ranges 
were found in each of these areas (Tables 5.3-5.24). 

In 1974, for the central basin of Georgian Bay, Warry (1978b) found S/E 
chloride and sulfate concentrations of 4.7-4.9 mg/L chloride and 15.2-15.9 mg/L 
sulfate. In 1980, the S/E chloride ranges were 4.7-5.1 mg/L and sulfate ranged 
from 14.8 to 15.7 mg/L. In CBGB, chloride and sulfate levels showed no change 
from 1974 to 1980. 

Figure 5.9 shows area sulfate concentrations in S/E and L/H layers. 
The L/H layer in the North Channel, St. Marys River, and CBGB tended to have 
higher sulfate concentrations than the S/E layer. This appears to be due to 
intrusion of water with higher dissolved solids from Lake Huron into the hypo- 
limnion layer in these areas. 

A discharge event was observed in the St. Marys River during cruise 4. 
Elevated levels of chloride (Fig. 5.10) (200% above average), sulfate (about 
300% above average), and specific conductance (about 70% above average) were 
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observed. These measurements were associated with the lowest pH measurement 
of the 1980 season at this site. Similar chemical responses have been observed 
in areas where pickling liquor is discharged during steel manufacturing. 

Alkalinity — Alkalinity differences between the S/E and L/H layers were 
significantly different at the 95% confidence level at one or more times during 
the stratified period in 22 of the 27 areas (Appendix A, Table A-8). This 
occurred most frequently during cruise 4, The increased levels in the hypolim- 
nion are associated with the accumulation of CO2 (Hutchinson 1957). 

Large variations in areal alkalinity concentrations are due to the differ- 
ent sources of water rather than due to seasonal variations. The lowest alka- 
linity measurements were found at the St. Marys River station with annual mean 
levels of 43.2 mg CaC03/L in the S/E layer and 47.7 mg CaC03/L in the L/H layer. 
The hypolimnetic alkalinity value of 71 mg CaC03/L was associated with unusually 
high sulfate, chloride, and conductivity values (Fig. 5.11). Lake Superior 
waters had alkalinity levels between 42-45 mg CaC03/L (Schelske and Roth 1973). 
The North Channel's lower alkalinity values reflect Lake Superior source waters 
and drainage from the Canadian Shield. The highest alkalinity values were found 
in the open lake area (A5LH) with values ranging from 85 to 92 mg CaC03/L with 
an annual mean cruise average of 88.3 mg (Tables 5.3-5.24). A5LH is most 
directly affected by Lake Michigan water. Northern basin Lake Michigan alka- 
linity was reported around 110 mg CaC03/L (Rockwell £t al . 1980). 

Nearshore areas G and I had similar high alkalinity levels. Nearshore 
area G is also influenced by Lake Michigan water. Nearshore area I is influ- 
enced by Saginaw Bay outputs. Nearshore area G mean annual alkalinity concen- 
trations were 80.5 mg CaC03/L in the S/E layer and 79.7 mg CaC03/L in the L/H 
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layer. Nearshore area I mean annual alkalinity concentrations were 80.3 mg 
CaC03/L in the S/E layer and 79.8 mg CaC03/L in the L/H layer (Tables 5.3-5.24). 
The open lake (NBLH, SBM, and SBLH) mean cruise alkalinity concentrations in the 
S/E and L/H layers ranged between 74 and 81 mg CaC03/L with SBM values slightly 
higher (Fig. 5.3, Tables 5.3-5.24). 

Central basin Georgian Bay (CBGB) mean cruise alkalinity concentrations 
ranged between 69 and 75 mg CaC03/L in the S/E and L/H layers (Fig. 5.3, 
Tables 5.3-5.24). L/H alkalinity values tended to be higher in the deep water 
areas and the North Channel. In areas influenced by high alkalinity source 
waters (A5LH, NS-I, SBM), the reverse was generally true (Fig. 5.3). Warry 
(1978b) reported alkalinity between 51.5 and 73.8 mg CaC03/L (WNC), between 
57.9-75.7 mg CaC03/L (CNC) , and between 60.4-64.8 mg CaC03/L (ENC). Warry 
(1978a) reported alkalinity between 68.2-72.6 in the epilimnion of CBGB. 
These levels were similar to those found in 1980. 

pH — The pH values were statistically different between the S/E and L/H 
layers during several cruises in 1980 (Fig. 5.12). Cruise 4 accounted for the 
majority of these occurrences (Appendix A, Table A-9). The pH values were 
usually above 7.8 for cruise means in the S/E and L/H layers in all areas. 
Slightly lower S/E pH values were observed in WNC during cruises 1 and 2. 
This may reflect acid spring runoff from the Canadian Shield. The pH values 
in the L/H layer tended to be more acidic (lower pH values) in some areas 
(Fig. 5.4). This pattern was observed in Lake Michigan (Rockwell et al. 1980) 
and is probably due to higher CO2 concentrations in the colder L/H layer (Hutch- 
inson 1957). 
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In 1974, Warry (1978a) reported mean pH values between 7.8 and 8,2 with 
a maximum of 8,4 and a minimum of 7.6 in the North Channel. In CBGB, Warry 
(1978b) reported mean pH to range between 8.0 and 8.2 with a high of 8.6. Areal 
pH extreme values reported in 1980 were slightly lower in CBGB with a hypolim- 
nion low of 7.85, but overall, pH levels were similar. 

Dissolved Oxygen — Dissolved oxygen (DO) was seldom depleted in 1980, and 
the overall average percentage saturation was 99% or greater on every cruise for 
all areas of Lake Huron and Georgian Bay (Table 5.28). The lowest observed DO 
saturation (70%) occurred during cruise 2 at Lake Huron station 27 in the bottom 
sample. Figure 5.13 shows DO concentrations in the S/E and L/H layers in ten 
selected areas. 

Major Metals — The alkali metals, potassium and sodium, and the alkaline 
earth metals, calcium and magnesium, were measured during cruise 1 in all areas 
except part of CNC and ENC where ice prevented operations. Table 5.29 contains 
the surface and lower layer concentrations for these metals. Concentrations of 
these metals are lower in Lake Superior water than in Lake Huron (Schelske and 
Roth 1973). As a result, the St. Marys River had the lowest concentration of 
these parameters. The average concentrations of all samples at the St. Marys 
River station were .70 mg/L potassium, 1.5 mg/L sodium, 2.8 mg/L magnesium, and 
13.4 mg/L calcium. Weiler and Chawla (1969) reported Lake Superior calcium 
concentration at 13.2 mg/L and Schelske and Roth (1973) reported Lake Superior 
calcium concentration at 14.5+1.5 mg/L. North Channel areas had metal concen- 
trations intermediate between those of the St. Marys River and those of Lake 
Huron. The average value of all samples in area CNC was .93 mg/L potassium, 
3.0 mg/L sodium, 22.6 mg/L calcium, and 6.0 mg/L magnesium. In most areas, 
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14.1-14.4 


102-106 




GB 


104 


95-109 


SBLH 




14.4 


14.2-14.7 


97-108 




LH 


103 


79-113 


NBLH 




13.6 


13.2-14.1 


82-104 










Cruise 


2 
















CBGB 




13.5 


13.3-13.7 


98-105 




GB 


104 


98-112 


SBLH 




14.0 


13.8-14.2 


102-110 




LH 


105 


70-118 


NBLH 




13.6 


13.0-14.1 


70-116 










Cruise 


3 
















CBGB 




13.2 


13.1-13.5 


103-111 




GB 


106 


100-115 


SBLH 




13.7 


13.6-13.8 


104-111 




LH 


106 


98-116 


NBLH 




13.7 


13.6-13.7 


102-108 










Cruise 


4 
















CBGB 




12.2 


11.8-12.5 


89-107 




GB 


100 


86-115 


SBLH 




11.7 


11.0-12.4 


85-100 




LH 


103 


85-121 


NBLH 




12.5 


12.3-12.7 


93-102 










Cruise 


5 
















CBGB 




11.3 


11.0-11.6 


85-119 




GB 


100 


85-122 


SBLH 




10.9 


10.4-11.3 


79-133 




LH 


102 


79-133 


NBLH 




11.9 


11.7-12.2 


88-121 










Cruise 


6 
















CBGB 




11.8 


11.7-12.0 


95-107 




GB 


100 


85-115 


SBLH 




11.0 


10.9-11.1 


89-105 




LH 


99 


89-111 


NBLH 




12.0 


11.9-12.2 


90-111 











169 



tS 



ia 



CO 

t 

5 



8S 
11 





tZi 



>- 



1/OU 



1' 



zl 



s 




1/OU 



oc g -- 
en H 8S 



M 



iSS 



s 



■SI 



tS 



I 






o 

00 



to 
c 

•H 

T3 

C 
(U 

>, 

X 

o 

0) 
> 

T—l 

o 

CO 
CO 

»H 

O 

C 

O 



> 

cu 



cu 

13 
GO 



n/ou 



i/ou 



170 



i^ 






r 

i 




t3I 



I 

>■ 



»: 



^s Si 



8i 



<n <n 




KS 



s 



■ s 



V3U 



1/3U 



Si 










"H: 



Ttii 



LlI 



s 

8y 




1/OU 



2i 


• 


h- 

u 


^3 


^ 


Q) 


«^ 





4J 


>- 


c 


1^ 


o 


3? 


CJ 


•-» 




^ 




a| 


CO 


M 




LO 


91 






3 




W) 



n/ou 



171 






So 



S? o 8 8 



9i 






l/QU 



S 

fas 



«Z 



3 S(S 

O in in 

uj — — 

5 ^ ^g 



§ 




Jl 



f^S 



>- 
i 



— P 22 

§i 11 

i° si 




25 



s 



•« s 

1/OU 



•s 



rH 

o 



0) 
•H 



"I/OU 



172 



CQ CO 



c 


(U 


CU 


Xi 


w 




(V 




M 


c 


ex 






V 


0) 


sz 




■u 


OJ 




5 


? 


CO 


n 






3 (U 



■• c 

3 VM 



- a 



:J 



1) 

3 
-5 



a> CO 

SI Qi 

H CD 

0) 



^3 



8 



m -« 


r^ • 


STo 


oi O '-^ 


l/> -^ ^s 


a* 1 '-v 


ON O ^ 


a\ \ m 


O 1 o 


^ :2 ^ 


• 1 f>J 


• 9\ ^ 


• sO -< 


• -^ -^ 


O 30 >-^ 


O r^ >-/ 


— « OS >-• 


O 00 w 


-* O >-/ 




OS 




• 




o 


• 


O 


o 


.^ 


— ^ 


o 


— * 




_ii« 



!G 9 ^ - 

30 I '^ OS I /-s 

• 30 r^o • 30 vO 

O P^ w O 30 >^ 



30 -< 
O^ I '^ 

• 30 ^O 
O 30 w 



a -^^ 

3 .-I 

CO N^ 



o 


3-^ 


Xi 


•H hJ 


1 


09 ^^ 




Sf 


0) 


00^^ 



3 

•H 


J 


O 


'i*^ 


CO 


f 



a 

00 >^ 

00 bO 

(0 a 



§3 

cn N^ 



a 

CO >^ 

c a 

bOv-/ 






CO g- 



-4 


• 


00 "-^ 


sO ^ 


00 — N 


• /— ^ 


-• vO ^^ 


"* -: ^ 


^ • <n 


>o • o 


■* m 


-! ' ^ 


• <n -< 


• .-n -♦ 


. CO z; 


1 Vi^ 


en ON ^-^ 


en 1 >^ 


m 1 x-' 


en 1 >-• 


■n 


!>• 


O 




en 



-^ 


en 


O ^^ 


en .'^ 


• /"^ 


ON -• ^ 


OS • o 


CM • m 


r^ m 


• 1 CM 


• rv -^ 


• r>. -^ 




\o aa ^~^ 


vfi 1 ^^ 


r>> 1 >»/ 


O 


r- 


30 


O 



r>. I >.^ 



O 



ON 
ON • /'-S 

-* I w 
O 



o • ^ 

• en \0 

en I ^r 



o • ^ 



^ . ^ 


NO • 


. en r^ 




en 1 w 




m 


vO 



CM 

.n — • . 

ON I 
O sO ' 



O ' 
e>4 I > 



m 



in 


O 00 


en . .^s 


• rsi 


• r>» r^ 


00 1 


r^ 1 >_• 


CM O 


CM 




• 


00 



nO 


, 


nO 


CM CM 


00 CN 


ON CM /^ 


• 1 ^ 


• ^ ^ 


• 1 O 


-^ 00 m 


•^ 1 CM 


in en -^ 


CM • ^-/ 


CM !>. s-' 


CM • N-/ 



CM 



CM CM 



OS CM /-N 

• I O 

vO '^ -< OS CM r>. 

CM • >-x w^ • >,^ 

CM .^ 



>^ 


00 


m 


m CM 


-* CM 


en — « 


• 1 ^v 


• 1 ^ 


• 1 


CM OS vO 


r** r^ r>. 


m <n 


CM . ^ 


CM . >^ 




O 


vO 


en 



30 

• I 



vO 
■M CM . 
• I 



g 


OS 
OS 


o 

CM 


s 


Qtf 


>* • 




8-? 


3 


OO O ^ 


o -• ^ 




• 1 -^ 


• 1 -3- 


• 1 ^ 




O OS ^^ 


-» vO >-^ 


-^ in s-*' 


s 


vO 


X 


ON 



O ON 

O OS 

>* • vO • 

OS ^-4 ^-s ON O >-> 

• I 00 • I 1^ 

O r»>» v-^ O CM s.^ 

00 OS 



3 ° 



m en 

OS OS 

O • en • 

00 O -^^ OS O ''^ 

• I \0 • I en 

O r«» >i^ O OS >*• 

m Qo 



CM • ^^ o • 

O — * 00 P^ o ■ 

•1-4 • I 

--^ in -w ' o o ' 

OS i>^ 



>* 




< 


— ^ 


OQ 


^ 




OS 


Z 


00 . ^ 


< 


CO O CM 


M 


• 1 -* 


a 


O CM w 


c^ 


00 


o 






C5 


u 





en • />-\ 

• m ^ 

en I s-» 

00 

CM 



en 


sO ^ 


r<. 


00 


vO • >-N 


en . >3. 


m • ^^ 


m • ^-s 


^ "f ^ 


• en -• 


• en 00 


• en r^ 


en 1 >»• 


en 1 ^ 


en 1 v^ 


en 1 v-/ 


00 


»^ 


CM 


CM 



J- . ^ 


o • ^ 


— « . ^ 


^ • ^^ 


—4 • ^N 


« 'T* ^ 


• OS ^ 


• r>» — * 


• r^ 00 


• r^ p^ 


1 N-' 


00 1 w 


r^ 1 N-^ 


r>. 1 >-• 


r* 1 s^ 


sj- 


en 


O 


CM 


00 



CM vO 



OS . ^s 

• in vo 
en I w 



' en en 



O • "^ 
• vO en 



en CM 



-d- . 



' en — • 

I I ^^ 

X 



en 



sj- ^^ o c^ 



vO I 



CM CM 

• 1 ^ 


vO en 

• 1 ^-^ 


sr >^ -* 

CM • «w 


r^ O ^ 

CM • N-/ 


CM 


S 



00 m ,^ 

• • • 

^ r^ r^ 

>^ CM ^-N O CM m CM 

• I ^ .],*-> . I >-N 

vOOO-H r>»inoo sOf^r^ 

eM«N^ CM.s^ CM»>-/ 

m NO m 

<^ CM CM 



— • 


,~^ 


u 


_ 


en 


sr 


> 


sO 


^ 


en 


8 


CM 


-* <N 


fN. CM 




en CM /-N 


^ ' ^ 


• 1 ^N 


Ex] 


• 1 — • 


vO CM so 


CM OS en 


U 


r^ O CM 



en in r>» 




en >^ in 


m 


1 1 - 


m 


-. m ^ 




<n >* m 


* 


• • 1 


CM 


OS 00 -• 


—4 c^ 


CM en m 


• ON 


•> 1 


ON 1 


r^ r^ 00 


• o 


CM en ^ 


VO ON 



>-^ CM • s-/ M 



vO 
00 

I 



-a- — • 
• I ' 

en ^ . 



o so 



00 m 
vO p^ 



OS 
00 

I 



—4 


en 




-^ 






»* 


CM O 




CM en 


— 


— . -H 


f^ 


1 1 

-• 30 


— 


CM CM 










—4 


30 -3- 


-^ 


— ♦ CM 



32 



PQ 



s 



3 






173 



Georgian Bay metal concentration levels were 1 to 2 mg/L lower for sodium, 
magnesium, and calcium than for the rest of Lake Huron. CBGB metal levels were 
.91 mg/L potassium, 2.9 mg/L sodium, 6.8 mg/L magnesium, and 25.5 mg/L calcium. 
In Lake Huron, the highest metal levels were in areas influenced by Lake Michi- 
gan (A5LH) and Saginaw Bay (SBM, NS-I) waters. Lake Michigan waters contain 
1.1 mg/L potassium, 4.8 mg/L sodium, 10.8 mg/L magnesium, and 34.9 mg/L calcium 
(Rockwell et al. 1980). NBLH is representative of open lake metal concentra- 
tions. For all NBLH samples the concentrations were 1.0 mg/L potassium, 
3.3 mg/L sodium, 7.0 mg/L magnesium, and 26.1 mg/L calcium. 

Comparison of 1980 to 1974 Segmentation Results — Comparison of nearshore 
areas A, B, C, D, E, F, G, and I results from 1980 with those of 1974 shows 
chemical changes suggestive of improved water quality (Table 5.30) between these 
2 years. Secchi disk measurements in 1980 were greater in six of eight areas. 
Total phosphorus concentrations were lower in 1980 in seven areas. Total nitro- 
gen levels were higher in five areas, with one area showing no change. 
Dissolved reactive silica concentrations were higher in six areas, with one area 
showing no change. Only two areas had decreased chlorophyll a_ measurements. 
Dissolved substances changes were mixed, with conductivity and sulfate showing 
increased concentration levels in a majority of the areas and chloride showing a 
decrease in a majority of the areas. (Data tables are found in Appendix B.) 

These observations should be used with caution because they are based on an 
annual comparison. The results can be biased between the years due to different 
cruise frequencies and areal station patterns, especially in the nearshore 
areas. There appears to be a general reduction in nearshore phosphorus levels, 
with some Canadian areas having large TP reductions between 1974 and 1980. 
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Table 5.30. Comparison of 1974 and 1980 epilimnetic/surface waters 
(1-10 m) in nearshore areas for selected parameters. 
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In 1974, one season was omitted in areas A, B, C, D, G, and I. For some of 
these areas, decreases in total phosphorus are quite large: areas A (approx. 
21.5 ppb) , B (approx. 15.5 ppb) , C (approx. 5 ppb), and F (approx. 7.5 ppb). 
There appears to be a general reduction in nearshore phosphorus levels, with 
some Canadian areas having large TP reductions between 1974 and 1980. 

Comparison of 1980 open lake segments with the 1974 areas for five para- 
meters (Table 5.31) also shows a decrease in total phosphorus concentrations. 
The 1980 results suggested that total phosphorus levels decreased in 18 of the 
19 areas, and dissolved nitrate + nitrite increased in 14 of the 19 areas. 
These changes were consistent with improved trophic status. The dissolved 
reactive silica (DRS) levels were lower in 13 of the 19 areas. The conservative 
ions were evenly split between increased levels, no change, and decreased 
levels. 

Comparison of areas A4LH, A5LH, NBLH, SBM, SBLH, and CBGB for nutrient 
depletion in 1980 with that found in 1974 showed about 50% less depletion in 
1980 for DRS and DNN. This is both statistically and environmentally signifi- 
cant. In these areas, DRS depletion in 1980 was 41+2% while in 1974 these areas 
had an average depletion of 65+4% (IJC 1976). DNN depletion in 1980 was 22+2%, 
while in 1974 the areas reported had an average depletion of 43+4% (IJC 1976). 
Depletion of silica and nitrate in the epilimnion increases with increasing 
degrees of eutrophication (Schelske 1975). Between 1974 and 1980, the opposite 
effect occurred, suggesting improved trophic conditions. 

Comparison of the epilimnion 1974 to 1980 areal analysis results with the 
1968 to 1980 year time series analysis (Chapter 8 - Table 8.7) provides some 
useful inferences. For the areal results (Table 5.31), chloride and conduc- 
tivity were about evenly split between increases and decreases. The long-term 
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Table 5.31. Comparison of 1974 and 1980 epilimnetic/surface (1-10 m) 
in open lake areas for selected parameters. 
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trend analysis showed increasing and decreasing trends within the homogeneous 
water mass, depending on the season. The areal results taken together suggested 
decreasing total phosphorus and dissolved reactive silica trends and increasing 
nitrate-nitrogen trends. Except for total phosphorus, these results are sup- 
ported by the long-terra trend analysis. 

This suggests that an individual area result may be indicative of long-term 
changes. Further, if an individual area result is supported by a majority of 
areas, then more sophisticated statistical analysis (see Chapter 8) must confirm 
the direction of change. 

Total Phosphorus-Ice Cover Relationship — What may be the cause (or 
causes) of these changes? Because the collection of changes suggests improve- 
ments in water quality, the apparent improvements in water quality may not be 
just normal environmental variation in an oligotrophic system. Lower total 
phosphorus levels have also been observed in Lakes Michigan (Bowden et al. 1981) 
and Ontario (Kwiatkowski 1982) during this time period. The phosphorus 
decreases in these lakes, as well as in Lake Huron, cannot be explained solely 
by cultural controls on phosphorus loadings. Thus, an additional mechanism 
which can influence the entire Great Lakes basin is suggested. 

The reduced total phosphorus levels found during 1980 follow a series of 
colder winters. Figure 5.14 shows several measures of the extent of the cold 
winters in the late 1970s and the spring cruise total phosphorus concentrations. 
Cumulative days of percent ice cover in the winters of 1976-77 through 1978-79 
were more than twice that which occurred in five of the six winters during 
1971-76. Maximum areal ice covers occurred in the winters of 1976-77 and 
1978-79. A more complete settling of particulate matter occurred in the winters 
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of 1977-79 than otherwise would have occurred. Wind stress effects would have 
been reduced by the more extensive ice covers. A statistical relationship 
between maximum ice cover and spring total phosphorus concentration levels has 
been shown to exist in Lake Michigan, Lake Ontario, central and eastern Lake 
Erie, as well as Lake Huron (Rockwell 1981). This study showed that ambient 
phosphorus levels tended to decrease after winters with extensive ice covers 
and tended to increase after mild winters. The decreases were larger than 
the increases with a resulting net reduction in total phosphorus. 

Rodgers and Salisbury (1981) used the hypothesis of a time variable setting 
velocity in their model study of the effect of winter ice cover on the water 
quality in Lake Michigan. They found that they could simulate the reduction 
in total phosphorus concentration observed in Lake Michigan between 1976 and 
1977 (Rockwell et al. 1980) by assuming that the rate of particle removal from 
the water column was dependent on the extent of lake ice cover. 

Lesht (1984a) examined the effects of this assumption on a long-term 
simulation of phosphorus concentration in Lake Michigan using a simple mass 
balance model, forced with historical load estimates (Chapra 1977), and found 
that the time variable settling velocity hypothesis is consistent with observed 
phosphorus concentrations . 

Lesht (1984b) applied a time-dependent solution to a multi-segment model 
of the Great Lakes basin. The time-dependent solution is used to evaluate year- 
to-year variations in the system response to time-variable forcing. The average 
annual total phosphorus in Lake Huron decreased over the period 1976 to 1980, 
with further total phosphorus decreases projected through 1982. Even greater 
decreases were projected in Lake Huron when Lake Michigan ice-cover-dependent 
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setting velocities were applied, because this further reduced the loads from 
Lake Michigan into Lake Huron. 

Direct measurenient of ice cover effects on setting rates is needed to 
confirm the validity of the hypothesis that setting rates are dependent on ice 
cover extents. Annual surveillance of Lake Huron prior to stratification and 
after fall overturn would be needed to quantify the extent of ice-cover-induced 
effects on natural environmental variation of water quality. Annual surveil- 
lance would be useful for following total phosphorus trends and confirming the 
simple mass balance model's projected decreases based on loads and ice cover 
effects. 

Conclusions 

For periods of thermal stratification, a number of chemical differences 
between the epilimnion and hypolimnion were noted. TP, TDP, SRP, DRS, DNN, DA, 
and alkalinity were higher in the hypolimnion; whereas, conductivity and pH were 
higher in the epilimnion. The mean TDP, SRP, DRS, and DNN maximum depletions 
for 1980 were 30%, 56%, 35%, and 18.9%, respectively. During 1980, the average 
dissolved oxygen saturation was 99%. 

For nearshore areas, 1980 Secchi disk depths, TN, DRS, sulfate, and con- 
ductivity were greater than those reported for 1974. Conversely, the 1980 TP 
and chloride concentrations were less than those reported for 1974. 

Within the open lake, 1980 TP and DRS were less and DNN was greater than 
those reported for 1974. The decrease in TP followed a series of relatively 
colder winters. There was a 50% decrease in mean DRS and DNN maximum depletions 
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in 1980 relative to those of 1974, Thus Lake Huron's trophic status appears 
to have improved between 1974 and 1980, 

Regions of Lake Huron impacted by Lake Superior water had relatively low 
alkalinity, sodium, calcium, magnesium, and potassium concentrations. Con- 
versely, regions impacted by Lake Michigan water had relatively high alkalinity, 
sodium, calcium, magnesium, and potassium concentrations. Regions measurably 
impacted by Lake Superior were confined primarily to the North Channel and Lake 
Huron waters adjacent to the western North Channel. Regions impacted by Lake 
Michigan were confined primarily to northern Lake Huron adjacent to the Straits 
of Mackinac and a region along the western shore of northern Lake Huron. 
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CHAPTER SIX 

HOMOGENEOUS WATER MASSES - STATISTICAL APPROACH 

by Russell Moll 

Lake Huron is generally considered oligotrophia (IJC 1977, Schelske et al> 
1974). Low levels of dissolved nutrients and chlorophyll are found throughout 
most of the open lake. But, extensive lake-wide limnological surveys have shown 
large areal differences in many variables. Furthermore, the areal distribution 
of nutrients and chlorophyll changes considerably across different seasons. 
Moll et al. (1976) used a multivariate statistical approach to identify regions 
of homogeneous water quality in the Straits of Mackinac region of Lake Huron. 
Berry (1980) used somewhat similar statistical techniques to define regions 
of homogeneous water quality in southern Lake Huron during the spring, summer, 
and fall of 1974. These two studies confirmed the same general hypothesis: 
although Lake Huron is an oligotrophia Great Lake, many regions of distinctly 
different water masses can be identified. These homogeneous water masses do not 
have a static areal location, but rather change position across seasons. 

The motivation behind identifying homogeneous water masses in Lake Huron 
is twofold. First, the identification of the homogeneous water masses and their 
locations can broaden an understanding of how different water masses mix. 
By locating gradients between homogeneous water masses, inferences about water 
mass and nutrient exchanges can be developed. Second, seasonal or long-term 
trends of water quality in Lake Huron cannot be identified unless homogeneous 
water masses are compared across time. 

The approach used in this chapter is to consider the identification 
of homogeneous water masses from a multivariate perspective. The areal 
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distribution of variables considered one at a time was presented in Chapter 
Four. Lake Huron was then divided into many segments, and variables within 
individual segments were considered one at a time. The results from this 
segmented approach are presented in Chapter Five. Inferences about overall Lake 
Huron water mass distribution are difficult when multiple plots are considered. 
The multivariate approach has the benefit of identifying trends in the data 
which are reinforced by the fact that the trends occur in several variables as 
opposed to one variable (Pielou 1977). 

Analytical Statistical Techniques 

The statistical problem of developing homogeneous water masses comes 
under the general topic of statistical classification (Sneath and Sokal 1973). 
A variety of analytical methods is available to perform the classification, 
but they all require some degree of subjective interpretation. No statistical 
method yields a clear and precise separation between homogeneous water masses. 
As a consequence of this imperfect separation, several statistical techniques 
were used and the results compared for consistency. The underlying assumption 
of this approach was that, if true patterns exist in the areal distribution of 
water masses, the patterns should emerge with several analytical methods. 
Cluster analysis, principal components analysis, and discriminant functions 
analysis were the three analytical procedures used to define the homogeneous 
water masses. The analyses were conducted on a cruise-by-cruise basis for two 
reasons: 1) homogeneous water masses were not expected to have the same areal 
distribution among the six cruises, and 2) the large size of the database 
prevented analysis of all cruises at once. The analytical procedures are 
briefly described below. 
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Cluster Analysis — Cluster analysis attempts to develop hierarchial group- 
ings of items in a multivariate mode (Sneath and Sokal 1973). The analytical 
approach used in cluster analysis does not require the data to follow any mathe- 
matical distribution. The items grouped (in this case, stations in Lake Huron) 
are simply compared for overall level of similarity and grouped according to 
some specific clustering algorithm. The choice of the clustering algorithm is 
somewhat arbitrary as are the interpretations of the results. One advantage of 
cluster analysis is that the technique will produce a logical, hierarchical 
grouping to the data even if such a grouping is not initially apparent. 

The complete data set is displayed in the results, providing an evaluation 
of how all the items fit together into one hierarchical scheme. The major 
drawback to cluster analysis is that the technique will cluster any items into 
a hierarchical pattern irrespective of their true relationship. Random numbers 
can be successfully run through a clustering algorithm and produce reasonable 
results. A high level of critical judgment is required to interpret cluster 
analysis results. This disadvantage often requires a second statistical analyt- 
ical procedure to confirm the cluster results. 

The cluster analysis used to determine homogeneous water masses in Lake 
Huron was conducted as follows: Surface samples ( 1 m) were analyzed using five 
variables (temperature, conductivity, dissolved reactive silica, nitrate plus 
nitrite, and chlorophyll). The data were standardized before analysis, which 
requires subtracting the mean and dividing by the standard deviation for each 
of the five variables. Standardization has the effect of placing each variable 
in the same numerical range which makes them more intercomparable. The Euclid- 
ean distance between each variable was determined and the distance matrix was 
then used in the cluster analysis. The average-distance cluster algorithm was 
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used to form clusters (Sneath and Sokal 1973). The results of the cluster 
analysis were checked by computing the cophenetic correlation coefficient, which 
is a measure of the association between the cluster results and the original 
data (Sneath and Sokal 1973). A high cophenetic correlation indicated the data 
conformed to the cluster results very well, while a low cophenetic correlation 
indicates a poor cluster result. In this study, all the cophenetic correlations 
were above 0.80. The results of the cluster analysis were displayed graphically 
by circling large, well-defined clusters of stations on Lake Huron (Figs. 6.1- 
6.6). Each cruise's water masses are labelled independently. Thus "a" for one 
cruise is not necessarily the same water mass as "a" for other cruises. 

Principal Components Analysis — Principal components analysis is a type of 
statistical ordination (Morrison 1978). This method begins with a correlation 
matrix among variables and reduces the number of original variables to one or 
two new variables. These new variables are called principal components and are 
composed from the original variables. The objective of this analysis is to 
produce two or three principal components which explain the correlation among 
these variables. The variable weights indicate how much each variable contrib- 
uted to each principal component. 

Principal components analysis has been extended in ecological applications 
to allow identification of groups of sampling sites (Pielou 1977). This group- 
ing is achieved by the following computations: The principal component weights 
are multiplied by the value of each variable at each station and added at each 
station. The result is one number for each station for each principal compo- 
nent. This number is called a principal component score. Scores for principal 
component 1 (the most important component and the one that explains the most 
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total variation) are plotted against scores for component 2 (second most impor- 
tant component). Stations which are highly similar will group together on this 
plot. A subjective decision is then made to determine if the groupings deter- 
mined by principal components match the clusters determined by cluster analysis. 
This approach was used for each of the six Lake Huron 1980 cruises separately. 

Discriminant Functions Analysis — Discriminant functions analysis is used 
to confirm that identified groups or segments of data are statistically differ- 
ent from one another (Morrison 1978). Unlike the methods described above, 
discriminant functions analysis will not identify or suggest groups of data. 
Rather, the technique uses a multivariate method to compare variance-covariance 
structure among already identified groups. Thus, discriminant functions analy- 
sis is a logical next step beyond cluster analysis or principal components 
analysis; this analysis tests the hypothesis that the variance-covariance 
matrices for each identified group or cluster of stations are all equal to one 
another. If the hypothesis is rejected, the discriminant functions can then be 
applied to the original data and a new grouping determined. This new grouping 
is then compared to the original grouping, and differences are identified. The 
stations within the new grouping created by the discriminant function analysis 
are not identified, but rather, the number of stations which do not concur 
between the new (discriminant function result) and original groupings are 
tallied. A high (80% or better) concurrence between the new and original group- 
ings indicates the classification supplied by the investigator was a fairly 
accurate appraisal of the true structure present in the data. 

An analysis of the success of the original groupings scheme can be deter- 
mined through calculation of the Goodman-Kruskal Tau statistic (Huberty 1975). 
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This statistic represents the reduction in overall within-group error from use 
of the original grouping scheme compared to chance (random) assignment. The Tau 
statistic is expressed as a fraction representing the proportion of reduction. 
A Tau value of 0.00 represents no reduction in the original grouping assignment 
over random assignment. 

Non-Statistical Identification of Water Masses 

Although a variety of statistical techniques are available to aid in the 
identification of homogeneous water masses, previously the most common approach 
has been intuitive geographic segmentation. This technique, which has proven to 
be surprisingly accurate, relies on one or more trained scientists to divide the 
lake(s) into homogeneous geographic regions. Despite the lack of statistical 
methods, this technique can substantially reduce within-group variance when 
conducted by an experienced limnologist. Lake Huron was the subject of several 
such intuitive segmentation schemes, one of which was used to evaluate lake-wide 
water quality from a 1974 database (IJC 1976) (Fig. 6.7). The results presented 
below compare the 1974 segmentation scheme to a 1980 geographic segmentation 
scheme (determined by the intuitive approach) and to the statistical segmenta- 
tion determined by the cluster/principal component analysis. 

Homogeneous Water Masses 

The segmentation scheme developed from the 1974 study of Lake Huron and 
published in 1976 is shown in Figure 6.7. This scheme determined 19 open lake 
segments. Five segments covered the open lake regions as follows: northern 
basin; eastern southern basin; western southern basin-Saginaw Bay mouth; 
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Georgian Bay; and North Channel. The remaining 14 segments are primarily near- 
shore. The advantages of this segmentation scheme are that it provides a high 
degree of resolution in determining spatial differences in water quality, and 
the segments have fixed boundaries based on the most reasonable separation of 
water masses. The disadvantages of this scheme are that the boundaries between 
segments are fixed and hence inflexible among seasons, and 19 segments may 
unnecessarily divide the lake into too many regions. 

In an effort to counter the last criticism, that Lake Huron does not really 
have 19 distinct water masses, a second geographical segmentation scheme was 
devised (Fig. 6.8). This scheme has only 10 segments, with some stations not 
included in any segment (excluded from the analysis). The advantages of this 
scheme over the IJC scheme are a more logical division of Lake Huron based on 
geographic considerations, and fewer segments. 

Both the IJC and second (future references to the second segmentation 
scheme will call this the EPA scheme) segmentation scheme divide the lake into 
segments based on distinct geographic regions: northern basin; southern basin; 
Saginaw Bay mouth; Georgian Bay; etc. This geographic definition of water mass 
boundaries is reasonable because it takes into account sources of water input 
into Lake Huron (see Chapter Four). 

The IJC segmentation scheme (Fig. 6.7) identifies water masses in the 
Straits of Mackinac, western North Channel, western southern basin, and northern 
nearshore Georgian Bay (four segments). These segments are associated with in- 
puts to Lake Huron from: Lake Michigan, Lake Superior, Saginaw Bay, and Canadian 
Shield rivers, respectively. The EPA segmentation scheme identifies segments 
in the same general areas as the IJC scheme, except nearshore Georgian Bay 
stations are not included. The major difference between the two segmentation 
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schemes is that the EPA scheme (Fig, 6.8) has fewer near shore segments than the 
IJC scheme (Fig. 6.7). 

The other approach used to determine homogeneous water masses was that of 
cluster analysis and principal components analysis. This multivariate approach 
has an advantage over the geographic method in that larger homogeneous water 
masses were identified, and hence Lake Huron was divided into a smaller number 
of segments than with either the IJC or EPA schemes. A further advantage is 
that the cluster analysis scheme does not set fixed water mass boundaries , but 
rather allows the location of boundaries to change among cruises. These chang- 
ing boundaries appear to reflect more accurately the true water mass distribu- 
tion than do the fixed boundaries. 

The cluster analysis segmentation scheme for cruise 1 (Fig. 6.1) determined 
Lake Huron should be divided into eight segments. The open lake was divided 
into two major segments, a southern and a northern basin. Most of Georgian Bay 
formed the third major segment during the April cruise. The remaining five 
segments were all in the nearshore zones and included: eastern southern Lake 
Huron; western southern Lake Huron: northwestern Georgian Bay; central North 
Channel and adjacent Lake Huron; and northwestern North Channel. 

The second cruise (May) analysis also yielded a segmentation scheme with 
eight clusters (Fig. 6.2). This analysis produced only one major open lake 
segment which covered the entire lake. Georgian Bay formed a second major 
segment. The remaining six segments were all nearshore and located in the 
following areas: western northern Lake Huron and eastern Georgian Bay; 
western southern Lake Huron; northwestern Georgian Bay; southeastern North 
Channel; northwestern North Channel; and DeTour and False DeTour passages. 
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The third cruise was conducted in early June 1980 and the analysis produced 
nine homogeneous segments (Fig. 6.3). Similar to the May results, one large 
homogeneous water mass was identified in open Lake Huron, and another segment 
covered most of Georgian Bay. The remaining seven segments were composed of 
several disjointed nearshore segments. These nearshore segments were as fol- 
lows: the mouth of Saginaw Bay and eastern and western southern Lake Huron; 
open southern Lake Huron and north of Saginaw Bay mouth; northeastern Georgian 
Bay; northwestern Georgian Bay; northern Straits of Mackinac and central North 
Channel; western North Channel; and northern, eastern, and western northern Lake 
Huron. 

The fourth research cruise was conducted in July 1980 and yielded a seg- 
mentation scheme (Fig. 6.4) that differed considerably from the first three 
cruises. Eight homogeneous water masses were identified, but most of these 
eight were offshore segments rather than nearshore segments. The only nearshore 
segments are: the mouth of Saginaw Bay and adjacent western southern Lake 
Huron; western North Channel; and northern and western Georgian Bay. The re- 
maining five segments included open Lake Huron as follows: southern Lake Huron; 
eastern northern Lake Huron; western northern Lake Huron; North Channel and 
adjacent Lake Huron; and Georgian Bay and nearshore western Lake Huron. 

The analysis of the fifth (September) cruise yielded a segmentation scheme 
with only six homogeneous water masses (Fig. 6.5). The only nearshore segment 
was extreme eastern and western North Channel. Central North Channel and adja- 
cent Lake Huron composed another homogeneous water mass. The remaining four 
segments roughly divided Lake Huron into geographic quadrants. These quadrants 
were: southern and central Lake Huron; eastern Lake Huron; Georgian Bay and 
adjacent Lake Huron; and northwestern Lake Huron. 
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The November cruise segmentation scheme divided Lake Huron into seven 
homogeneous water masses (Fig. 6.6). Similar to the previous cruise, only one 
nearshore segment was identified in northwestern Lake Huron. The remaining six 
open lake segments were distributed as follows: North Channel and adjacent Lake 
Huron; northwestern Lake Huron and eastern North Channel; western Georgian Bay 
and central Lake Huron; eastern Georgian Bay and south central Lake Huron; and 
southern Lake Huron and the mouth of Saginaw Bay (constitutes two complex seg- 
ments) . 

The representation of homogeneous water masses is only approximate because 
several water masses are not easily represented on a map of Lake Huron. Figures 
6.1-6.6 show only the approximate location of the boundaries of homogeneous 
water masses. Table 6.1 gives the station numbers which compose each homogene- 
ous water mass. This table shows that occasionally one Georgian Bay station may 
group with a geographically distant water mass (cruise 1, segment eight, for 
example). But, the vast majority of homogeneous water masses were geographi- 
cally compact. 

Discriminant Functions Analysis 

The results of the cluster analysis and principal components analysis 
yielded substantially different homogeneous water masses than the geographic 
procedure used in the past. Linear Discriminant Function Analysis was used 
to confirm if the segmentation scheme developed by the cluster analysis was sta- 
tistically significant, i.e., the differences among segments were statistically 
significant using five variables. The Discriminant Analysis was applied to the 
IJC, EPA, and cluster analysis segmentation schemes in the same fashion so the 
three schemes could be compared, 
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Table 6.1- Homogeneous water masses determined by cluster analysis/principal 
components analysis. 



Water Mass Stations 

Cruise 1 

1 1, 8, 10, 18-20, 22-23, 25, 34-36, 41, 

44-49, 51-56, 61, 62, 65, 66 

2 2, 6, 9, 12, 15-16, 21, 24, 26, 27, 42, 

43, 90-93 

3 50, 101-124, 126, 128-137, 139-140 

4 13, 14, 17, 63 

5 64, 67, 125, 138, 141-144 

6 3-5, 11 

7 68-72, 127 

8 57-58, 60, 75-77 

Cruise 2 

1 1, 11, 13, 22-23, 25, 34-36, 41, 46-47, 

55-56, 63-64, 105, 108, 110, 115-116, 
128, 139, 141-142 

2 2, 6, 8-9, 12, 15-16, 18-21, 24, 26-33, 

37-39, 43-45, 48-54, 57, 59, 61-62, 
65-66, 90-93 

3 101-104, 106-107, 111-114, 117-122, 

129-136, 140 

4 58, 73-78, 80, 82-83, 85-86 

5 60, 67, 79, 87, 89 

6 7, 10, 14, 40, 94 

7 68-72, 81, 84 

8 88, 109, 123-127, 138, 143, 144 



(Continued) 
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Table 6.1. Continued . 



Water Mass 



Stations 



1 
2 
3 



1 

2 

3 

4 



Cruise 3 

1-2, 7-8, 10-11, 13-14, 16-20, 22, 94 

6, 21, 23-25, 50, 90-91 

101-103, 106-107, 111-114, 117-122, 
129-135 

104, 105, 108, 110, 115-116, 123-124, 
128, 136-137, 139-142 

9, 12, 15, 26-27, 29, 31-33, 37-38, 43, 
45, 48-49, 51-54, 57, 61, 66, 92-93 

28, 30, 35, 39-42, 46-47, 55-56, 58-59, 
62-64 

65, 67, 70, 72-78, 80-86 

68, 69, 71 

87-89, 109, 125-127, 138, 143-144 

Cruise 4 

6-7, 9-12, 15-16, 24, 26, 28, 33, 90-94 

23, 25, 34-35, 41, 46-47, 55-56, 62, 
101-108, 111-112, 114, 117, 119-120, 122, 
128-130, 137, 139 

13-14, 17-18, 20, 22, 113, 118, 121 

17, 19-32, 37, 39-40, 42, 45, 49-54, 59 
59, 133-135 

109-110, 115-116, 123-127, 136, 138, 
140-142 

57, 60, 73-87 

38, 68, 69, 70 
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(Continued) 



Table 6.1. Concluded. 



Water Mass 



Stations 



1 
2 

3 
4 

5 
6 



Cruise 5 

1-21, 24, 26-28, 30, 33, 37, 39, 90-94 

25, 34-35, 46-47, 54-57, 61, 63, 65-66, 
125, 142 

29, 31-32, 38, 40-41, 43 

42, 44-45, 48-49, 51-53, 59, 101-108, 
110-124, 126-141 

58, 60, 67, 70, 72-79, 81-85 

71, 80, 86-89 



1 
2 
3 



5 
6 
7 



Cruise 6 

1-6, 8-11, 18, 20-22, 90-92 

7, 13-14, 19, 23, 25, 34-36, 109 

15-16, 24, 26, 28, 30-31, 33,39, 40-42, 
49, 93, 101-108, 110-124, 131 

27, 29, 32, 37-38, 43-45, 52-53, 125-130, 
132-142 

46-47, 55-56, 62, 65 

48, 50-51, 54, 57, 59-61, 66-68 

58, 70-86 
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The initial hypothesis tested in Discriminant Analysis was that all the 
different segments have the same variance-covariance matrix versus the variance- 
covariance matrices were not equal. For all cruises for each of the three seg- 
mentation schemes this hypothesis was rejected and the alternate hypothesis that 
the variance-covariance matrices were not equal was accepted. The inference 
from these results was that distinctly different homogeneous water masses were 
identified. This analysis was carried an additional step by examining the 
amount of reduction in the variance-covariance matrices for each segmentation 
scheme for each cruise. The results are shown in Table 6.2 as represented by 
the Goodman-Kruskal Tau statistic. The larger the statistic, the larger the 
reduction of the within-group (within-segment) variance and hence statistically 
more distinct homogeneous water masses. Table 6.2 confirms that all three 
segmentation schemes yielded some reduction of the within-group variance over 
random assignment of stations to segments. The cluster analysis/principal 
components analysis approach consistently produced the most reduction in within- 
group variance, while the IJC scheme yielded the least reduction of variance. 



Table 6.2. Goodman-Kruskal Tau statistics for each cruise for the three 

different homogeneous water mass schemes. Number of water masses 
are shown in parentheses behind each Tau statistic. 







EPA 






IJC 


Cluster 


April 


.4794 (18) 


.5318 (10) 


.7994 (8) 


May 


.4470 (18) 


.6576 (10) 


.9091 (8) 


June 


.5925 (18) 


.6301 (10) 


.8527 (9) 


July 


.4092 (18) 


.5846 (10) 


,6234 (7) 


September 


.4852 (18) 


.6738 (10) 


.7515 (6) 


November 


.5450 (18) 


.7025 (10) 


.8951 (7) 
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Homogeneous Water Masses - Ecological Inferences 

The primary motivation behind identifying homogeneous water masses is to 
determine if the water body has spatial structure across several variables, and 
what is the nature of that spatial structure. For a body of water such as Lake 
Huron, identification of spatial structure is not surprising since such a large 
lake cannot be completely mixed both vertically and horizontally. The input of 
several different water types into Lake Huron is a major source of spatial 
variability (IJC 1976). Anthropogenic loadings in the nearshore zone contribute 
further to spatial inhomogenities (Davis e£ ad. 1980). 

The IJC and EPA segmentation schemes (Figs. 6.7 and 6.8) recognized these 
different inputs to Lake Huron and developed a segmentation scheme accordingly. 
But, these schemes lack flexibility. Anthropogenic inputs, nutrient loadings, 
nutrient utilization, etc. all have seasonal components based on physical and 
biological events. These events can lead to large changes in the spatial dis- 
tribution of a suite of variables and thus change the shape, size, and distribu- 
tion of homogeneous water masses. Furthermore, long-term changes in nutrient 
loadings can have a major effect on spatial water mass distribution. Examples 
of physical events which can affect water mass distribution would be formation 
and migration of the thermal bar, thermocline formation, upwelling, and large 
amplitude internal waves. Biological events which readily affect water mass 
composition include intense primary production in algal blooms, vertical migra- 
tion by zooplankton, and anaerobic decomposition in the benthos. Mitigating 
measures (sewage treatment) also have important effects on nutrient loadings and 
distribution. The net result is that the distribution of homogeneous water 
masses in any large lake will not be static but rather dynamic based on a 
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combination of physical and biological events. The geographic approach used in 
the past cannot accommodate the dynamic nature of the lake because water mass 
boundaries are not static. The statistical approach described above does accom- 
modate this dynamic aspect of Lake Huron. 

The segmentation scheme determined by the cluster analysis/principal compo- 
nents analysis yielded results which indicated the influence of thermal bar 
conditions and thermal stratification. The first three cruises were conducted 
during thermal bar conditions. These conditions are characterized by a well- 
mixed, homothermous open lake, and a warmer, thermally stratified nearshore zone 
(Boyce 1974). Because thermal bar conditions occur during periods of high 
spring runoff, nutrient concentrations near shore are often elevated (Davis 
et^ al . 1980). The segmentation scheme developed by the statistical approach 
reflected the development and progression of the thermal bar. 

The April segmentation scheme (Fig. 6.1) shows only a few nearshore seg- 
ments, with the two largest of these in the southern basin. During April 1980, 
thermal bar conditions were just beginning in the southern end of Lake Huron. 

The thermal bars rapidly developed between the first (April) and second 
(May) cruises. This development includes both the movement of the bar offshore 
as well as development of thermal bars in the northern basin (Boyce 1974, Scavia 
and Bennett 1980). The segmentation scheme accurately reflects these physical 
conditions; the segmentation scheme from the May cruise included long nearshore 
segments in both the northern and southern basins (Fig. 6.2). Some of the May 
nearshore segments were much wider than the corresponding April segment, i.e., 
northern Georgian Bay. 
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Thermal bar conditions continued through the third cruise (June). The seg- 
mentation scheme reflects these physical conditions with large, wide nearshore 
segments around all of Lake Huron (Fig. 6.3). In addition, the central southern 
basin is in a separate segment from the remaining open lake. This segment 
accurately reflects the onset of thermal stratification throughout the lower 
southern basin. 

Lake Huron became completely thermally stratified between the third and 
fourth (July) cruises. This change in the physical status of the lake induced 
a large change in the spatial location of homogeneous water masses. Rather than 
a well-mixed, uniform open lake, large water masses were found throughout the 
entire lake (Fig. 6.4). Nearshore segments did not account for most of the 
segments as in the first three cruises. Furthermore, the open lake segments 
were not strictly divided by east-west or north-south boundaries as are often 
used in geographic segmentation, but rather were complex and somewhat meandering 
boundaries . 

The boundaries between the July water masses were less distinct than for 
the three earlier cruises. In July, subtle differences in nutrient concentra- 
tions were the basis of the different water masses. Earlier in the year, sharp 
contrasts in nutrient concentrations across the thermal bar were the basis for 
distinguishing different water masses. As a result, the homogeneous water 
masses identified by the multivariate techniques were a more accurate represen- 
tation of true conditions in Lake Huron for the first three cruises than for 
cruise 4. The Tau statistics in Table 6.2 reflect the change from cruise 3 to 
cruise 4; the statistics for cruises 1 through 3 were .799, .909, and .853, 
respectively, while the statistic for cruise 4 was .623. 
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These results suggest the rather complex distribution of nutrients, chloro- 
phyll, and conductivity in Lake Huron. Physical factors such as currents, depth 
of thermocline, etc. were probably an important mechanism in determining the 
distribution of the homogeneous water masses. This hypothesis is further sup- 
ported by the concept that the plankton, which are passive floaters, are the 
major biological group affecting open lake nutrient concentrations. Therefore, 
the results indicate that physical processes appear strongly to influence the 
location of homogeneous water masses while the composition and abundance of the 
plankton affect nutrient concentrations within homogeneous water masses. 

The fifth cruise was conducted in September with Lake Huron remaining 
thermally stratified. The segmentation scheme (Fig. 6.5) yielded the lowest 
total number of segments of the six cruises (six segments) and the lowest number 
of nearshore segments (one). These results also support the hypothesis that the 
shape and distribution of homogeneous water masses originates primarily from 
physical mixing events. In September, Lake Huron appeared to be well mixed 
in a horizontal plane, but highly variable (stratified) in the vertical plane. 

The sixth cruise was conducted in mid-November when Lake Huron was in the 
transition between thermally stratified and winter homothermous . Under these 
conditions, a deep thermocline exists only in the middle of the lake. Nearshore 
zones are relatively warm and homothermous (see Chapter Three). The segmenta- 
tion scheme for cruise 6 (Fig. 6.6) resembles a mixture of the pre- and post- 
thermal bar schemes; there are some nearshore segments, but also some irregu- 
larly shaped open lake segments. 

Two segments deserve particular attention because of previous research. 
In 1973, a long homogeneous water mass was observed along Lake Huron's western 
shore from the Straits of Mackinac to Thunder Bay (Moll et al. 1976). This 
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water mass was interpreted to represent Lake Michigan water entering Lake Huron 
through the Straits of Mackinac. This same homogeneous water mass was observed 
in three of the six cruises (2, 3, and 6). The other homogeneous water mass 
of historical interest was also a long segment on Lake Huron's western shore. 
In 1974 this segment was observed running from the mouth of Saginaw Bay around 
the "Thumb" of Michigan and into the southern basin. This water mass was inter- 
preted to represent Saginaw Bay water flowing into Lake Huron. This segment was 
only observed during the first three cruises of 1980, and was only partially 
present during cruise 1. 

Homogeneous Water Masses - Statistical Inferences 

The segmentation scheme developed by the cluster analysis/principal compo- 
nents analysis consistently yielded fewer segments with a larger reduction in 
variance than the geographic segmentation scheme (Table 6.3). This result can 
be used to advantage in analysis of data and development of experimental designs 
for further monitoring. The analysis of data is improved with fewer segments by 
way of decreasing the within-group variance. Since the smallest number of 
homogeneous water masses were identified with the cluster analysis/principal 
components analysis approach, each water mass, on the average, will have a 
larger number of stations. This larger number of stations will reduce the 
variance. As a result of the smaller within-segment variances, the statistical 
power in determining differences among water masses is increased (Lindgren 1976, 
Sokal and Rohlf 1969). Table 6.3 shows that the confidence intervals were 
usually smallest for all cruises using the cluster analysis/principal components 
analysis approach versus the geographic approach. The reduction in confidence 
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Table 6.3. Number of samples (N) , means, and 95% confidence intervals for 
soluble silica, nitrate plus nitrite, and chlorophyll for each 
of the three homogeneous water mass schemes. 



Mean 



Confidence Interval 



Northern Basin — July 



IJC 

Silica 

Nitrate 

Chlorophyll 

EPA 

Silica 

Nitrate 

Chlorophyll 

Cluster 
Silica 
Nitrate 
Chlorophyll 



131 

130 

78 

84 
83 
42 

103 

102 

56 



1.301 

.268 

1.305 

1.428 

.277 

1.505 

1.226 

.266 

1.200 



.066 
.004 
.105 

.080 
.005 
.119 

.068 
.004 
.083 



Southern Basin — April 



IJC 

Silica 

Nitrate 

Chlorophyll 

EPA 

Silica 

Nitrate 

Chlorophyll 

Cluster 
Silica 
Nitrate 
Chlorophyll 



36 


1.481 


36 


.332 


19 


2.347 


28 


1.510 


28 


.283 


8 


1.738 


62 


1.491 


62 


.281 


19 


1.658 


Georgian Bay - 


- November 


44 


1.247 


44 


.277 


12 


1.108 


44 


1.247 


44 


.277 


12 


1.108 


146 


1.248 


146 


.273 


50 


1.232 



.014 
.027 
.297 

.018 
.005 
.161 

.012 
.004 
.100 



IJC 

Silica 

Nitrate 

Chlorophyll 

EPA 

Silica 

Nitrate 

Chlorophyll 

Cluster 
Silica 
Nitrate 
Chlorophyll 



.031 
.004 
.096 

.031 
.004 
.096 

.038 
.003 
.052 
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intervals is especially important in the nearshore segments. In these segments, 
the geographic approach yields too many segments with very large variances 
because of the small sample sizes. The cluster analysis/principal components 
analysis tends to remove the problem by having fewer segments with larger sample 
sizes. 

The inference for improved experimental design originates from more effi- 
cient sample station location. Parsimonious placement of sampling stations for 
a lake-wide monitoring program would dictate at least three stations per homoge- 
neous water mass. The IJC segmentation scheme would therefore require at least 
54 stations, while the cluster analysis scheme would require from 18 to 27 sta- 
tions. In practice, it is unlikely a sampling grid could be tailored precisely 
to any segmentation scheme; the identification of segments by the cluster/prin- 
cipal component analysis approach cannot be achieved until after the samples 
have been collected. But, if a sampling grid of approximate uniform station 
density were developed (such as for 1980), the segmentation scheme may aid in 
judicious choice of sample analysis. For example, phytoplankton cell counts are 
normally counted long after the samples are collected. Given a desire to count 
three samples from each segment, fewer segments means less counting. 

Finally, identification of homogeneous water masses is an important step 
for long-term trend analysis. Only by comparing similar water types can trends 
be examined across years. Because the cluster/principal components analysis 
approach shows that homogeneous water masses change considerably among seasons, 
long-term trends should be inferred only from similar water types in similar 
seasons. Without such an approach, the statistical power in resolving long-term 
trends will be very low. These statistical concepts are used in Chapter Eight 
to develop long-term trends in Lake Huron water quality. 
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CHAPTER SEVEN 

STATISTICAL INFERENCES: LAKE HURON, 1980 

By Russell Moll 

The surveillance and study of Lake Huron in 1980 yielded a large, extensive 
database. The geographical area studied during the six cruises from April to 
November produced a large range in many of the measured variables (Table 7.1). 
Complete ecological interpretation of the Lake Huron data was not possible by 
non-statistical methods alone, i.e., data plots and contour maps. For example. 
Chapter Four, the descriptive physical-chemical limnology of Lake Huron, was 
restricted to comparisons among surface samples and vertical profiles of only 
one station for each of the major basins. The purpose of this chapter is to 
invoke several statistical analytical methods to expand the inferences about 
Lake Huron from the 1980 database. The techniques presented below have been 
used before in ecological studies (Langbein and Lichtman 1978). 

The results presented in Chapter Four had two shortcomings which were 
somewhat removed by the statistical approach used below. These shortcomings 
were: the data were analyzed only at one depth or for one station at a time, 
and comparisons among cruises were limited. The multivariate techniques 
presented below were used to take large, complex databases and reduce the data 
into a smaller database or into one meaningful statistic. By reducing the 
database, ecological inferences were more apparent. 

Statistical Methods 

Two groups of statistical techniques were used to yield broader inferences 
from the 1980 Lake Huron data: correlation analyses and tests of multivariate 
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Table 7.1. Complete list of chemical and physical variables measured during 
1980 Lake Huron study. 



Alkalinity 

Ammonia 

Chlorophyll 

Dissolved Oxygen 

Kjeldahl Nitrogen 

Particulate Organic Carbon 

Particulate Organic Nitrogen 

pH 

Secchi Depth 

Soluble Nitrate-Nitrogen 

Soluble Reactive Phosphorus 

Soluble Reactive Silica 

Specific Conductance 

Temperature 

Total Filterable Phosphorus 

Total Phosphorus 

Transmissivity 

Turbidity 



hypotheses. The computations were performed using The University of Michigan 
statistical software package MIDAS (Fox and Quire 1976). Correlation, unless 
otherwise specified, refers to the simple Pearson Product-Moment zero-order 
correlations ( Sokal and Rohlf 1969). The data set is not totally balanced and 
has some missing observations. As a result of these empty or missing cells, all 
zero-order correlations were computed with a missing-data algorithm. In effect, 
all pairs of observations were used to calculate each correlation coefficient, 
not just pairs from complete cases (Fox and Guire 1976). Computation of higher- 
order correlations and factor analysis required complete cases. 
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Factor analyses were conducted for each cruise using the complete-case 
correlation matrices. The number of factors extracted per cruise was based on 
the Scree Test (Kim and Mueller 1978). Five factors were extracted for cruises 
1, 2, 5, and 6, four factors for cruise 3, and three factors for cruise 4. 
The values of the coramunalities were estimated using the iterative principal 
axis factor solution as described in Harman (1967). Convergence of coramunali- 
ties was deemed sufficient when succeeding estimates differed by less than 
0.001, or after twenty iterations. An orthogonal varimax rotation was performed 
on the factor matrix (Kim and Mueller 1978). 

Results 
Correlation Analyses 

The correlation analyses were conducted in a three-step procedure. 
First, zero-order correlations were calculated among all pairs of variables. 
The results of calculations were used to reduce the number of variables analyzed 
by removing highly correlated, and thus somewhat redundant, variables. Second, 
a reduced data set was composed of 10 variables (Table 7.2) which represented a 
non-redundant data set. A factor analysis was conducted on the reduced correla- 
tion matrix for each cruise. Third, results from each cruise were compared to 
one another to determine seasonal changes in relationships among variables. 

The correlation matrices show a high, positive correlation for all six 
cruises among chloride (CI), conductivity (Cond), and alkalinity (Alkal). 
These correlations were all highly significant at the P <.01 level. The values 
of these correlations ranged from 0.821 to 0.966. A further indication of the 
somewhat redundant nature of these three variables was the similar correlations 
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Table 7.2. Reduced set of variables used in multivariate analysis. 

Temperature 

Specific Conductance 

pH 

Soluble Nitrate-Nitrogen 

Soluble Reactive Phosphorus 

Soluble Reactive Silica 

Ammonia 

Total Phosphorus 

Total Filterable Phosphorus 

Kjeldahl Nitrogen 

among chloride, conductivity, and alkalinity to other variables. For example, 
the CI to total phosphorus correlations ranged from -0.081 to 0.155 across 
cruises; the conductivity to total P correlations ranged from -0.148 to 0.043; 
the alkalinity to total P correlations ranged from -0.186 to 0.025. Likewise, 
the correlations among ammonia (NH3) and chloride, conductivity, and alkalinity 
were similar in magnitude and seasonal response. Table 7.3 gives the correla- 
tions among chloride, conductivity, alkalinity, total P, and NH3. These results 
indicated that the three variables provided approximately the same ecological 
inference. 

A similar type of analysis determined a high level of redundancy between 
(or among) the following sets of variables: chloride to alkalinity to conduc- 
tivity to sulfate; particulate organic carbon to particulate organic nitrogen; 
unfiltered ammonia to filtered ammonia; unfiltered soluble nitrate plus nitrate 
to filtered soluble nitrate plus nitrite. The highly correlated pairs or sets 
of variables listed above do not become decoupled across seasons. These results 
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Table 7,3. Cruise 1 missing data correlations for 1980. 

























VARIABLE 


MEAN 




STD DEV 




N 


CORR 


T-STAT 


SIGNIF 


10 


. ALKALINI 


72.898 




8.3331 




435 


.9749 


91 .025 


0. 


1 1 


.CONDUCT I 


193.67 




21 .352 












10 


ALKALINI 


72. 137 




8.6551 




378 


-.6373 


-16.034 


.oooo 


17 


. NH3 


.32354 


-2 


.50183 


-2 










10 


ALKALINI 


72.718 




8.38 12 




403 


-.0445 


-.89263 


.3726 


26 


TOTAL P 


.57385 


-2 


.41385 


-2 










1 1 


.CONDUCT I 


191 .89 




22.262 




378 


-.6321 


-15.816 


.0000 


17 


. NH3 


.32354 


-2 


.50183 


-2 










1 1 


.CONDUCT I 


193.33 




2 1 . 590 




403 


-.0226 


-.45365 


.6503 


26 


.TOTAL P 


.57385 


-2 


.41385 


-2 










17 


. NH3 


.33526 


-2 


.52251 


-2 


346 


.2989 


5.8100 


.0000 


26 


.TOTAL P 


.56289 


-2 


.36328 


-2 











Table 7.3. Cruise 2 missing data correlations for 1980. 







MEAN 



















VARIABLE 


STD DEV 




N 


CORR 


r-STAT 


SIGNIF 


10 


.ALKALINI 


72.453 




8.5115 




502 


.9600 


76.615 


0. 


1 1 


CONDUCT I 


190.57 




22. 166 












10 


.ALKALINI 


72.802 




8.5160 




456 


.8416 


33 . 205 


0. 


21 


.CHLORIDE 


5.0923 




.84531 












10 


.ALKALINI 


72.374 




8.5828 




491 


-.4698 


-1 1 .770 


.0000 


16 


. NH3 


.25418 


-2 


.34262 


-2 










10 


.ALKALINI 


72.415 




8.5032 




500 


-. 1856 


-4.2152 


.0000 


23 


TOTAL P 


.51774 


-2 


. 14927 


-2 










1 1 


.CONDUCT I 


191 .39 




22.554 




457 


.8885 


4 1 . 302 


0. 


21 


.CHLORIDE 


5.0919 




.84443 












1 1 


CONDUCT I 


190.36 




22.329 




492 


- . 4030 


-9.7475 


.0000 


16 


.NH3 


.25407 


-2 


.34228 


-2 










1 1 


.CONDUCT I 


190.46 




22. 110 




501 


-. 1483 


-3.3509 


.0009 


23 


.TOTAL P 


.51756 


-2 


. 14917 


-2 










21 


.CHLORIDE 


5.0837 




.85147 




447 


-.2492 


-5.4282 


.0000 


16 


.NH3 


.25414 


-2 


.35520 


-2 










21 


.CHLORIDE 


5.0868 




.84089 




455 


.0049 


. 10337 


.9177 


23 


.TOTAL P 


.52121 


-2 


. 151 16 


-2 










16 


. NH3 


.25429 


-2 


.34296 


-2 


490 


.2410 


5.4849 


.0000 


23 


.TOTAL P 


.5174 1 


-2 


. 15015 


-2 
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Table 7.3. Cruise 3 missing data correlations for 1980. 





VARIABLE 


MEAN 





STD DEV 


r 


N 


CORR 


T-STAT 


SIGNIF 


10 


ALKALINI 


73.785 




8.3941 




543 


.9192 


54.288 


0. 


1 1 


.CONDUCT I 


194 .46 




21.671 












10 


ALKALINI 


73 . 809 




8.3947 




544 


.8209 


33.461 


0. 


19 


.CHLORIDE 


5. 1206 




.76946 












10 


.ALKALINI 


73.694 




8.3280 




540 


-.3010 


-7.3210 


.0000 


14 


.NH3 


.27759 


-2 


.28251 


-2 










10 


.ALKALINI 


73.807 




8.3871 




545 


-.0697 


-1 .6288 


. 1039 


21 


.TOTAL P 


.46552 


-2 


. 11504 


-2 










t 1 


CONDUCT I 


194.47 




21.691 




542 


.9334 


60.437 


0. 


19 


CHLORIDE 


5. 1192 




. 77053 












11 


.CONDUCT I 


194. 19 




2 1.550 




538 


-. 1688 


-3.9658 


.0001 


14 


.NH3 


.27732 


-2 


.28299 


-2 










1 1 


CONDUCT I 


194.46 




21.671 




543 


.0072 


. 16818 


.8665 


21 


.TOTAL P 


.46552 


-2 


. 11525 


-2 










19 


CHLORIDE 


5. 1124 




.76732 




539 


-.0496 


-1 . 151 1 


.2502 


14 


. NH3 


.27774 


-2 


.28275 


-2 










19 


CHLORIDE 


5. 1206 




.76946 




544 


.0602 


1 .4048 


. 1607 


21 


TOTAL P 


.46575 


-2 


. 11502 


-2 










14 


NH3 


.27759 


-2 


.28251 


-2 


540 


.2140 


5.0822 


.OOOO 


2 1 


TOTAL P 


.46567 


-2 


. 11545 


-2 











Table 7.3. Cruise 4 missing data correlations for 1980. 

























VARIABLE 


MEAN 




STD DEV 


f 


N 


CORR 


T-STAT 


SIGNIF 


10 
1 1 


.ALKALINI 
.CONDUCT I 


73.459 
190.69 




8.2832 
19. 144 




673 


.9656 


96. 144 


0. 


10 
19 


.ALKALINI 
CHLORIDE 


74.042 
5. 1 174 




7.9599 
.75999 




478 


.8782 


40.057 


0. 


10 
14 


.ALKALINI 
. NH3 


73.459 
.48945 


-2 


8.2827 
. 36309 


-2 


673 


.0556 


1 .4428 


. 1495 


10 
21 


.ALKALINI 
TOTAL P 


73.459 
.49565 


-2 


8.2827 
. 16008 


-2 


673 


.0245 


.63587 


.5251 


1 1 
19 


CONDUCT I 
.CHLORIDE 


191 .83 
5. 1171 




18.953 
.75921 




479 


.9130 


48.883 


0, 


1 1 
14 


CONDUCT I 
.NH3 


190.67 
.49004 


-2 


19. 138 
.36278 


-2 


673 


.0313 


.8 1165 


.4173 


1 1 
21 


CONDUCT I 
TOTAL P 


190.67 

.49572 


-2 


19. 138 
. 16013 


-2 


673 


-.0427 


-1 . 1066 


.2689 


19 
14 


CHLORIDE 
NH3 


5.1171 
.51754 


-2 


.7592 1 
.37482 


-2 


479 


.0637 


1 .3943 


. 1639 


19 

21 . 


CHLORIDE 
TOTAL P 


5. 1171 
. 50305 


-2 


.7592 1 
. 17356 


-2 


479 


-.08 10 


-1 . 7740 


0767 


14 , 
21 


NH3 
TOTAL P 


.48947 
.49583 


-2 
-2 


.36283 
. 16003 


-2 

-2 


674 


.2989 


8. 1189 


.OOOO 
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Table 7,3. Cruise 5 missing data correlations for 1980. 

























VARIABLE 


MEAN 




STD DEV 


N 


CORR 


T-STAT 


SIGNIF 


10 


. ALKALINI 


74.4 18 




6.2588 




593 


.9175 


56.094 


0. 


1 1 


.CONDUCT I 


187 . 1 1 




14.397 












10 


.ALKALINI 


74.353 




6.3230 




597 


.8501 


39.368 


0. 


19 


.CHLORIDE 


5. 1802 




.54706 












10 


.ALKALINI 


74.356 




6.3181 




598 


.0929 


2.2771 


.0231 


14 


.NH3 


. 16689 


-2 


. 19085 


-2 










10 


ALKALINI 


74.356 




6.3181 




598 - 


.0527 


-1 .2877 


. 1984 


21 


TOTAL P 


.47669 


-2 


. 12344 


-2 










1 1 


CONDUCT I 


187.08 




14.403 




594 


.8789 


44.831 


0. 


19 


CHLORIDE 


5. 1850 




.53752 












1 1 


CONDUCT I 


187.08 




14.391 




595 


. 1173 


2.8754 


.0042 


14 


.NH3 


. 16622 


-2 


. 19003 


-2 










1 1 , 


CONDUCT I 


187.08 




14.391 




595 


.0431 


1 . 0500 


.2941 


21 , 


TOTAL P 


.47555 


-2 


. 12301 


-2 










19, 


CHLORIDE 


5. 1793 




.54641 




599 


.092 1 


2.2595 


.0242 


14, 


. NH3 


. 16678 


-2 


. 1907 1 


-2 










19 


CHLORIDE 


5. 1793 




.5464 1 




599 - 


.0800 


- 1 . 96 1 7 


.0503 


21 , 


TOTAL P 


.47661 


-2 


. 12333 


-2 










14 , 


.NH3 


. 16667 


-2 


. 19057 


-2 


600 


. 1770 


4.3967 


.0000 


21 


TOTAL P 


.47638 


-2 


. 12335 


-2 











Table 7.3* Cruise 6 missing data correlations for 1980. 



VARIABLE 

10. ALKALINI 
1 1 CONDUCT I 

10. ALKALINI 
19 CHLORIDE 

10. ALKALINI 
14 .NH3 

10. ALKALINI 
21 .TOTAL P 

11 .CONDUCT I 
19. CHLORIDE 

1 1 .CONDUCT I 
14 .NH3 

1 1 CONDUCT I 

2 1 -TOTAL P 

19. CHLORIDE 
14 .NH3 

19. CHLORIDE 
2 1 .TOTAL P 

1 4 . NH3 

21 . TOTAL P 



MEAN 

75. 192 
201 .04 

75.265 
5.2712 

75.274 
.25696 

75.274 
.52241 

201 .02 
5.2634 

201 .04 
.25359 

201 .04 
.52287 

5.2712 
.25708 

5.2712 
.52256 

.25696 
.52241 



STD DEV 

7. 1858 
17.402 

7 . 1 704 
.65656 

7. 1661 
-2 . 16720 -2 

7. 1661 
-2 . 14682 -2 

17.413 
.65784 

17.402 
-2 . 16552 -2 

17.402 
-2 . 14712 -2 

.65656 
-2 . 16736 -2 

.65656 
-2 . 14693 -2 

-2 . 16720 -2 
-2 . 14682 -2 



N 
474 

480 

48 1 
481 

473 
474 
474 
480 
480 
481 



CORR 
.9479 



T-STAT 
64.674 



.9176 50.475 



1231 2.7141 



.0325 



.71204 



9593 


73.699 


0510 


1 . 1 101 


0135 


-.29421 


1081 


2.3766 


0244 


.53449 


3496 


8. 1655 



SIGNIF 
0. 



.0069 



.4768 



.2675 



.7687 



.0179 



.5932 



.0000 
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indicate that a more efficient sampling program could be designed in the future 
based on eliminating redundant measurements. 

Correlations among variables cannot be used to infer cause-and-ef f ect 
relationships (Sokal and Rohlf 1969). But, the results of correlations analysis 
do indicate which variables covary together and if that covariance is positive 
or negative. 

The correlations from cruise 1 reflect the well-mixed conditions typical 
of late winter conditions (Table 7.4). There were few large correlation coeffi- 
cients, and none of these large correlations involved depth or temperature. 
Inverse thermal stratification was present and some vertical stratification 
of nutrients accompanied these conditions (see Chapter Four). But, there was 
little areal change in nutrient concentrations and chlorophyll. Table 7.5 shows 
relatively small ranges for all variables indicating uniform spatial distribu- 
tion of all nutrients. The inference from these results is that spring algal 
growth was insufficient by cruise 1 to produce spatial patterns in nutrient 
levels below normal late winter maximums. Large negative correlations among 
Secchi depth and ammonia, total P, and chlorophyll were observed, as were nega- 
tive correlations among pH, conductivity, soluble silica, and ammonia. These 
negative correlations suggest that some low nutrient values were found in 
regions of high conductivity water, e.g., the mouth of Saginaw Bay. Although 
outer Saginaw Bay was homothermous during cruise 1, the inner bay supported 
sufficient algal growth to influence nutrient concentrations of water flowing 
from the bay into open Lake Huron (Smith et al. 1983). 

The second cruise (9-21 May) was conducted during thermal bar conditions. 
Spring phytoplankton blooms, common during thermal bar conditions (Davis et al. 
1980) were evident from elevated chlorophyll levels between the thermal bar and 
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shore. Similar to results from 1974 (Davis ejt al . 1980), high nitrate-nitrogen 
concentrations were observed inshore of the thermal bar along with high chloro- 
phyll concentrations (Table 7.6). Large nitrate-chlorophyll (+.524) and temper- 
ature-chlorophyll correlations (+.450) indicated this association. The cruise 2 
correlation matrix (Table 7.7) yielded several other high correlations among 
temperature and dissolved nutrients. These correlations suggest that thermal 
stratification in the form of a thermal bar and its associated algal bloom had 
a major influence on nearshore nutrient concentrations. 

Thermal bar conditions continued through the third cruise (27 May to 
7 June) although the thermal bar was much further offshore than in cruise 2. 
The more advanced stages of the thermal bar in early June allowed some degree 
of vertical thermal stratification to develop inshore of the thermal bar. As a 
result, three distinct limnological conditions were present simultaneously in 
Lake Huron during cruise 3: homothermous (in the center of the lake), thermal 
bar (well offshore), and thermally stratified (nearshore). This mixture of 
limnological conditions yielded large differences in nutrient concentrations 
between the nearshore zone and offshore zone. Table 7.8 shows a larger range 
for most of the variables during cruise three than was found in cruises 1 or 2. 
Analysis of the entire cruise 3 data base produced relatively few high correla- 
tions (Table 7.9). Large correlations were found among temperature, depth, 
conductivity, POC, and Secchi depth. The high nitrate-chlorophyll association 
inside the thermal bar observed during cruise 2 was not present during cruise 3. 

By cruise 4 (18-30 July) Lake Huron was thermally stratified in both the 
offshore and nearshore zones. The formation of a thermocline restricted ver- 
tical mixing and resulted in surface nutrient depletion from algal growth 
(Wetzel 1975). As a consequence, the major source of variation in nutrient 



222 



! - 
I E 

i i 



r* ^ n in r^ 
00 00 (0 ^ ex 
U) CM r^ CN "^ 



tn 

00 



00 



^ O 
^ n 
00 n 



00 

in 



n 






f^ 
^ 



in -^ 
0) o 

OJ 00 



i/) 
tn 



Z 
3 



I/) 



o 

10 

m 



^ ic 

O o) 

O CN 



o 



10 



in 

00 

n 



in 

n 

CN 



00 00 
CN CN 

TT O 



^ — to 



n 
n 
n 



n 



CN O 

n CN 
rr CM 



rj 00 

00 CM 



O 

00 

ON 



o 



CO 

a 



r* m 

— 00 

n CM 

00 



>- rj -^ 



in 'ir ID 

00 ^ n 

in «p- '- 

n ID 

•^ CM (0 

CM 



00 
CM 
CM 



n 
in 

00 



I 

u> in 

rr en 

CM *- 

1^ CM 

CM n 



! uj 



O t- 
00 -^ 
in 0) 

CM 
00 
CM ^ 



to 



to 00 

in in 

So" 
O) n 

00 "^ ^ 



O 
in 

CM 



O) 

CM 



'^ O 

0) n 
in 



O) 0) 



t 

n to 
(D in 

*- m 



00 n 
^ O 
(0 



n 

00 



en 

CM 



CM 

in 

00 

0) 



— 00 

in ^ 

CM Ui 

CM 



^ CM 

O CM 

*- CN 

C\ 

cn 



CO 
CO 



> 



X 

< 

z 



8 8 

CM 
00 
10 '- 



O 

o 

CN 



8 8 

in 

00 

n in 

CO CM -^ 



888888888 



in 
n 



in 
n 



O 
0) 
n 



CM *- 

10 10 



O 
(0 



o 

00 



(0 

CM 



V4 

o 

CO 



i 

•■4 

z 



« ^ CM 







8 

CN 



O 

o 

O 
in 



CO 



n 

O 
in 



cn 
O 
in 



CM n n 
O O C0 
in in -^ 



CN 
Oi 



tn 
O 
in 



cn 
O 
in 



CM 

O 
in 



8 o 

in in 



rr to 
O) in 

^ CM 



in 

CSi 



in 



u 



CO 



< 



•- K a 



< 
> 



UJ 

O 



•-« > 

u o 

D 

Q 00 

Z in n 

X o - z 

0. U Q Z 



O CN O X 

a o c a CL c 

^ a o ^ 

O I/) f- -i -J Q< X 

X ^ < ^ C (J 

nt-_jZH-u-«JCjcj 

O Ql O O X C UJ 

ZOcoi^K»-cja.oo 



Ui -^ 



in 



CD 



05 



in 



ID 



00 



223 



3 



Wl 


n 


(0 


O 


in 


f^ 


n 


0) 


CM 


'T 


r* 


n 


TT 


o 


0) 


^ 


»i« 


O 


O 


KD 


00 


n 


f^ 


•v 


f^ 


^• 


r^ 


CN 


0) 


(0 


n 


in 


\n 


in 


r» 


o 


CM 


n 


"^ 


r* 


O 


ID 


O 


n 


o 


CM 


^ 


n 


o 
































H- 


n 


1 


1 


CM 


1 


^ 


0) 


(D 




r- 


•• 


00 


CM 


VD 




cr 












in 


^ 








(0 










D 
































^ 
































un 


C 


^ 


n 


r? 


in 


n 


r* 


in 


in 


O 


(0 


o 


n 


^ 


rj 


in 


n 


CM 


n 


n 


^ 


CM 


in 


CM 


*• 


^ 


ID 


n 


in 


0) 


O 


Ui 


f^ 


(13 


in 


*» 


IT 


rr 


*- 


in 


u) 


0) 


^ 


CM 


00 


h' 


in 


































3 


*- 






^» 


1 


in 


rr 


^ 




*- 


in 


^- 




<^ 




UJ 








1 
























^ 
































wo 

































o 

00 
ON 



o 



0) 



o o 



00 '- 

0) 0) 

CM (0 

CM 



n 


w 


•^ 


p^ 


CD 


n 


TT 


00 


05 


05 


r^ 


00 


in 


in 


00 


00 


O 


o 


^• 


in 


^ 


CM 


n 


«N 


00 


o 


CD 


O 


in 


00 


r* 


CM 


n 


^ 


00 


CO 


n 


^ 


n 


»- 


^ 


00 


in 


O 


r* 


CM 


n 


CM 


P3 


n 


^ 


^ 


tn 


U) 





CO 

C 

o 



ID 
00 

(D 

CM 



00 0) 

ID O) 

CN in 

n CM 



^ 



CM 

in 



(D 

o 

01 



BO -p- — 



0) 

in 



in 

CM 

O) 



(D 
00 

CM 



o 

CM 



ID 

CM 



CM 

in 
in 

(D 



O) n 

rr CM 

05 ^ 

O) r- 

CM 



in 

CM 
00 



u 

o 
a 



CO 
T3 

bO 

c 

•H 
CO 
CO 

•H 

H 

CVJ 

<w 

CO 



z 


8 


8 


8 


8 


8 


g 


8 


g 


8 


8 


8 


Z 




in 


(D 






o 


in 


o 


00 


CM 


05 




n 




r* 


»» 




00 


(D 


^ 


(D 


CM 


05 


X 


r** 


n 




in 


TT 


n 


r- 


CM 




n 


^ 


< 


♦» 


^ 


00 


CM 


^ 








CM 






Z 












CM 

1 




D 




1 


CM 


1 


§ 


§ 


8 

05 


8 


8 

00 


1 


8 


O 


§ 


§ 


§ 


kii« 


o 


r^ 


05 






(D 


o 


*- 


o 


00 


z 








00 


O 


^- 


^- 


CM 


tn 


(D 


CM 




^» 


CM 


f^ 


05 
















z 


























f^ 


in 


^ 


in 


n 


O 


in 


r^ 


r^ 


in 


in 


z 


rr 


rr 


^ 


^ 


n 


n 


T 


n 


'^ 


TT 


^ 




in 


in 


in 


in 


in 


in 


tn 


in 


in 


in 


in 



o 


Q 


Q 


g 


g 


g 


h- 


(D 


CM 



O 


8 


g 


Q 


rj 


o 


Q 


^ 


in 


0) 


*- 





^ ^ n 



1—1 

CO 
H 



< 



QC 

< z 

^ »- a 



a. 

< 
> 






D 

Z 

o 
u 





O 


CM 


CJ 






I 








a 


o 


c 


a 


a 


O 




,^ 




o 


I/) 


h- 


~i 


-J 


o: 




I 








kM 


< 


*m 


o 




u 


r) 


H- 


^ 


2 


p» 


U. 


^ 


u 


u 




cr 


o 




o 




I 


o 


UJ 


2 


o 


m 


iii 


>— 


h- 


(J 


c 


tn 



224 



a> — 


•• mm 


« 


^ 


XD - 




fN .^ 




f^ m 


Q ID 
O « 


eo o» 


■ w 


"-^ « 


CM in 


'» 


• *i» 


• -mr 




IN >- 



o 

00 



u 

o 



CO 

o 



CO 
CO 

a; 

> 



XJ ' 



O 
CO 



CO 



3 



CX5 
H 



c o 



in (\ 

O O 
— in 





?l 


(0 «M 


CO in 


(0 in 


o — 


< 






(0 0) 


(^ m 


tf in 


n n, 








n ^ 


n CN 


CM CM 


(0 in 


n o 




- 










CM ♦- 

S 

CK 




n — 


w-* 


f^ — 


n -* 


(T 






0) 0) 


0) — 


CN tf 


OO (C 


O) -» 






1^ o» 


in 01 


OS in 


- in 


r- ^ 


' z 




— ^ 


•• n 


— CN 


— CM 


- in 


O 


— 












CM 5i 

CM 

o 


o -» 


W — ' 


n — 


CM — 




O 


10 


- 0) 


(0 Q 
(0 O 


f* P5 


V XS 


0» <£ 


in^ 




2 « 


0) 01 


- in 


o in 


m T 


•^ 


O "« 




— n 


CM CM 


n CM 


r* m 


0) a 
— «/i 

o 

0. 


0-' 


«• «» 


o- 


nr — ' 


o — 


r*' 


9 


28 


« — 


m «T 


in vo 


C 10 


n •-' 




(?) o 


CM 0» 


T in 


V in 


V r- 




OiP 


o in 


CM ^ 


O CM 


O CM 


- in 


a ce 


■ "-^ 


• "-» 


• V 


• '.^ 


■ "— 


• w 


- o 



r» •-» 


gr 


o- 


T — 


o> — 


(0 




S8 


o ^ 


«y vo 


in vc 


ij) -. 




o O 


- 0) 


CM in 


n in 


in V 


■ o 


o m 


n m 


o -» 


m CM 


CM CM 


<N in 


r* O 


• -^ 


• «» 




• ■— 


• ■— ' 




- Z 



Q 


r) — 


0» — 


w «» 


V «» 


o *-» 


CM -» 


n -• 


CD •-« 


^ 


Q 


r> VH 


in CM 


CM CN 


n OJ 


-o 


(0 n 


0) CM 


m CM 


CN -» 




- OJ 


V o» 


CM 0) 


CM a 


^ 0) 


r* a 


- in 


tr in 


10 O 


j 


— V 


CM n 


m ^ 


— n 


CN n 


•«- V 


°2 


®2 


c^ in 


f» --^ 




• -" 


eo — . 


0-* 


ot^ 


m-^ 




V 


n 


10 0» 


i0 r» 


CM n 


n CM 


CM CM 


(0 CM 


«0 — 


8 O 




n 


CM in 


n (0 


r* iO 


8s 


r^ (0 


- (0 


m (0 


" O 


V — ' 


•" ^ 


CM — 


O - 


O- 


CM — 




n vo 


- (0 


in « 



iS2 



eo -* 
r* CM 
^ O 
CM in 



I 


0) -» 

CM •• 


OO) 


o» o 
— m 


•T ■— • 

(N n 

r- O 

Oin 


n CM 


if 




m -' 

— 0) 
CM ^ 


^ 10 

in in 

CM CN 


n "^ 
n vo 
n in 
en CN 


O 

O T 


5 

i 

- o 


CM —■ 


0»-' 
V0 CM 

O (0 




CO -• 

in CM 
0) o 
•^ m 


00 CM 


m — « 

^ O 
in in 


CM 0» 


m Q 

CN O 


«0-* 

O ^ 


CO -» 

in (0 
CM m 

V (N 


O 10 
vo in 

V CM 


en ^ 
CM m 


.£ 


CN O 

o o 
CM in 


0) n 
n vo 
m - 


^ 0) 

— * 


V0 — 
O r) 


5S 


"- CN 


CN Q 
O O 


li 


li 


0) — 
0) (0 

V m 

^ CN 


00 — . 

T vo 

a in 

n CN 


S- 

m in 


< 
tt 

1 

■ UJ 



n — 


tn -* 


vo -» 


^ ^ 


in -' 


o -. 


fM — « 


m -' 


~ ^ 


iO -^ 


CM ■— ■ 


r» "^ 


in -» 


5 Zl 


P» (N 


m n 


vo en 


10 CM 


Q cn 


in n 


OJ CM 


CM Q 
CM O 


n — 


0) "^r 


in vo 


n vo 


c o 


o o 


- C 


00 10 


^ 0) 


2 o 


O in 


8£ 


vo O 


W 0) 


i£ in 


00 in 








CM — 


O -v 


O in 


— in 


- in 


— ^ 


OCM 


O CM 



to c^ u) m IT z 



f ^ 



225 



CM — n — o 

in — — — r* P5 • k 
CM rj on ^ ^» CN 



u) —» in ** o ** «D 
<fltn «>» v'- to-^ 
mv in^ ^^ ipn 



^in ••« O)** o>** O"^ 
oDv n^ rt^ ^•" wn 
om ••« on OCT — p* 



z 

eo 



«>«» n-* ««-' r*-^ r-—' 0) 

— in •-tfJ — f* V— 01— — — ' 

(pv «n o»n »•• in— <i)n 

cNin — m Om — ct c*ct mr* 



V p* 0) in 



j-» 00^ — -» r*-' V0-' 



O 
00 

ON 



o 

CO 

c 
o 



0) 

o 
a 



CO 
T3 

bO 

C 

•H 

cn 

CO 

•H 

B 

m 

0) 
CO 

•H 
O 

u 



ON 

a; 

r-l 
CO 



— V ip «? IP n 

— in O in O «n 



CN m OCT 



r^A «<« ^0— » vn-* tp^ "T^ v-» -r 

CTir Oin <Tin «in CTin t— o— O-^ 

O) ■n V n CTv CM^ mv «— ^— — CT n 

om om — in — m — in — ct ^ct nr* mo 



a»0 •fO — O ^"0 ^O coo fNP* CTr* 

Cn «Pn <Pn «V — n <NV CTO — O 

cNm — m vm — m nm c>*in — ct — ct 



I V) m Ul — I J — i-i %■* r- ^ .^ 



8-» tf)«i^ ot«-> CT-* - -» — -* o— ;-* o— w 

u) o»— on — n po— CM— ipCT 6— O— U)-» 

o»CT mn con n^T o»n mn cn^ O— tp— ^ct 

— m c^m — m — m — m — m — m cmct — ct vp» 



d-* o-' — -» 0-* 



O V-* eo-^ rM-— » vp-^ O)'-' in— r<i-^ n^- t— ^ »•' -* 

O <NCT neo OCT mm n in ipct r-CT f*»m v— O— t-* q 

O vn VPCT mv P^n — T — 'W Qn CT^ m— n— r^CT z 

O nm — m wm — m mm om Om — m on vn np» — o 

O m** CM-" o— — -* c^»-» O) -> - -^ m-* «p-» - -• - -* m 

O on CTc« 0)f*> oo«M CTn on OJP* Q« ^n <nO p-O m-» 

O cin a>n — CT On mn On On On — n m— «n— nCT 

O nm om om «Nin om om oim o»m — m ctct vpct nr* x 

< 
a 

(>i ^ O-* — -* f^ -* w^ f^ •* Q -^ n— ' c« — ' n— m — • vp — co iw 

mn - m nn mco «pct mm om — ct p*ct cim «p— f»— o^ ^ 

\tt n on wn vpct nn a»n — n On CTn ctv O— oo— chct z 

CM m CH m om om cm m — m om cm m — m — m « n cm ct n p* ^ i*^ 



m m 

CN or 
IP m 



in -— IP -* CT ^ 

on n m on 

(NT m n CN T 

CN m — m mm 



o -* p* — ' m — 

o o o m CT p» 

CM T n n an 

cm — m o m 



o m 


m r* 


CO - 


T — 


^p n 


(N "«r 


n — 


IP — 


o m 


- m 


O n 


O CT 



r* CT 


• U 


CT f* 


n O 




CN a 


' 


X 




c 




o 


n 


c 


CN — 


o 


— n 


-i 


IP r- 


n z 




CN u 



226 



concentration shifted from areal (during thermal bar conditions) to vertical 
(during thermally stratified conditions). The cruise 4 correlation matrix 
(Table 7.10) indicates the high level of vertical change in that largest corre- 
lations were among depth, temperature, and pH. The vertical stratification also 
induced a larger range in most variables compared to cruises 1 to 3 (Table 7.11), 
This larger range primarily resulted from lower values of soluble nutrients as 
they were utilized by near-surface algal blooms, typical of mid-summer condi- 
tions (Moll et al. 1980). 

Limnological conditions during cruise 5 (8-21 September) remained thermally 
stratified throughout Lake Huron, but a small amount of fall cooling had 
occurred. The fall cooling was indicated by a slightly lower mean temperature 
(Table 7.12). The cruise 5 correlation matrix (Table 7.13) was similar to the 
cruise 4 results with many high correlations among depth, temperature, and pH. 
Similar to cruise 4, the largest correlations among these three variables were 
associated with depth (because of thermal stratification) and less areal varia- 
bility. 

Cruise 6 (26 October to 4 November) was conducted during conditions of late 
thermal stratification. These conditions were represented by a rapidly cooling 
epilimnion. The epilimnion cooled both from atmospheric processes and a deepen- 
ing thermocline which entrained hypolimnetic water. As hypolimnetic water was 
entrained into the expanding epilimnion, nutrient concentrations increased in 
the epilimnion. These rapidly changing conditions promote decoupling of the 
normal phytoplankton-to-nutrient relationships observed during the mid-summer. 
Furthermore, the thermocline was deep enough to place the majority of the water 
in Lake Huron into the epilimnion rather than the hypolimnion. The cruise 6 
correlation matrix (Table 7.14), and descriptive statistics (Table 7.15) reflect 
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the increase in dissolved materials in the epilimnion. Mass of nutrients was 
conserved throughout the lake and simply mixed from the hypolimnion upward into 
the epilimnion. But, the sampling program was somewhat biased toward sampling 
in the epilimnion, i.e., more samples were taken in the epilimnion than hypo- 
limnion. As a result, the descriptive statistics (Table 7.15) show a uniform 
increase in dissolved nutrient concentration from cruise 5 to cruise 6. 

The results from the southern Lake Huron 1974 study (Moll 1980) showed an 
annual seasonal cycle of soluble silica depletion. This cycle, originally 
defined by Schelske and Stoerraer (1971), was evident in the 1974 results by a 
large negative correlation between soluble silica and chlorophyll from late June 
through October. The 1980 results presented here do not yield these high nega- 
tive correlations. Rather, the correlations among most soluble nutrients and 
chlorophyll were generally small (statistically non-significant from 0.00) and 
negative. These results should not be taken to imply a change in nutrient- 
phytoplankton relationships, but rather more as a consequence of the experi- 
mental design. The large area of the Lake Huron survey includes several differ- 
ent water types (see Chapters Five and Six), each with a somewhat different 
nutrient regime. These nutrient regimes promote different phytoplankton assem- 
blages (Stoermer 1978) which have different nutrient requirements. Under these 
conditions, the lake-wide correlation matrix has limited inferences other than 
general trends such as an increase in variability with depth after the onset of 
thermal stratification. 
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Factor Analyses 

Factor analyses carry the correlation analyses further in that a few under- 
lying sources of covariance are sought among the variables sampled (Kim and 
Mueller 1978) • These sources of variation are exposed in several ways. 
One method is the calculation of communalities or the proportion of variation 
in one variable which is shared by all other variables included in the analysis. 
Additionally, factor analysis reduces the number of variables in a data set 
by composing factors or linear combinations of the original variables (Harman 
1967). The order in which the analysis extracts the individual factors indi- 
cates where the highest source of covariance lies. Thus, the first factor con- 
tains more covariance than the second which explains more than the third factor 
and so on. 

The cruise 1 factor analysis yielded five factors (Table 7.16). The first 
and major factor had high loadings^ for soluble silica, ammonia, pH, and con- 
ductivity (last two opposite in sign from first two). This factor was taken to 
represent the high silica-low conductivity water mass covering much of northern 
Lake Huron. The remaining four factors each had high loadings for only one or 
two variables. Temperature had high loadings on the fifth and least important 
factor. These results were indicative of the well-mixed, early spring condi- 
tions typical of homothermous water masses. Strong inter-variable relationships 
were lacking throughout Lake Huron in April 1980. This hypothesis was supported 
by the low communalities of the cruise 1 factor solution. 



^The reader is cautioned about the confusing use of the term loadings in this 
section. In a limnological context, loading is often used to indicate a high 
nutrient input or nutrient loading. In a statistical context, factor loading 
refers to the coefficient in a factor pattern matrix. Throughout the remainder 
of this chapter, the term loading(s) indicates a statistical score. High 
inputs of nutrients are referred to as nutrient loading(s). 
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The cruise 2 factor analysis produced five factors (Table 7.17) and rela- 
tively large communal i ties. The first and most important factor had high load- 
ings for pH, conductivity, ammonia, and soluble silica, with the last two oppo- 
site in sign to the first two variables. This was the same result as from the 
cruise 1 factor analysis. Again, this factor was interpreted as representing 
the high silica-low conductivity water covering most of the open northern Lake 
Huron. The second factor had high loadings for temperature, total phosphorus, 
and total filtered phosphorus. With the onset of thermal bar conditions during 
cruise 2, temperature became a more important influence than it was during 
cruise 1; the temperature communality increased between cruise 1 and 2. The 
remaining three factors had high loadings for only one variable, implying that 
no single underlying mechanism was controlling the soluble nutrients together. 

The cruise 3 factor solution yielded four factors and a large range in the 
value of the communalities (Table 7.18). With the development of the thermal 
bar well offshore, thermal stratification inshore of the bar accounted for a 
large portion of Lake Huron water. The temperature communality was large, while 
the communalities associated with the various forms of phosphorus (total, total 
filtered, and soluble) were small. The first factor had high loadings for phos- 
phorus (total and total filtered) and ammonia. Highest total phosphorus levels 
were found offshore of the thermal bar and in the vicinity of river mouths. 
This factor implied that several regions of uniformly high nutrients were found 
throughout Lake Huron, particularly offshore of the thermal bar. The second 
factor had high loadings for conductivity and soluble silica (opposite in sign). 
Most of northern Lake Huron remained covered by high silica-low conductivity 
water, while the mouth of Saginaw Bay had low silica-high conductivity water. 
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The factor analysis for cruise 4 produced only three factors • The first 
and clearly dominant factor was interpreted as representing variability due to 
depth (Table 7.19), This factor had high loadings for NH3, NO3, soluble silica, 
pH, and temperature, with the last two variables opposite in sign to the first 
three. These five variables were observed to change rapidly with depth during 
periods of thermal stratification. The second factor had high loadings for 
conductivity alone. This factor was taken to represent high conductivity water 
from Saginaw Bay mixing into open Lake Huron. During periods of thermal strati- 
fication, Saginaw Bay water has been observed to flow far out into southern Lake 
Huron along the surface (Moll et al. 1980). The third factor had high loadings 
for total phosphorus and total filtered phosphorus. This factor was interpreted 
as representing high phosphorus water in most of the nearshore zones of Lake 
Huron. 

The cruise 5 factor solution (Table 7.20) produced five factors. Similar 
to the cruise 4 results, the first factor represented variability associated 
with increasing depth. This factor had high loadings for temperature, pH, NO3, 
and soluble silica, with the first two opposite in sign to the last two. The 
fourth factor in the cruise 5 solution was similar to the second factor from 
cruise 4, with high loadings for conductivity. The interpretation of this 
factor was the same as before, representing high conductivity in Saginaw Bay 
water. Factors two, three, and five had high loadings for ammonia, Kjeldahl 
nitrogen, and total filtered phosphorus, respectively. These three factors were 
taken to represent different types of nearshore water enriched in a particular 
nutrient. 
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The cruise 6 factor solution was similar to the cruise five results. Five 
factors were extracted (Table 7.21), with the first factor representing the var- 
iability associated with depth. Factor three is the high conductivity factor, 
while factors two, four, and five have high loadings for or tho phosphorus- total 
filtered phosphorus, Kjeldahl nitrogen, and total phosphorus respectively. 
These three factors were again considered to represent nearshore waters which 
had high nutrient loadings. 

Discussion 

The statistical analyses presented above yielded three major ecological 
inferences about Lake Huron. First, it is an oligotrophic ecosystem with low 
concentrations of nutrients and small annual ranges in nutrient concentrations. 
Second, physical processes have a major effect on the distribution of water 
masses in Lake Huron, while biological processes affect the concentration of 
nutrients within water masses. Third, a major transition in limnological condi- 
tions occurs when the lake changes from thermal bar conditions to thermally 
stratified conditions. 

The first inference, that Lake Huron is an oligotrophic ecosystem, comes 
from the means, ranges, and variances. Small ranges and coefficients of varia- 
tion imply small extremes in nutrient conditions throughout the lake. Further- 
more, small ranges suggest an absence of high nutrient loadings which stimulate 
large algal crops. The ranges for the major nutrients in Lake Huron during 1980 
are small when compared to more eutrophic lakes, i.e.. Lake Erie (Burns 1976). 
The mean value of the nutrient concentrations throughout Lake Huron also exem- 
plifies oligotrophic conditions when compared to other lakes. However, 
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comparison of means among different ecosystems almost always implies comparison 
of means among different studies and/or analytical laboratories. Because 
analytical procedures vary among laboratories, results are rarely comparable 
without some biases. The comparison of scale-free ranges (coefficients of 
variation or relvariances) provides inferences among different lakes without 
bias. These statistics from Lake Huron underscore the oligotrophic character- 
istic of the ecosystem when compared to Lakes Michigan and Erie (Rockwell et al. 
1980, Burns 1976). 

The second inference, that nutrient distributions in Lake Huron were pri- 
marily controlled by physical events, was developed by considering the results 
of the correlation and factor analyses. These results show major associations 
among the conservative variables and many of the non-conservative nutrients. 
For example, temperature-nutrient correlations were large and negative for 
cruises 4 to 6. These correlations, which are commonly observed in many lakes, 
are the consequence of nutrient depletion by near-surface phytoplankton. 
Although the algae are responsible for the nutrient reduction, the physical 
process of thermal stratification separates the epilimnion from the hypolimnion. 
Thus, physical processes were responsible for the distribution of water masses 
while biological processes controlled the concentration of nutrients within each 
water mass. 

Other examples of this general mechanism were made more visible with the 
aid of the multivariate analyses. The factor analyses showed that the major 
water masses of Lake Huron had distinct characteristics which often persisted 
across several months. Although the concentration of many nutrients changed 
within these water masses during the biologically important months, the rela- 
tionships among the variables remained relatively constant so that the same 
water masses from several different cruises were identified by factor analyses. 
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This third inference, that Lake Huron underwent a major limnological change 
during the transition from thermal bar to thermally stratified conditions, was 
also derived from the factor analysis. The six factor analyses considered 
together show a clear change in limnological conditions between cruises 1 to 3 
and 4 to &• The first three cruises were conducted during homothermous and 
thermal bar conditions. These are periods of intense spatial and vertical 
mixing (Scavia and Bennett 1980). As a result, the factor analyses identified 
only one or two major factors that had straightforward ecological interpreta- 
tion. Major factors in Lake Huron during cruises 1 to 3 were primarily due to 
large geographic differences in water quality. The major factor observed was 
the high silica-low conductivity of northern Lake Huron. 

The factor solutions from cruises 4 to 6 showed the large influence of 
thermal stratification. Surface temperatures varied considerably among these 
three cruises. However, the first and by far most dominant factor in each of 
the three analyses was variation with depth. The results indicate that once 
thermal stratification becomes established, a series of events takes place to 
alter the vertical distribution of many variables. The transition from homo- 
thermous to stratified conditions occurred rapidly in Lake Huron between early 
June and mid- July. 
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CHAPTER EIGHT 
LONG-TERM TRENDS IN NUTRIENT CONCENTRATIONS IN LAKE HURON 
by Russell Moll, Ronald Rossmann, and William Y. B. Chang 

A major objective of the Great Lakes Intensive Surveillance Program is 
to identify long-term trends in the water quality of each of the Great Lakes. 
These long-term trends are one of the few mechanisms used to judge the overall 
impact man has made on the Great Lakes. Lake-wide measurements made within one 
year or across two or three adjacent years rarely indicate an overall trend in 
water quality. Only through analysis of data collected over an extended period 
(such as a decade or more) can more reliable general trends in water quality, 
fish stocks, etc. be identified. 

Most of the long-term changes observed in the Great Lakes have been a 
consequence of accelerated rates of eutrophication from anthropogenic nutrient 
loadings (Beeton 1969). The consequences of the increase in eutrophication have 
been shifts in species composition and abundance (Stoermer and Yang 1970, Davis 
1964, Christie 1974). Although the results of long-term species shifts are 
evident from the flora and fauna, the primary mechanism responsible for those 
shifts is changes in the dissolved nutrient pools (Schelske and Stoermer 1971). 
Because of this, the surveillance of long-term changes in the nutrient concen- 
trations is of primary importance in the Great Lakes. 

The 1980 study of Lake Huron presents an opportunity to extend the histori- 
cal database by another six years; the previous lake-wide study of Lake Huron 
was completed in 1974 (IJC 1976). Lake Huron, along with Lake Superior, is 
considered one of the most oligotrophic of the Great Lakes (IJC 1976). Nutrient 
loadings to the Lake Huron ecosystem are primarily confined to the southern 
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basin, and more specifically to Saginaw Bay and the extreme southern Ontario 
shore (Schelske je£ al . 1980, Bierman and Dolan 1981). Over the past five years, 
nutrient loads in these regions decreased considerably due to mitigating actions 
taken by local municipalities (Bierman jet al. 1982). As a result, long-term 
trends have changed over the past decade. The results presented below show 
curvilinear long-term trends for several variables over the last 30 years. 

The term "long-term trend" can include a variety of different processes 
depending on the specific trend of interest. This chapter deals with three 
types of long-term trends in Lake Huron: 1) changes in mean nutrient concentra- 
tions across a period of 10 years or more holding the effect of season constant, 
2) changes in the range of variables among years, and 3) changes in the areal 
extent of homogeneous water masses among years. The first type of trend is used 
most commonly to evaluate the overall status of a Great Lakes ecosystem. 
A simple comparison of means across years reveals whether or not nutrient con- 
centrations have changed, which in turn may cause species shifts. Extreme care 
must be taken to make this comparison for the same season of each year. Sea- 
sonal effects can have a major influence on mean nutrient concentrations and 
thus either mask significant long-term trends or give the false impression that 
a trend exists (Richards 1981a). The second type of long-term trend is used to 
evaluate the magnitude of seasonal fluctuations of non-conservative variables. 
A large annual range of dissolved nutrients indicates high nutrient loadings 
stimulating large algal crops (Vollenweider et al. 1974). A large annual range 
can also be the result of a high degree of spatial heterogenity (see Chapter 
Five). The third type of annual trend follows from above in that the areal 
extent of high-nutrient or low-nutrient water masses can be compared across 
years. This type of comparison is used to assess the extent that point-source 
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nutrient inputs affect the lake. This chapter deals primarily with the first 
two types of long-term trends in Lake Huron: 1) comparison of means across 
years and 2) changes in annual ranges. 

Historical Database 

The modern historical database for Lake Huron begins in 1954 (Ayers et^ al . 
1956) and proceeds through 1980 (see Chapter Three). The historical database 
includes studies in the following years: 1954 (Ayers et al. 1956); 1956 (Allen 
1964); 1966, 1967 (Kramer 1968); 1968, 1969 (CCIW 1968, 1969); 1970 (Schelske 
and Roth 1973); 1971, 1972, and 1973 (CCIW STAR data); and 1974 (Schelske e£ al. 
1980). Most of these studies were conducted in the mid-summer through late fall 
(Table 8.1); but the 1971, 1972, 1973, and 1974 studies also included spring 
cruises. Because the nutrient dynamics and biological characteristics of this 
lake are considerably different in spring and summer, comparisons of means 
across years for the purpose of identifying long-term trends are made for each 
season (Table 8.1). Summer is not a good season for annual long-term trend 
analysis, because during this period thermal stratification is already well 
developed (Boyce 1974). As a result of this, surface nutrients are depleted by 
algae (Moll et al. 1980), and vertical nutrient profiles are markedly different. 
The experimental design of a research cruise conducted during the summer can 
very much affect the mean and range of many variables. 

The types of variables measured over the past 30 years have changed along 
with analytical methodology. Table 8.2 shows that relatively few dissolved 
nutrient analyses were conducted in the late 1950s (mostly silicon). By 1968, 
a fairly complete suite of analyses were conducted as a routine basis on most 
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research cruises. Analytical methodology changed considerably across the 
30 years of study. However, corrections for bias cannot easily be made across 
years. The long-term trends developed in this chapter have not adjusted the 
data for inter-lab or inter-year differences. 

Statistical Methods 

A variety of methods have been used to evaluate long-term trends in the 
Great Lakes. These methods frequently include graphical techniques, determin- 
istic models, and statistical models. The techniques presented below are pri- 
marily statistical models. The major reasons behind the use of statistical 
models is that they take into account the stochastic nature of the ecosystem 
and provide probability statements to attach to inferences concerning long-term 
trends. 

The statistical procedures used to determine long-term trends were part of 
the University of Michigan statistical package MIDAS. This is a large, inter- 
active data management and analysis package which provides software for univari- 
ate and multivariate analyses (Fox and Guire 1976). The protocol in analyzing 
the long-term trends was as follows: A database with 12 years of data was 
created. Data were available from 1954, 1956, 1966, 1967, 1968, 1969, 1970, 
1971, 1972, 1973, 1974, and 1980. Eight analytical variables were included in 
this database: temperature, conductivity, chloride, sulfate, soluble silica, 
nitrate-nitrogen, dissolved reactive phosphate, and total phosphorus. 

Each variable was analyzed to determine whether or not it was normally dis- 
tributed within each year. This was accomplished by examining the skewness and 
kurtosis of each variable and comparing the mean to the median. Most of the 
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variables had relatively low skew, generally below 2.0. However, some parame- 
ters in the data showed deviation from normality, especially soluble phosphate. 
In those cases, comparisons between medians and means were made. A large dif- 
ference between medians and means suggests a lack of normality and is evidence 
for the need to use non-parametric analysis. The results showed little differ- 
ence between the means and medians. Although in a few cases the data were not 
normally distributed, the conclusion from this analysis was that the data were 
close enough to being normally distributed for the use of parametric analyses. 

Analyses for long-term trends in Lake Huron were conducted as a two-step 
statistical procedure. The first hypothesis tested was to determine if the 
variable(s) changed among years. This hypothesis was tested using ANOVA with 
the null hypothesis that there is no difference among years versus the alterna- 
tive hypothesis that at least two different years are statistically different. 
The statistical procedure used was a one-way ANOVA with unequal cell sizes 
(Winer 1971). 

After confirming that there were significant differences across years, 
the analysis proceeded to determine if a linear or curvilinear trend was present 
across years. The null hypothesis tested for each variable was that there was 
no linear trend across years (slope is 0.0) versus the alternative that there 
was a significant slope (Netter and Wasserman 1974). The linear and curvilinear 
trends were analyzed by a least squares stepwise polynomial procedure. In this 
technique, a linear model was fit to the data followed by a quadratic model. 
An F-test was conducted to determine if the quadratic fit significantly improved 
the linear fit. The final results of the ANOVA and regression analyses were 
combined into one ANOVA table with a linear decomposition of the treatment sums 
of squares for the regression model (Netter and Wasserman 1974). In some cases, 
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however, polynomial fitting was extended beyond the quadratic term if a signifi- 
cant gain in R^ resulted from the inclusion of higher order terms, and nonzero- 
regression coefficients were obtained for these terms after such an inclusion. 
The best polynomial equation for describing a data set was then chosen to repre- 
sent the curvilinear trend of the data. The results from these polynomial fits 
are presented along with linear fits in the figures which follow. For the data 
considered, the significance of nonlinearity as indicated by quadratic terms 
is not different from that shown by the inclusion of the higher order terms. 
The authors recognize that significant higher order polynomial fits to the data 
are not useful for predicting trends but are useful for describing changes 
within parameters over a designated period of time. Such curves are useful for 
describing long-term variations of parameters. Superimposed on these variations 
may be a long-term linear trend. 

Long-term trends were also analyzed by comparing homogeneous water masses 
across years. Using the same analytical techniques as in Chapter Six, homo- 
geneous water masses were identified from surface samples of one research cruise 
for each year. The statistical procedures used were cluster analysis and dis- 
criminant functions analysis. Once the homogeneous water masses were identified 
for each year, water masses of approximately the same location were compared 
across years for long-term trends. The long-term trends using homogeneous water 
masses were calculated using ANOVA and regression techniques similar to the 
long-term trends based on all the data from each research cruise. 
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Annual Means and Ranges 

Table 8.3 presents the annual means and standard deviations for eight vari- 
ables sampled in twelve different years spanning 1954 to 1980. Although some 
long-term trends are evident from these means, this type of analysis is inappro- 
priate for identifying trends. Composite analysis of annual data has two major 
short-comings. The first is that research cruises conducted in each year are 
unlikely to cover the same seasonal span. The temperature data in the Lake 
Huron historical database indicate that this has occurred in the Lake Huron 
sampling. Conceptually, long-term trends in temperature are not realistic. 
Thermal budgets of the Great Lakes are primarily controlled by insolation, wind 
stress, and water column depth (Sundaram 1973). Seasonal fluctuations in the 
thermal regime of Lake Huron could be expected, but on the average the thermal 
cycle should not have a long-term trend (Sundaram 1973). The annual differences 
observed in mean temperatures in Table 8.3 are a reflection of the experimental 
design changing among years. Warm annual averages with small standard devia- 
tions (1967) indicate sampling programs concentrating in the summer months and 
collecting most of the samples from the epilimnion. Cold annual temperatures 
with large standard deviations (1980) reflect a broader sampling program, both 
seasonally and spatially. 

The second problem associated with identifying long-term trends from annual 
means is the areal inconsistency of the sampling programs across the years. 
Some studies have been confined mostly to the nearshore (1970), while others 
include the entire lake (1968, 1969, 1971, 1972, 1973, and 1980). An example 
of this effect can be seen in Table 8.3. The mean 1966 conductivity reading 
is 135 umho/cm, which is considerably lower than the other annual means. 
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This lower mean is more a function of the location sampled (near Manitoulin 
Island) than an unusual year. The spatial and seasonal incompatibilities across 
years prevented any analysis for long-term trends using annual means. 

Annual ranges for nine variables were fairly similar between 1974 and 1980 
(Table 8.4). These two sampling programs covered the same seasons (April to 
November) and all of offshore Lake Huron. Changes in the annual range suggest 
changes in processes occurring within Lake Huron. For example, a large chloro- 
phyll range indicates large algal blooms occurred sometime during the year. 
Large nutrient ranges imply high nutrient loadings at some time during the year 
and periods of nutrient depletion at other times. 

The ranges in Table 8.4 provide several ecological inferences in regard to 
differences between 1974 and 1980. The temperature ranges were similar for 1974 
and 1980 at 20.6°C and 21.9°C, respectively. This suggests the two surveys were 
conducted across similar seasons. The conductivity and chloride ranges were 
different between 1974 and 1980. The 1974 study tended to sample more in the 
mouth of Saginaw Bay and thus had higher conductivity and chloride maximum 
values. If one assumes that similar station patterns and densities were sampled 
in 1974 and 1980, the higher ranges in the 1974 study for all dissolved nutri- 
ents (chloride, sulfate, silica, and soluble phosphate) and chlorophyll could 
suggest a reduction in nutrient loadings and algal growth in 1980 as compared 
to 1974. However, different stations were sampled in these two years, and it is 
difficult to infer long-term trends by comparing ranges for only two years. 
Annual ranges in this case are of little use; only if more data on annual ranges 
are collected from the same locations can these ranges be useful in inferring 
the nutrient conditions in Lake Huron. 
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Table 8.4. Minimum, maximum, and ranges from 1974 and 1980 Lake Huron surveys, 









1974 






1980 




Variable 


Min. 


Max. 


Range 


Min. 


Max. 


Range 


Temperature 


1.0 


21.6 


20.6 


0.7 


22.6 


21.9 


Conductivity 


89. 


322. 


233. 


86. 


253. 


167. 


Chloride 


1.3 


20.2 


18.9 


1.3 


9.0 


7.7 


Sulfate 


3.5 


27.0 


23.5 


2.8 


21.2 


18.4 


Silica 


0.1 


3.5 


3.4 


0.6 


2,7 


2.1 


Nitrate 


0.181 


0.838 


0.657 


.033 


.765 


.732 


Sol. Phosphorus 


<0.1 


17.5 


17.4 


0.1 


9.0 


8.9 


Total Phosphorus 


<0.1 


42.0 


41.9 


2.0 


58.1 


56.1 


Chlorophyll 


0.1 


16.3 


16.2 


0.8 


10.0 


9.2 
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Long-Term Trends 

The effect of seasonality on nutrient concentrations in Lake Huron is large 
and statistically significant (Chapter Four, Appendix A, Schelske et al> 1980). 
This effect can be removed either statistically or by using the database for 
each season (Richards 1981b). This latter approach was taken in the trend anal- 
ysis; spring and summer data were analyzed individually. The spring data 
include those collected in April and May while the summer data include those 
from the cruises during July, August, and September. 

The initial null hypothesis tested for each of the eight variables sepa- 
rately was that there was no difference among years versus the alternative that 
significant differences exist among years. This hypothesis was tested with one- 
way ANOVA with unequal sample sizes (Winer 1971). The results (Table 8.5) show 
that all variables were significantly different among years at the less than 
0.1% probability limit. The statistical significance of temperature varying 
across years implies differences in experimental design for the different 
studies. Seasonal factors should be fairly small, because this analysis consid- 
ered studies from only April and May. The differences in temperature across 
years suggests that some studies may have taken a larger proportion of samples 
from the nearshore zone than other studies. The total sums of squares are 
compared in Table 8.6 for cases where the Lake Huron historical database was 
analyzed using all the data, the spring data only, and one water mass from the 
spring months only. The total sums of squares were reduced significantly when 
only the spring months were analyzed compared to the full-year data. Therefore, 
the long-term trends developed from the spring data have excluded a large por- 
tion of variation due to seasonality. 
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Table 8. 5- ANOVA tables for Lake Huron spring data. Each analysis included 
all data collected from the months of April and May. The ANOVA 
tables are decomposed into linear and quadratic trends. 







Sum of 


Degrees of 


Mean 


F- 




Squares 


Freedom 


Square 


Statistic 


TEMPERATURE 










Years 


569.3 


4 


142.3 


43.4 


Linear Trend 


330.2 


1 


330.2 


98.4 


Lack of fit 


239.1 


3 


79.7 


24.2 


Quadratic 


1.7 


1 


1.7 


0.5 


Lack of fit 


237.4 


2 


118.7 


36.0 


Error 


9558.4 


2917 


3.3 




Total 


10128 


2921 







CONDUCTIVITY 



Years 

Linear trend 
Lack of fit 
Quadratic 
Lack of fit 

Error 

Total 



204070 

136880 

67190 

9238.2 
57951.8 
1414100 
1618200 



4 


51017 


91.3 


1 


136880 


234,1 


3 


22396.7 


40.1 


1 


9238,2 


16.3 


2 


28975.9 


51.8 


2530 


559.0 




2534 







CHLORIDE 



Years 

Linear Trend 
Lack of fit 
Quadratic 
Lack of fit 

Error 

Total 



146.5 
68.3 
78.2 
12.7 
65.5 
5162.9 
5309.4 



4 
1 
3 
1 
2 
1313 
1317 



36.6 


9.3 


68.3 


17.4 


26.1 


6.6 


12.7 


3.2 


32.8 


8.3 


3.93 





SULFATE 



Years 

Linear trend 
Lack of fit 
Quadratic 
Lack of fit 

Error 

Total 



79.7 

0.1 

79.6 

79.6 



5277. 
5357. 



2 
1 
1 
1 

909 
911 



39.8 

0.1 

79,6 

79.6 

5.8 



6.9 

0.0 

13.7 

13.7 



(Continued) 
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Table 8.5. (Concluded). 



Sum of 
Squares 



Degrees of 
Freedom 



Mean 
Square 



Statistic 



SOLUBLE SILICA 



Years 


40.2 


Linear trend 


0.2 


Lack of fit 


40.0 


Quadratic 


37.1 


Lack of fit 


2.9 


Error 


238.8 


Total 


279.0 



4 


10-0 


98.6 


I 


0.2 


1.7 


3 


13.3 


133.3 


1 


37.1 


360.0 


2 
2346 


1.5 
0.1 


14.5 


2350 







NITRATE-NITROGEN 



Years 


0.470 


3 


Linear Trend 


0.370 


1 


Lack of fit 


0.100 


2 


Quadratic 


0.0999 


1 


Lack of fit 


0.0001 


1 


Error 


3.49 


2011 


Total 


3.96 


2014 



0-157 


90.2 


0.370 


213.2 


0.050 


28.9 


0.0999 


57.6 


0.0001 


0.1 


0.00173 





SOLUBLE PHOSPHATE 



Years 


1113.7 


Linear trend 


3.1 


Lack of fit 


1110.6 


Quadratic 


251.4 


Lack of fit 


859.2 



Error 
Total 



7870.5 
8984.2 



4 


278.4 


80.4 


1 


3.1 


0.8 


3 


370.2 


105.8 


1 


251.4 


71.8 


2 


429.6 


122.7 


2272 


3.5 




2276 







TOTAL PHOSPHORUS 



Years 


1453 


Linear trend 


2.1 


Lack of fit 


1450.9 


Quadratic 


940 


Lack of fit 


510.9 



Error 
Total 



37177 
38630 



4 


363-3 


22.5 


1 


2.1 


0-1 


3 


483.6 


29-9 


1 


940 


57.4 


2 


255.5 


15.8 


2299 


16.2 




2303 
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Spring - All Data 

Least-squares regression was used to determine whether or not a significant 
linear trend existed for any of the eight variables between the years 1971 and 
1980. Temperature, conductivity, chloride, and nitrate had a significant linear 
trend (0.1 level of significance) (Table 8.7 and Fig. 8.1). Sulfate, soluble 
reactive silica, soluble reactive phosphorus, and total phosphorus had no sig- 
nificant long-term linear trend, but each had a significant curvilinear trend 
(Fig. 8.1). Mean concentrations of sulfate decreased between 1974 and 1977 and 
increased after 1977. Those for silica increased between 1971 and 1976 and 
decreased after 1976. Mean concentrations of soluble reactive and total phos- 
phorus increased between 1971 and 1976 and decreased after 1976. 

Spring - One Homogeneous Water Mass 

The problem of inconsistent experimental designs across years is both 
seasonal and spatial. The approach of using only spring data removes many 
of the seasonal differences, but it does not remove any spatial differences. 
The technique of identifying homogeneous water masses developed in Chapter Six 
was applied to the historical database. The objective was to generate one or 
more homogeneous water masses that were compared across years. The technique 
used was cluster analysis of surface samples to identify the water masses. 
Discriminant functions analysis was then used to confirm whether or not the 
identified clusters were statistically distinct from one another within each 
year. The final result was that one water mass per year was identified as 
representing open southern Lake Huron water quality. In this manner, spatial 
differences in the sampling program were reduced across the study years. 
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The identification of homogeneous water masses was less than perfect in 
that a few years actually had their water masses located more toward the north- 
ern end of Lake Huron (1966). Table 8.6, nonetheless, shows a considerable 
reduction in the sums of square between the spring month analyses and the one 
water mass analyses. The objective of decreasing the overall variance in the 
long-term trend analysis was achieved. 

The same hypotheses were tested for the one water mass data as for the 
spring month data. First, do variables differ significantly across years 
(ANOVA)? Second, is there a long-term linear trend (Regression)? The results 
of these analyses in Table 8.8 show that the inferences using the one water mass 
data are almost the same as using all of the spring data. All eight variables 
changed significantly across the 1971 to 1980 period (Table 8.8). 

The regression analysis using the one water mass data results were somewhat 
similar to the results using all of the spring data (Table 8.7). Nitrate and 
temperature continued to show an increasing trend, and conductivity showed a 
decreasing trend for both data sets. Unlike all of the spring data, one water 
mass data for sulfate and total phosphorus showed an increasing trend. Only 
chloride showed a major change in the regression coefficient using the one water 
mass data. In this case, a decrease was noted for all spring data, while an 
increase was found for the spring one water mass data. The difference may 
perhaps be due to the inclusion of data collected from the nearshore zone or 
from an emphasis on southern Lake Huron data in the second case. Six of the 
eight variables showed significant linear trends over the 10 years, with soluble 
phosphate and soluble reactive silica displaying no linear trend (Fig. 8.2). 
However, both displayed a curvilinear variation (Fig. 8.2). Both increased 
between 1971 and 1976 and decreased after 1976. All curvilinear fits displayed 



267 



Table 8. 8, ANOVA tables for Lake Huron spring data using one homogeneous 
water mass from each cruise. The ANOVA tables are decomposed 
into linear and quadratic trends. 



Sum of 
Squares 



Degrees of 
Freedom 



Mean 
Square 



F- 
Statistic 



TEMPERATURE 



Years 


181.9 


4 


Linear Trend 


128.3 


1 


Lack of fit 


53.6 


3 


Quadratic 


8.3 


1 


Lack of fit 


45.3 


2 


Error 


943.8 


504 


Total 


1125.6 


508 



45.5 


24.3 


128.3 


68.6 


17.9 


9.6 


8.3 


4.4 


22.7 


12.1 


1.9 





CONDUCTIVITY 



Years 40263 

Linear trend 1195 

Lack of fit 39068 

Quadratic 3956 

Lack of fit 35112 

Error 49702 

Total 89965 



4 


10066 


86.9 


1 


1195 


10.3 


3 


13023 


112.4 


1 


3956 


34.1 


2 


17556 


151.5 


429 


116 • 




433 







CHLORIDE 



Years 

Linear Trend 
Lack of fit 
Quadratic 
Lack of fit 

Error 

Total 



9.65 
1.70 
7.95 
5.98 
1.97 
28.3 
38.0 



4 
1 
3 
1 
2 
204 
208 



2.41 


17.4 


1.70 


12.1 


2.65 


18.9 


5.98 


42.7 


0.99 


7.1 


0.14 





SULFATE 



Years 

Linear trend 
Lack of fit 
Quadratic 

Lack of fit 

Error 

Total 



18. 

6e 

11. 
11. 

90. 
108. 



2 
1 
1 
1 

151 
153 



9.2 


15.4 


6.5 


10.8 


11.9 


19.8 


11.9 


14.5 



0.60 
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(continued) 



Table 8.8. (Concluded). 



Sum of 
Squares 



Degrees of 
Freedom 



Mean 
Square 



F- 
Statistic 



SOLUBLE SILICA 



Years 


7.73 


Linear trend 


0.09 


Lack of fit 


7.64 


Qxiadratic 


7.62 


Lack of fit 


0.02 


Error 


10.8 


Total 


18.6 



4 
1 
3 
1 
2 
428 
432 



1.93 


76.3 


0.09 


3.6 


2.55 


102.0 


7.62 


304.8 


0.01 


0.4 


0.025 





NITRATE -NITROGEN 



Years 


0.48 


3 


Linear Trend 


0.47 


1 


Lack of fit 


0.1 


2 


Quadratic 


0.0005 


1 


Lack of fit 


0.0005 


1 


Error 


1.18 


387 


Total 


1.66 


390 



0.16 


53.1 


0.47 


156.7 


0.05 


16.7 


0.0005 


0.2 


0.0005 


0.2 


0.003 





SOLUBLE PHOSPHATE 



Years 


1509.2 


Linear trend 


1.2 


Lack of fit 


1508.0 


Quadratic 


133.0 


Lack of fit 


1375.0 


Error 


5968.0 


Total 


7477.3 



4 


377.3 


25.9 


1 


1.2 


0.1 


3 


502.7 


34.4 


1 


133.0 


9.1 


2 


687.5 


47.1 


410 


14.6 




414 







TOTAL PHOSPHORUS 



Years 


415.8 


Linear trend 


142.6 


Lack of fit 


273.2 


Quadratic 


56.9 


Lack of fit 


216.3 


Error 


1737. 


Total 


2152.8 



4 


104.0 


25.5 


1 


142.6 


34.8 


3 


91.1 


22.2 


1 


56.9 


13.9 


2 


108.2 


26.4 


426 


4.1 




430 
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in spring for one homogeneous water mass resemble those shown in spring for the 
entire lake. 

Summer 

The same approach and hypotheses used in analyzing the spring data were 
utilized for the summer data. The results showed that all variables were sig- 
nificantly different among years at the less than 0.1% probability limit 
(Table 8.9). The statistical significance of temperature varying across years 
for the summer data was again attributed to sample design as discussed in the 
preceding section. The analysis by season has greatly reduced variation in the 
data, as shown by the reduction of the sum of squares (Table 8.6). 

Least-squares regression was also used to determine if a significant linear 
trend existed for the 1954 to 1980 period. All eight variables displayed a sig- 
nificant linear trend (Table 8.7, Fig. 8.3). Conductivity displayed a signifi- 
cant 2nd order variation (Fig. 8.3). It increased between 1954 and 1970 and 
decreased after 1970. Soluble reactive phosphorus displayed a significant 3rd 
order variation. It decreased between 1966 and 1970, increased between 1970 and 
1977, and decreased after 1977. 

Summer - One Homogeneous Water Mass 

The analysis of nutrient data by season can remove many of the seasonal 
differences, but it does not reduce any spatial heterogeneity. In order to 
select one or more homogeneous water masses that could be compared across years , 
the technique of identifying homogeneous water masses used in connection with 
the spring data was applied to reduce the heterogeneity in the summer data set. 
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Table 8.9. ANOVA tables for Lake Huron summer data. Each analysis included 

all data collected from the months of July, August, and September. 
The ANOVA tables are decomposed into linear and quadratic trends. 







Sum of 


Degrees of 


Mean 


F- 




Squares 


Freedom 


Square 


Statistic 


TEMPERATURK 










Years 


6848.1 


11 


622.6 


19.5 


Linear Trend 


1066.0 


1 


1066.0 


33.4 


Lack of fit 


5782.1 


10 


578.2 


18.1 


Quadratic 


83.6 


1 


83.6 


2.6 


Lack of fit 


5698.5 


9 


633.2 


19.8 


Error 


137110 


4304 


31.9 




Total 


143960 


4315 






CONDUCTIVITY 










Years 


777080 


11 


70643 


334.3 


Linear trend 


281100 


1 


281100 


1330.3 


Lack of fit 


495980 


10 


49598 


234.7 


Quadratic 


326770 


1 


326770 


1546.5 


Lack of fit 


169210 


9 


18801.1 


89.0 


Error 


851800 


4031 


211.3 




Total 


1709100 


4042 






CHLORIDE 










Years 


458.8 


8 


57.4 


42.4 


Linear trend 


245.7 


1 


245.7 


175.5 


Lack of fit 


213.1 


7 


30.4 


21.7 


Quadratic 


15.0 


1 


15.0 


10.7 


Lack of fit 


198.1 


6 


33.0 


23.6 


Error 


2474.7 


1829 


1.4 




Total 


2933.6 


1837 






SULFATE 










Years 


2250.1 


7 


321.4 


95.6 


Linear trend 


876.8 


1 


876.8 


257.9 


Lack of fit 


1373.3 


6 


228.9 


67.3 


Quadratic 


120.7 


1 


120.7 


35.5 


Lack of fit 


1252.6 


5 


250.5 


73.7 


Error 


3699.2 


1100 


3-4 




Total 


5949.3 


1107 










Ol-K 




(Continued) 



Table 8,9. (Concluded) • 



Sum of 
Squares 



Degrees of 
Freedom 



Mean 
Square 



F- 
Statistic 



SOLUBLE SILICA 



Years 


553.6 


Linear trend 


365.1 


Lack of fit 


188.50 


Quadratic 


123.3 


Lack of fit 


65.2 


Error 


650.0 


Total 


1203.6 



11 


50.3 


264.2 


1 


365.1 


1825.5 


10 


18.8 


94.2 


1 


123.3 


616.5 


9 


7.2 


36.2 


12 


.2 





NITRATE -NITROGEN 



Years 


2.7926 


6 


Linear trend 


1.6642 


1 


Lack of fit 


1.1284 


5 


Quadratic 


0.0052 


1 


Lack of fit 


1.1232 


4 


Error 


4.3083 


2075 


Total 


7.1009 


2081 



0.46544 


224.2 


1.6642 


792.5 


0.2257 


107.5 


0.0052 


2-5 


0.2808 


133.7 


0.0021 





SOLUBLE PHOSPHATE 



Years 


1346.9 


Linear trend 


36.3 


Lack of fit 


1310.6 


Quadratic 


116.5 


Lack of fit 


1194.1 


Error 


3603.8 


Total 


4950.7 



9 


149.7 


101.9 


1 


36.3 


24.2 


8 


163.8 


109.2 


1 


116.5 


77.7 


7 


170.6 


113.7 


2453 


1.5 




2462 







TOTAL PHOSPHORUS 



Years 


3438.3 


Linear trend 


275.4 


Lack of fit 


3162-9 


Quadratic 


18-6 


Lack of fit 


3144.3 


Error 


28610 


Total 


32048 



6 


573.05 


43.3 


1 


275.4 


20.9 


5 


632-6 


47.9 


1 


18-6 


1.4 


4 


786.1 


59.6 


2162 


13.2 




2168 
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However, spatial heterogeneity occurs not only horizontally but also vertically 
in summer. Vertical differences often are larger and more distinct than the 
horizontal differences. In order to identify a water mass which is homogeneous 
both vertically and horizontally, an epilimnetic portion of the homogeneous 
water mass was selected for use in the analysis. The epilimnion is defined 
as the surface to 5 m in depth. The comparison shown below, however, included 
the results from both one water mass and only the epilimnetic portion of the 
water mass. The result showed a considerable reduction of the sums of squares 
between the summer month analysis and the one water mass analysis and between 
the one water mass and the epilimnetic portion of the water mass (Table 8.6). 
This indicates that the technique applied has been successful in reducing the 
variation in the data. However, it is also noted that such a selection has 
reduced the sample size in the data set of the epilimnetic portion of the water 
mass. 

The hypotheses used for the spring data were tested for the summer one 
water mass data. All eight variables changed significantly between 1954 and 
1980 (Table 8.10). The regression coefficients derived from both the one water 
mass and the epilimnetic portion of the water mass were quite similar to those 
from the whole summer data set (Table 8.7). For the one water mass subsettiag 
of the data, all but total phosphorus displayed significant linear trends. 
The signs of the regression coefficients generated from either one water mass 
or the epilimnetic portion of the water mass were identical to those from all 
of the stations' summer data. The magnitude of the coefficients was also very 
similar between these data sets. Like those of spring, the trends are small 
but statistically significant. These changes were smaller than those reported 
from Lakes Erie and Michigan (Bee ton 1969). For the epilimnion one water mass 
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Table 8.10. ANOVA tables for Lake Huron summer data using one homogeneous 
water mass from each cruise. The ANOVA tables are decomposed 
into linear and quadratic trends. 







Sum of 


Degrees of 


Mean 


F- 




Squares 


Freedom 


Square 


Statistic 


TEMPERATURE 










Years 


5875.0 


11 


534.1 


15.1 


Linear trend 


888.7 


1 


888.7 


25.1 


Lack of fit 


4986.3 


10 


498.6 


14.1 


Quadratic 


1445.7 


1 


1445.7 


40.8 


Lack of fit 


3540.6 


9 


393.4 


11.1 


Error 


35046 


990 


35.4 




Total 


40921 


1001 






CONDUCTIVITY 










Years 


215990 


11 


19635 


515.9 


Linear trend 


141780 


1 


141780 


3721.3 


Lack of fit 


74210 


10 


7421 


194.8 


Quadratic 


12640 


1 


12640 


331.8 


Lack of fit 


61570 


9 


6841. 


179.6 


Error 


35700 


938 


38. 




Total 


251690 


949 






CHLORIDE 










Years 


65.4 


8 


8.2 


18.1 


Linear trend 


5.2 


1 


5.2 


11.6 


Lack of fit 


60.2 


7 


8.6 


19.1 


Quadratic 


1.0 


1 


1.0 


2.2 


Lack of fit 


59.2 


6 


9.9 


22 


Error 


180.2 


398 


0.45 




Total 


245.6 


406 






SULFATE 










Years 


835.1 


7 


119.3 


81.0 


Linear trend 


262.1 


1 


262.1 


178.3 


Lack of fit 


573.0 


6 


95-5 


65.0 


Quadratic 


7.2 


1 


7.2 


4.9 


Lack of fit 


565.8 


5 


113.2 


77.0 


Error 


363.8 


247 


1.47 




Total 


1199 


254 














(continued) 
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Table 8.10. (Concluded). 



Sum of 
Squares 



Degrees of 
Freedom 



Mean 
Square 



F- 
Statistic 



SOLUBLE SILICA 



Years 


128.5 


Linear trend 


94.8 


Lack of fit 


33.7 


Quadratic 


17.1 


Lack of fit 


16.6 


Error 


132.8 


Total 


261.3 



11 


11.7 


74.1 


1 


94.8 


592.5 


10 


3.4 


21.1 


1 


17,1 


106.9 


9 


1.8 


11.5 


842 


0.16 




853 







NITRATE -NITROGEN 



Years 


.526 


Linear trend 


.271 


Lack of fit 


.255 


Quadratic 


0.0 


Lack of fit 




Error 


.803 


Total 


1.33 



6 
1 
5 
1 

499 
505 



0.877 


54.5 


0.271 


169.4 


0.051 


31.9 


0.0 


0.0 



0.0016 



SOLUBLE PHOSPHATE 



Years 


395.6 


Linear trend 


30.8 


Lack of fit 


364.8 


Quadratic 


23-1 


Lack of fit 


341.7 


Error 


692.0 



Total 



1087.7 



9 


44.0 


32.2 


1 


30.8 


22.0 


8 


45.6 


32.6 


1 


23.1 


16.5 


7 


48.8 


34.9 


507 


1.4 




516 







TOTAL PHOSPHORUS 



Years 


987.3 


Linear trend 


0.1 


Lack of fit 


987.2 


Quadratic 


2.6 


Lack of fit 


984.6 


Error 


2595.4 


Total 


3582.7 



6 


164.6 


31.6 


1 


0.1 


0.0 


5 


197.4 


38.0 


1 


2.6 


0-5 


4 
498 


246.2 
5.2 


47.3 


504 
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of summer, linear trends were significant for temperature, conductivity, chlo- 
ride, sulfate, dissolved silica, and nitrate (Table 8.7, Fig. 8.4). Four vari- 
ables were best described by curvilinear fits (Fig, 8.4). Temperature decreased 
between 1954 and 1966 and increased after 1966. Sulfate increased between 1956 
and 1962, decreased between 1962 and 1975, and increased after 1975. Soluble 
reactive phosphorus decreased between 1966 and 1970, increased between 1970 and 
1977, and decreased after 1977. Total phosphorus decreased between 1968 and 
1977, and increased after 1977, 
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Conclusions 

Long-term changes in the chemistry of Lake Huron waters have been identi- 
fied using the results of twelve studies spanning a 26-year period of observa- 
tions. As might be expected of a dynamic system lake. Lake Huron long-term 
trends are complex. For roughly one-half of the parameters investigated, 
changes over the 26 years were curvilinear or oscillatory. Of these, one-half 
had a linear trend superimposed over the oscillations. For each of the various 
ways we manipulated the data base, sulfate and nitrate were found to increase. 
For the majority of cases, dissolved silica and chloride decreased. Sulfate 
increased at a rate of 0.05 to 0.13 mg/L/yr or 0.3 to 0.8% of the 1980 mean 
concentration per year. Nitrate increased at a rate of 0.0030 to 0.0089 mg/L/yr 
or 1.1 to 3.1% per year. Dissolved silica decreased at a rate of 0.040 to 
0.049 mg/L/yr or 2.8 to 3.6% of the 1980 mean per year. Chloride decreased at a 
rate of 0.02 to 0.06 mg/L/yr or 0.4 to 1.2% of the 1980 mean per year. All four 
parameters should continue to be monitored. 

For those parameters which showed reversals in trend over the period of 
observation, reversals consistently occurred in the mid-1970s. It is unfortu- 
nate that data were not available for this critical period of time. Since the 
mid-1970s, dissolved reactive silica and soluble reactive phosphorus have 
decreased, and nitrate and sulfate have increased. Changes in long-term trends 
could be or have been attributed to factors such as changes in sampling and 
analytical methods, reduction in anthropogenic inputs from the major point 
sources, and differences in biological utilization of these elements. Despite 
the changes which have occurred in Lake Huron between 1954 and 1980, Lake Huron 
remains an oligo trophic lake. 
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APPENDIX A 
COMPARISON OF SURFACE/ EPILIMNION TO LOWER LAYER/ HYPOLIMNION 



A-1 



Table A-1 . Comparison of surf ace/ epilimnion total phosphorus with lower layer/ 
hypolimnion total phosphorus. Average values which differ between 
the two layers at the 95% confidence level are indicated by *. 
T indicates the presence of a thermal bar. S indicates thermal 
stratification. 
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Table A-2, Comparison of surf ace/epilimnion total dissolved phosphorus with 

lower layer/ hypo limn ion total dissolved phosphorus. Average values 
which differ between the two layers at the 95% confidence level are 
indicated by *. T indicates the presence of a thermal bar. 
S indicates thermal stratification. 
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Table A-3. Comparison of surf ace/epilimnion dissolved ortho phosphorus with 

lower layer/ hypolimnion dissolved orthophosphorus . Average values 
which differ between the two layers at the 95% confidence level are 
indicated by *. T indicates the presence of a thermal bar. 
S indicates thermal stratification. 
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Table A-4. Comparison of surf ace/ epilimnioti dissolved reactive silica with 

lower layer/ hypolimnion dissolved reactive silica. Average values 
which differ between the two layers at the 95% confidence level 
are indicated by *. T indicates the presence of a thermal bar. 
S indicates thermal stratification. 
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Table A-5. Comparison of surf ace/ epilimnion dissolved nitrate + nitrite with 

lower layer/ hypolimnion dissolved nitrate + nitrite. Average values 
which differ between the two layers at the 95% confidence level are 
indicated by *. T indicates the presence of a thermal bar. 
S indicates thermal stratification. 
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Table A-6. Comparison of surf ace/ epilimnion dissolved ammonia with lower 

layer/ hypolimnion dissolved ammonia. Average values which differ 
between the two layers at the 95% confidence level are indicated 
by *. T indicates the presence of a thermal bar. S indicates 
thermal stratification. 
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Table A-7 . Comparison of surf ace/epilimnion specific conductivity with lower 
layer/ hypolimnion specific conductivity. Average values which 
differ between the two layers at the 95% confidence level are 
indicated by *. T indicates the presence of a thermal bar, 
S indicates thermal stratification. 
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Table A-8, 



Comparison of surf ace/ epilimnion alkalinity with lower layer/ hypo- 
limnion alkalinity. Average values which differ between the two 
layers at the 95% confidence level are indicated by *. T indicates 
the presence of a thermal bar. S indicates thermal stratification. 
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Table A-9. Comparison of surf ace/ epilimnion pH with lower layer/hypolimnion 
pH. Average values which differ between the two layers at the 
95% confidence level are indicated by *. T indicates the presence 
of a thermal bar. S indicates thermal stratification. 
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APPENDIX B 
STATISTICAL SUMMARIES FOR 1974 AND 1980 
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